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К 90-летию co дня рождения E. IIT. Молотова 


В конце декабря 2019 года ветерану нашего предпри- 
ятия Евгению Павловичу Молотову исполнилось 90 лет. 


Евгений Павлович родился в 1929 году в г. Харькове. 
После окончания в 1952 году Киевского политехнического 
института по специальности «Радиотехника» Евгений Моло- 
тов пришел на работу в НИИ-885 (сегодня — АО «Российские 
космические системы»). В настоящее время занимает долж- 
ность главного научного сотрудника. 

Работал в должности инженера, старшего, ведущего ин- 
женера, руководителя сектора, начальника лаборатории, Ha- 
чальника отдела, ведущего научного сотрудника. 
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но за его вклад в разработку, создание и испытания назем- 
ной аппаратуры ракетно-космической техники. Участник со- 
здания и испытаний систем радиоуправления ракетами P5, 
P7, P9, P36, P100. Руководил созданием наземной аппарату- 
ры первого космического спутникового геодезического ком- 
плекса «Сфера», комплекса контроля навигационного поля системы ГЛОНАСС, наземной аппаратуры 
системы «Коспас», радиотехнических комплексов управления пилотируемыми аппаратами «Квант-П», 
«Квант-СП», «Квант-Р», «Уран-Ц», комплексов управления космическими аппаратами для исследования 
дальнего космоса «Сатурн-МСД», «Квант-Д». 

Е. П. Молотов был техническим руководителем наземного комплекса управления при выполне- 
нии космических программ по исследованию Луны и дальнего космоса «Mapc-3-7», «Марс-7-16», 
«Венера-7-16», «Вега», «Фобос», «Прогноз», «Астрон», «Гранат», «Интербол». И это далеко не весь 
перечень. 

За большой вклад в развитие отечественной космической техники Е. П. Молотову присвоено зва- 
ние «Лауреат Государственной премии СССР», он награжден орденами Октябрьской Революции, «Знак 
Почета», юбилейной медалью «За доблестный труд. В ознаменование 100-летия со дня рождения 
В. И. Ленина», нагрудным знаком «Почетный радист» и другими наградами. 

Е. П. Молотов — доктор технических наук, профессор, академик Российской академии космонавти- 
ки им. К.Э. Циолковского, автор более 400 научных трудов, статей, монографий и изобретений. 

В 2004-2008 преподавал на кафедре «Радиосистемы управления и передачи информации» Москов- 
ского авиационного института. 

Является членом научно-технического и диссертационного советов АО «Российские космические 
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Фазоразностная радиотомография ионосферы 
на базе кластера из двух космических аппаратов 


О. И. Барабошкин, baraboshkin oi(Qspacecorp.ru 
АО «Российские космические системы», Москва, Российская Федерация 
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АО «Российские космические системы», Москва, Российская Федерация 
C. A. Бобровский, bobroosky. sa(Qspacecorp.ru 
АО «Российские космические системы», Москва, Российская Федерация 


Аннотация. Рассмотрен способ получения двумерных распределений электронной концентрации в ионосфере с помощью кла- 
стера из двух космических аппаратов фазоразностным методом по двухчастотным измерениям. Приведено описание алгоритма, 
позволяющего построить проекционный оператор для рассматриваемой геометрии расположения передатчика и приемника 
и решить систему линейных уравнений с учетом неравномерного заполнения области реконструкции трассами радиосигналов. 
Проведено моделирование задачи радиотомографии при наличии неоднородностей в слое максимума электронной концентра- 
ции. Показано, что предлагаемый способ позволяет восстанавливать двумерные распределения электронной концентрации, 


в том числе крупномасштабные структуры и локальные неоднородности. 


Ключевые слова: радиотомография, ионосфера, ПЭС, кластер, космический аппарат 


Phase Difference Radio Tomography ої the Ionosphere 
Based on Two Spacecrait Cluster 


О. I. Baraboshkin, baraboshkin, oi(Qspacecorp.ru 
Joint Stock Company "Russian Space Systems", Moscow, Russian Federation 


S. V. Trusov, Cand. Sci. (Engineering), ігиѕоо su(Ospacecorp.ru 
Joint Stock Company "Russian Space Systems", Moscow, Russian Federation 
S. A. Bobrovskiy, bobrovsky sa(Qspacecorp.ru 
Joint Stock Company "Russian Space Systems", Moscow, Russian Federation 


Abstract. The paper studies the method to obtain two-dimensional electron density distributions in the ionosphere by means 
of a two spacecraft cluster by a phase difference method by double-Írequency measurements. An algorithm is described, which 
allows constructing a projection operator for the considered geometry of location of the transmitter and receiver and solving 
a linear equation system taking into account non-uniform filling of the reconstruction area with radio signal paths. The task 
of radiotomography was simulated in the presence of inhomogeneities in the layer of maximum electron concentration. It is shown 
that the proposed method allows restoring two-dimensional distributions of electron concentration including large-scale structures 


and local inhomogeneities. 


Keywords: radio tomography, ionosphere, total electron content, cluster, spacecraft 
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Введение 


Метод радиотомографии (РТ) ионосферы поз- 
воляет получать двумерные распределения элек- 
тронной концентрации в ионосфере. Широкий про- 
странственный охват, высокая точность, небольшое 
время накопления измерений позволяют обнаружи- 
вать и исследовать данным методом быстропроте- 
кающие процессы в ионосфере, которые могут по- 
влиять на работу космических систем связи, на- 
вигации и др. Наибольшее распространение по- 
лучил метод РТ ионосферы сигналами «спутник- 
Земля» [1,2], при котором использовались назем- 
ные приемные станции и сигналы навигационных 
спутников первого поколения. Однако ввиду окон- 
чания эксплуатации указанных космических аппа- 
ратов (КА) для применения данного подхода в бу- 
дущем требуется запуск новых спутников. Поэтому 
в настоящее время перспективным представляется 
переход к РТ ионосферы кластером КА, которая, 
в отличие от наземной, позволяет осуществлять 
глобальный мониторинг ионосферы. 

Спутниковая РТ ионосферы кластером КА бы- 
ла впервые предложена и рассмотрена в работах 
[3,4], в которых определена оптимальная конфигу- 
рация кластера, состоящего из четырех КА, а так- 
же предложено определять полное электронное со- 
держание на основе фазовых измерений, получа- 
емых четырехчастотным методом [5]; представлен 
предполагаемый облик аппаратуры для КА нано- 
класса. При этом не рассмотрены вопросы форми- 
рования проекционного оператора в случае радио- 
трасс «спутник-спутник», выбора метода решения 
систем линейных уравнений (СЛУ), не проанали- 
зирована возможность реализации схемы двухча- 
стотных измерений, применявшегося ранее для РТ 
ионосферы сигналами «спутник-Земля» [2]. 

В работе рассматривается способ РТ ионосфе- 
ры кластером КА, предполагающий использование 
двух КА в комбинации с двухчастотным фазораз- 
ностным методом. Способ обладает меньшей слож- 
ностью реализации в связи с тем, что производство 
и запуск двух КА дешевле, чем четырех, а энер- 
гопотребление двухчастотных передатчика и при- 
емника существенно ниже четырехчастотных, что 
позволяет разместить их на платформах, аналогич- 
ных CubeSat 1U или ТНС-1 [6]. 


Схема эксперимента и построение 
проекционного оператора 


Введем ряд параметров, характеризующих гео- 
метрию системы регистрации в полярной системе 
координат (CK) (рис. 1): (œj, №) — координаты КА 
передатчика, (o5, Ao) — координаты КА-приемника, 
О — центр Земли, ОО’ — перпендикуляр, опущен- 
ный из точки О на траекторию луча. Траекторию 
луча в полярной СК можно представить в виде 

I 
"n cos(a — В)’ 
где P — длина перпендикуляра ОО’, В — угол nep- 
пендикуляра к траектории луча, отсчитываемый от 
начала координат. Расстояние l от точки О’ до точ- 
ки на луче с координатами (о, ћ) 


l= Ptg(o — В). 


(1) 


(2) 


Рис. 1. Схема эксперимента по лучевой PT кластером 
из двух КА 


Для элемента длины луча справедлива формула 


Al - [xd 
cos" (а — B) 


тогда линейный интеграл набега фазы сигнала ф, 
обусловленного влиянием ионосферы, равен 


(3) 


Q9 Q9 
P = Are | ав. da = Ar, | (a) Е (а, h) da, 
eq eq 


(4) 


где F(o,h) — двумерная функция распределения 
плотности электронной концентрации в плоскости 
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радиопросвечивания, À — длина волны сигнала, 
т. — классический радиус электрона. 

Заменим непрерывную функцию F(o, h) на ее 
кусочно-планарную аппроксимацию на сетке Fyn 
с дискретом по углу Да и высоте Ah (рис. 2). 


Рис. 2. Кусочно-планарная аппроксимация двумерной 
функции на сетке реконструкции 


В верхнем треугольном элементе 


Eo я 
Plah) = Fapa (В — Һр) 
Fmtinti Етті (5) 
+ Am (o — ат). 


В нижнем треугольном элементе 


Е т E 


F(a, h) = Fugiat m Ah mib (h — h) + 
do 1% — Dn 
pma (a anyi). (6) 
В верхнем 


Jer = dob us uai ЕЯ V POET. = Pon 


t JU ciiin = Frangi) 


B нижнем 


MODE i Jo skin e JA тъ > Вет 
Е do E niin > T n: (8) 


Bases Jod dads dba 


а — ЭТО интегралы: 


Jg em um po do, (9) 


1 


n дя | 20601009) -h.u]de, (10) 


Jh = ag | OPa) — в. do 


Ja = да | в — o] do 


1 
Л = да [tllo — ossi] do. 


(11) 
(12) 


(13) 


После интегрирования по всем ячейкам сет- 
ки реконструкции, в которые попадает траектория 
7-го луча, в элемент матрицы Lyn n заносятся COOT- 
ветствующие Fm n коэффициенты и T. д.: 


Loin — —Jn > Ja 
T S = Jo + Jy Е Jo; 
афа — Jo > Jf, + ds 


Lonain = Л s Jy: 


(14) 
(15) 
(16) 
(17) 


Тогда интеграл (4) можно представить B виде 
суммы 
Q — м, Ly Fy, (18) 
M 

где произведена перенумерация (m, n) — М. Таким 
образом, задача нахождения неизвестных значений 
электронной концентрации в узлах сетки рекон- 
струкции Ем сведена к решению системы линей- 

ных уравнений вида 
Pi 


pv" LimFm- 


(19) 
B фазоразностной постановке задачи томогра- 
фии измеряется не фаза сигнала p, а ee производ- 
ная 14/11, поэтому систему (19) необходимо npo- 
дифференцировать по времени. Это возможно кор- 
ректно сделать в силу того, что применялась ку- 
сочно-планарная аппроксимация функции F, cje- 
довательно, производные р; не имеют разрывов: 


1 de; Гум — Гум 
Ac. s Ем = Р;мЕм, (20) 
где Lim — матрица Lim построенная для Tpa- 


екторий лучей в момент времени f, — At, Dim — 
проекционный оператор. 

Интегрирование (9)-(13) можно сильно упро- 
стить путем введения кусочно-линейной аппрокси- 
мации траектории луча. Пусть внутри 1-й ячейки 
сетки реконструкции траектория луча задается фор- 
мулами (рис. 3) 


h = Ка + b; Qil S Q € Os, (21) 
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(22) 
(23) 


К; = (hi — йа) (а — 0); 
b; = hig — Кодо, 


где (oj, hj) и (Qio Rio) — точки пересечения тра- 
ектории луча с 1-й ячейкой сетки реконструкции, 
найденные по формуле (1) для граничных значе- 
ний а и h в ячейке. 


Рис. 3. Аппроксимация луча внутри ячейки сетки рекон- 
струкции 


Диагональ ячейки сетки реконструкции зада- 
ется формулой 


(24) 
где (o; h;) — координаты нижнего левого угла 


1-й ячейки. Тогда координаты пересечения диаго- 
нали ячейки и траектории луча 


Ah 
" h; — b; — Tali (25) 
Е ГЗ Ah 
i^ Aa 


Если а; € Qie < ао, то траектория луча nepe- 
секает и нижний, и верхний треугольные элементы 
1-й ячейки. Координаты Qil, о, могут быть rpa- 
ничными значениями как для верхнего, так и для 
нижнего треугольного элемента, поэтому для Qil < 
X Q X Qic H Oc < Q X Qo необходимо определить, 
в какой треугольный элемент попадает траектория 
луча. Пусть в первом случае траектория попадает 
в верхний элемент, а во втором — в нижний, тогда 
интегралы (9)-(13) вычисляются по формулам: 


Јо = Ptg(o.— 8) |4 = Ptg(oj, — 8) - Рав (а; — В), 
(26) 
Ј = Ptg(o — В) (27) 


Л = a (klin [cos(a — 8) + (a — В) tgla — 8) 
(KB + b; — h; — Ah)tg(a 8)) |2 | (28) 

J, = 2. (вова — + (a — B) te(o — 8)+ 
T(k8 + b; — А) ів(а — 8)) | T (29) 

Ja = A (сова — )| (o. Bua - B 


- 8)) | (30) 
а в) un — В)+ 
- 8)) |=. (31) 


Qu 

Если Qj, < Qj; ИЛИ Qe > ао, то траектория лу- 
ча т. один треугольный элемент 1-й ячей- 
ки. Интегралы (26)-(31) рассчитываются на отрез- 
ке а, < а < а, только для треугольного элемента, 
в который попадает траектория луча. 

Так как расчет интегралов выполняется анали- 
тически, а не численно, процедура построения про- 
екционного оператора эффективна, скорость рас- 
чета не влияет на оперативность решения зада- 
чи томографии. К тому же при радиопросвечива- 
нии спутниками на круговых орбитах относитель- 
ная конфигурация источников и приемников сиг- 
нала неизменна, и проекционный оператор может 
быть построен единожды и не перестраиваться ли- 
бо изредка корректироваться. 


de A (In | eos(a — 8 | 


8) 
(8 — o4) (а 
1+ 
tg( 


— oj — 0)tg 


Решение системы линейных 
уравнений 


В работе [2] приведены результаты анализа 
различных алгоритмов решения систем линейных 
уравнений (СЛУ) применительно к задачам РТ — 
ART, SIRT, MART, DART и др. Эксперименты по 
наземной РТ ионосферы показали, что однозначно 
выбрать алгоритм решения СЛУ не предоставляется 
возможным, так как в разных условиях (метод ре- 
конструкции, величина ошибок в исходных данных, 
наличие определенных структур в ионосфере и др.) 
лучше работают разные алгоритмы. Преимущества 
и недостатки алгоритмов в целом распространяются 
и на случай спутниковой РТ ионосферы. 

В настоящей работе для решения СЛУ (20) был 
выбран алгоритм DART. Однако ввиду неравномер- 
ного заполнения области реконструкции трассами 
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радиосигналов алгоритм был дополнен усредняю- 
щим коэффициентом 1/P;: 


т; — (Ах 
АН а. BI 2| (» 

J J Tu кА?) 
где X — вектор- столбец неизвестных Tj ze E 
номер итерации, E; — число ненулевых ай-е ер 


B j-M столбце матрицы А,., А? — 1-я строка матри- 
цы А;,, т, (вектор-строка т) — данные измерений. 

Моделирование решения фазоразностной зада- 
чи РТ показало, что в данной геометрии распо- 
ложений КА при восстановлении небольших ло- 
кализованных неоднородностей происходит иска- 
жение реконструкции, проявляющееся в виде уве- 
личения электронной концентрации в областях 
пространства, через которые проходят лучи, пере- 
секающие неоднородность. Чтобы минимизировать 
данный эффект, а каждой итерации после нахож- 
дения поправки ps = = gt a находятся ee макси- 
мальное и минимальное значения, и если поправка 
в 7-ю ячейку удовлетворяет HepaBeHcTBy 


о хо) - вв), 
min(pj)-F ———,— << 
MED 
& max(p;) Dd (33) 


TO p; обнуляется. Такая операция позволяет на пер- 
вых итерациях вносить поправки только в обла- 
сти неоднородности, где ошибка электронной кон- 
центрации максимальна, что обеспечивает повы- 
шение приоритета восстановления неоднородностей 
при их наличии и приводит к улучшению качества 
реконструкции. 

Для регуляризации решения СЛУ в задаче РТ 
ионосферы в конце каждой итерации перед сло- 
жением тк с поправкой p, целесообразно приме- 
нять сглаживающий фильтр, что использовалось 
в ряде работ [2,7]. Данную процедуру можно ин- 
терпретировать как введение в решение априорной 
информации о гладкости распределения электрон- 
ной концентрации. В настоящей работе применял- 
ся КИХ-фильтр с коэффициентами [0,0236, 0,0927, 
0,2324, 0,3026, 0,2324, 0,0927, 0,0236] по обеим ко- 
ординатам. 


На малых высотах значение электронной кон- 
центрации равно нулю, поэтому поправка для ячеек 
сетки реконструкции, соответствующих этим высо- 
там, не вносится. Контроль сходимости итерацион- 
ного процесса производится путем вычисления нор- 
мы невязки || Az — m|| в пространстве {2. Если e(l?) 
на следующей итерации больше, чем на предыду- 
щей, то итерационный процесс прерывается: 


(34) 


Результаты моделирования 
и анализ результатов 


Моделирование задачи РТ ионосферы класте- 
ром КА состояло из следующей последовательно- 
сти функций. 

1. С использованием модели SGPA [8] и файла 
TLE спутника с подходящими параметрами орби- 
ты производился расчет положений двух спутников 
в рассматриваемый отрезок времени и радиотрасс 
сигналов без учета их преломления в ионосфере. 

2. Исходя из крайних положений КА в рассмат- 
риваемый отрезок времени и размера ячейки сетки 
реконструкции формировалась двумерная сетка то- 
мографической реконструкции. 

3. С использованием модели ионосферы 
NeQuick [9] формировалось начальное распреде- 
ление электронной концентрации на сетке рекон- 
струкции. При решении модельной задачи это рас- 
пределение использовалось в качестве начального 
приближения. 

4. Из начального распределения формирова- 
лось модельное распределение путем внесения ис- 
кажений двух типов: 

а) случайное гладкое изменение значений 
электронной концентрации в пределах 20% вдоль 
орбиты КА, имитирующее крупномасштабные воз- 
мущения в ионосфере в диапазоне погрешностей 
модели NeQuick; 

6) локальная неоднородность с характерным 
размером 500 x 50 км и амплитудой, равной фо- 
новому значению в данной точке. 
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5. IIo модельному распределению электронной 
концентрации для каждой радиотрассы по форму- 
ле (4) методом трапеций численно рассчитывались 
фазовая задержка сигнала и ее производная M. 

6. Для положений KA, полученных из п. 1 
и сетки реконструкции из п. 2, рассчитывался про- 
екционный оператор А,.. 

7. Итерационно решалась СЛУ DE Ах; = 
=m; 

8. Производилась оценка ошибки реконструк- 
ции путем вычисления норм Ô в пространствах І? 
и [°, которые характеризуют отклонение рекон- 


струируемой функции F от исходной функции F: 


(35) 


(36) 


Так как метод фазоразностный и фактиче- 
ски восстанавливается разность между начальным 
и модельным распределением, показательна нор- 
ма д для этой разности — она характеризует каче- 
ство восстановления добавленных к модели неод- 
нородностей. 

При прохождении радиосигнала через ионо- 
сферу сигнал преломляется и при низкой частоте 
сигнала может отразиться. Отклонение от прямо- 
линейного распространения сигнала максимально 
в области максимума ионосферы, особенно при ма- 
лых углах падения [10]. Поэтому для стабильно- 
сти приема сигнала целесообразно расположить КА 
так, чтобы радиотрассы проходили ближе к поверх- 
ности Земли, тогда угол падения на высоте мак- 
симума электронной концентрации будет больше. 
При этом при достаточной высоте орбиты радио- 
трассами будет охвачен весь необходимый для мо- 
ниторинга ионосферы диапазон высот. 

Расстояние между КА при фиксированном ми- 
нимальном расстоянии между радиотрассой и по- 
верхностью Земли т зависит от высоты орбиты, 
и чем высота меньше, тем быстрее формируется на- 
бор данных для каждой ячейки сетки реконструк- 
ции. Для высоты орбиты 500 км и минимального 


расстояния т = 60 км время накопления данных 
для каждой ячейки менее 11 мин, что позволяет 
считать ионосферу квазистационарной. 

Моделирование проводилось для двух КА на 
круговых приполярных орбитах высотой около 
500 км, период времени проведения измерений — 
30 мин, размер ячейки сетки реконструкции — 25 х 
х 25 км, ширина области реконструкции — около 
17 тыс. км (рис. 4). Исходное (модельное) распре- 
деление электронной концентрации и результат его 
томографической реконструкции представлены на 
рис. 5, 7. 


КА! 


KA2 0430 MHH 


KA 1 104-30 Ma 


КА? d 


Рис. 4. Область реконструкции и положения КА B на- 
чальный и конечный моменты времени 


Результат моделирования при наличии круп- 
номасштабных возмущений показан на рис. 5, 6. 
Ошибки реконструкции для разности модели ионо- 
сферы и начального приближения 6(1?) = 0,5, 
5(1?) = 0,60. Результат моделирования при Ha- 
личии дополнительно локальной неоднородности 
представлен на рис. 7, 8, ошибки реконструкции 
в этом случае 6(1?) = 0,55, 6(12°) = 0,67. 

Моделирование решения задачи РТ показало, 
что предложенный способ позволяет восстанавли- 
вать как крупномасштабные возмущения в ионо- 
сфере, так и локальные неоднородности. Абсолют- 
ные значения амплитуд небольших неоднородно- 
стей на реконструкциях в среднем составляют около 
половины от реальных значений, при этом круп- 
ные структуры восстанавливаются практически пол- 
ностью. Наибольшие искажения в реконструкциях 
возникают на краях, так как в этих областях ши- 
риной, равной расстоянию между КА, недостаточно 
проекций. Поэтому расчет необходимо производить 
для реконструкций шириной существенно больше 
двух расстояний между КА. При этом время меж- 
ду последовательными расчетами реконструкций 
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Рис. 5. Модель ионосферы (вверху) и реконструкция модели методом фазоразностной РТ (внизу) при наличии 
крупномасштабных возмущений 
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Рис. 6. Разность между моделью ионосферы и начальным приближением (вверху) и разность между реконструк- 


цией модели методом фазоразностной РТ и начальным приближением (внизу) при наличии крупномасштабных 
возмущений 
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Рис. 7. Модель ионосферы (вверху) и реконструкция модели методом фазоразностной РТ (внизу) при наличии 
крупномасштабных возмущений и локальной неоднородности 
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Рис. 8. Разность между моделью ионосферы и начальным приближением (вверху) и разность между реконструк- 


цией модели методом фазоразностной РТ и начальным приближением (внизу) при наличии крупномасштабных 
возмущений и локальной неоднородности 
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может быть небольшим, так как измерения могут 
использоваться в нескольких расчетах, тогда опе- 
ративность получения данных после пролета двух 
КА над исследуемой областью в случае наличия 
связи с Землей или использовании данных непо- 
средственно на КА будет порядка минуты. 


Выводы 


Рассмотрен способ получения двумерных рас- 
пределений электронной концентрации в ионосфере 
с помощью кластера из двух космических аппара- 
тов фазоразностным методом по двухчастотным из- 
мерениям. Такое сочетание количества КА в кла- 
стере и выбранного метода является наиболее про- 
стым способом получения двумерных распределе- 
ний электронной концентрации с точки зрения прак- 
тической реализации. Рассмотренный алгоритм фа- 
зоразностной РТ адаптирован к геометрии РТ ионо- 
сферы кластером КА, позволяет с минимальными 
вычислительными затратами построить проекцион- 
ный оператор и решить систему линейных уравне- 
ний с учетом неравномерного заполнения области 
реконструкции трассами радиосигналов. Проведен- 
ное моделирование показало, что предлагаемый спо- 
соб позволяет восстанавливать двумерные распреде- 
ления электронной концентрации, в том числе круп- 
номасштабные структуры и локальные неоднород- 
ности. Достоинством рассмотренной методики так- 
же является возможность применения частот 150 
и 400 МГц, которые в данный момент используют- 
ся в целях мониторинга ионосферы. В этом случае 
можно использовать сигналы КА-передатчика также 
для наземной томографии ионосферы на уже развер- 
нутых сетях томографических приемников. 
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Аннотация. Обработка координатно-трассовой информации о мобильных объектах от нескольких разнотипных источников 
на центрах обработки информации требует высокой достоверности формируемых результатов. В статье предлагаются методы 
расчета доверительного интервала локации объекта, методы расчета экстраполированного доверительного интервала локации 
объекта, зависящие от маневра объекта, и методы вычислений параметров и коэффициентов алгоритмов обобщения дан- 
ных от нескольких источников информации. Также предлагаются методы вычисления времени хранения трасс, потерянных 
источниками сопровождения. 

Разработанные методы позволяют гибко использовать данные алгоритмы обобщения, расширяя возможность их примене- 


ния в условиях разрывного поля видимости источников информации, существенной разноточности измерений и неопределен- 
ности гипотез движения объектов. 
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Обработка координатно-трассовой информа- 
ции, предназначенная для оценки текущего состо- 
яния объектовой обстановки и формирования в ре- 
жиме реального времени описания объектовой об- 
становки в виде единого потока обобщенных трасс 
объектов, осуществляемая на основе информации 
от нескольких источников, является задачей тре- 
тичной обработки информации. Третичная обра- 
ботка информации применяется на центрах обра- 
ботки информации, например на центрах обработ- 
ки радиолокационной информации, принимающих 
информацию от замкнутых на них радиолокаци- 
онных станций. Она реализуется путем объедине- 
ния потоков трасс объектов, поступающих от мно- 
жества различных источников информации [1,2], 
с исключением дублирования информации разных 
источников и повышением достоверности сформи- 
рованного описания объектовой обстановки. 

Для решения данной задачи привлекаются раз- 
нотипные источники информации, работа которых 
основана на различных физических принципах. Ин- 
формация, поступающая от них, различается: 

— темпом поступления сообщений; 

— способом наблюдения объекта — эхолока- 
ция, пеленгация, вторичная локация; 

— априорной и текущей точностью; 

— длительностью сопровождения объекта. 


Третичная обработка решает следующие задачи: 


— отождествление обобщенной трассы объекта 
с информацией, приходящей в сообщениях по тому 
же объекту от нескольких разнородных рассредо- 
точенных источников информации; 

— формирование новой обобщенной трассы; 

— обновление координат и вектора скорости 
(оценки положения трассы объекта в пространстве) 
обобщенной трассы; 

— сброс обобщенной трассы в случае прекра- 
щения поступления информации по объекту от ис- 
точников информации. 

Цель данной работы — разработка методов 
третичной обработки информации с использова- 
нием интервальных оценок положения объекта, 
позволяющих повысить достоверность формируе- 
мого описания объектовой обстановки в условиях 
разрывного поля видимости источников информа- 
ции, их существенной разноточности, а также при 
неопределенности гипотез движения объектов. 


Интервальная оценка положения 
трассы объекта 


В задачах оценки параметров при небольшом 
объеме выборки, к которым относится задача оцен- 
ки положения объекта в пространстве по данным 
локационных измерителей, более предпочтительна 
интервальная оценка [5-7]. В отличие от точечной 
оценки интервальная оценка (доверительный ин- 
тервал 4718.) представляется интервалом 3Haue- 
ний, внутри которого с задаваемой доверительной 
вероятностью (5) находится истинное значение (х) 
оцениваемого параметра: 


P (|£ —m,| < trusta) = В, 


где m, — математическое ожидание величины T. 
Из неравенства Чебышева 


P (|£ — m,| 2 k x o) S 


2 
t| ^7 


следует, что доверительный интервал приблизи- 
тельно пропорционален среднеквадратическому от- 
клонению (о). 

Например, для нормального закона это точная 


пропорция: 
trust 
o crd = f, 
(58) =0 


где Ф(2) — функция Лапласа. 
Для 8 = 0,9 


Вена = 1,7 +0. 


Неопределенность положения экстраполирован- 
ной обобщенной трассы можно оценить доверитель- 
ным интервалом (truste а), пересчитанным на Bpe- 
мя экстраполяции. Д 


Экстраполяция доверительного 
интервала 


В случае, когда объект движется прямолиней- 
но и равномерно 


т=з + У. 


РАКЕТНО-КОСМИЧЕСКОЕ ПРИБОРОСТРОЕНИЕ И ИНФОРМАЦИОННЫЕ СИСТЕМЫ т. 7 вып. 1 2020 


МЕТОДЫ ПРИМЕНЕНИЯ ДОВЕРИТЕЛЬНОГО ИНТЕРВАЛА ЛОКАЦИИ ОБЪЕКТА 15 


и гипотеза движения совпадает с движением объек- 
та, экстраполяция производится по формуле 


Ti — Ti] tg tg 
senti Hg s (Lex) -na 8). 
где т, и 2; — измеренные координаты на 1-й 


и (i — 1)-й моменты времени; 

At — время, за которое производится подсчет 
скорости (темп измерения источника); 

ig — время, на которое производится экстра- 
поляция. 

Дисперсия экстраполированного значения ко- 
ординаты равна сумме дисперсий независимых H3- 
мерений: 


D = D, (i18) «D . (=) 
BG ^t 1—1 At 


Дисперсии 1-го и (i— 1)-ro измерений равны 


В случае если trust, — доверительный интер- 
вал обобщенной трассы на момент последнего ее 
обновления, то экстраполированный доверительный 
интервал (ёгиѕё. ..,) подсчитывается по формуле 


(тазбЕ_сга = газа ' 


В случае, когда объект движется прямолиней- 
но и равноускорено: 


2 


экстраполяция должна производиться по формуле 


— И 
XE др NE "+0. = 
— ig їв ЇЕ 
= (S n (22 qm 


Однако ускорение (U) не рассчитывается. 
Предполагая, что ускорение, с которым может дви- 
гаться объект, является случайной величиной, рас- 
пределенной по нормальному закону с дисперси- 
ей Ду и нулевым математическим ожиданием, дис- 
персия экстраполированного значения координаты 
будет 


igY* 
Отсюда экстраполированный доверительный 
интервал (truste ..,) подсчитывается по формуле: 


2 
traste erd = (trasta: (1 +2. E +2. (25) ) + 


1/2 
+ trast? . by р . 
о 


На рис. 1 приведены результаты расчета зна- 
чения экстраполированного доверительного интер- 
вала в метрах (ось ординат) в зависимости от вре- 
мени (ось абсцисс) в секундах. Графики построены 
при параметрах, указанных в табл. 1. 


(1) 


Таблица 1. Параметры зависимостей экстраполиро- 
ванного доверительного интервала от времени 


Номер графика ) | trusty (м/с?) 


4 |с |м 


Траектория движения объекта состоит из от- 
резков, когда объект не маневрирует (прямолиней- 
ное равномерное движение), и отрезков, когда объ- 
ект совершает маневр. Пусть Pyy — вероятность 
того, что объект находится на отрезке траектории 
маневра. Тогда плотность распределения случайной 
величины U будет определяться выражением для 
смеси вероятностных распределений [3] 


F(U) = (1 — Pu): foU) + Ру: Рм(0), 


где fg(U) — плотность распределения случайной 
величины U на отрезке траектории без маневра 
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Рис. 1. 


с дисперсией Dyo, fy(U) — плотность распределе- 
ния случайной величины U на отрезке траектории 
с маневром Дум. 

Предполагая, что ускорение, с которым может 
двигаться объект, как при маневре, так и без ма- 
невра является случайной величиной, распределен- 
ной по нормальному закону с дисперсией DU и ну- 
левым математическим ожиданием, получим 


Dy = |001 Pu) fU) + Ри fu (U))dU = 


= va = Ри): 0) + | U* Pu fuU) 40, 
Dy = (1 - Py): Dyo + Py: Бим. 


Подставляя в (2) данные результаты, получим 


2 


2. 1/2 
n 


Или 
‚в (2) 
(1+2 E 2: (8) } + 


24 1/2 
+ trast?. . ів | 
U 2 d 


traStg erd = (ies : 


где 


trasta = \/(1— Py) - газо + Py, - газе y. 


На рис. 2 приведены результаты расчета значе- 
ния экстраполированного доверительного интерва- 
ла в метрах (ось ординат) в зависимости от време- 
ни (ось абсцисс) в секундах в случае обнаружения 
маневра (тапеуг = 1) и в случае, когда маневр не 
обнаружен (manevr = 0). 

Расчеты сделаны при следующих значениях 
параметров: trast = 250, trastpy = 0,2 · 9, 
trastp м =3:9. 


crd 


Использование доверительных 
интервалов в алгоритме 
отождествления 


Идентификация информации источника с o606- 
щенной информацией решается методами отож- 
дествления на основе сравнения информации, при- 
шедшей от источника (2.5, Утза), C информаци- 
ей обобщенной трассы (LyenTrcE> YgenTreE)» ЭКстра- 
полированной на тот же момент времени. В слу- 
чае, когда гипотеза движения объекта совпада- 
ет с реальным его перемещением, отклонение 
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информации источника от экстраполированной HH- 
формации обобщенной трассы зависит только от 
случайных ошибок измерения источниками ин- 
формации. В случае, когда гипотеза экстраполя- 
ции не соответствует действительному движению 
объекта, к случайной ошибке измерения прибав- 
ляется ошибка экстраполяции, зависящая от дина- 
мических характеристик объекта. Это обстоятель- 
ство чрезвычайно важно учитывать при отождеств- 
лении высокоскоростных, маневренных объектов. 
Поэтому достоверное решение задачи отождеств- 
ления возможно лишь в пределах ограниченно- 
го времени экстраполяции (fg), рассчитываемо- 
го на основе случайной погрешности измерений, 
определяемого маневра, скорости и типа объекта 
и др. 

При отождествлении трассы источника, сооб- 
щение о котором пришло впервые, в качестве кан- 
дидатов на отождествление могут быть отобраны те 
обобщенные трассы, доверительный интервал кото- 
рых, экстраполированный на время пришедшего со- 
общения, имеет пересечение с доверительным ин- 
тервалом отождествляемой трассы источника. До- 
верительный интервал отождествляемой трассы ис- 
точника определяется источником в зависимости от 
применяемых им методов и приходит в сообщении. 

Доверительный интервал отбираемых обобщен- 
ных трасс, экстраполированный на время пришед- 
шего сообщения, рассчитывается по формуле (2). 


Использование доверительных 
интервалов в алгоритме обновления 
координат 


Фильтрация 


Пусть 2.3, Ymsg — Координаты объекта, при- 
шедшие в сообщении от источника информа- 
ЦИИ, ТрептесЕ» UgenTrcE координаты обобщенной 
трассы, экстраполированные на время сообщения 
источника. 

Обновление координат обобщенной трассы 


и параметров движения Va 


T oenTrc , UsenTrc genTrc' 


V YgenTre производится методами фильтрации по 


формулам: 
TgenTre = T'genTrcE фах (аы > TgenTrcE); 


YgenTre = YgenTrcE tax (Ymse ni YgenTrcE); 


= Tmsg — TgenTrcE , 3 
V ZgenTre ad V sitios Ex At ту d ( 
Umsg EE UgenTrcE 
V YgenTre = V YgenTreS LE At + Ty ) 


где At — разница между временем последнего об- 
новления координат обобщенной трассы и време- 
нем сообщения источника; 

œa — коэффициент фильтрации координат 
(значение зависит от признака наличия маневра 
объекта); 


РАКЕТНО-КОСМИЧЕСКОЕ ПРИБОРОСТРОЕНИЕ И ИНФОРМАЦИОННЫЕ СИСТЕМЫ т. 7 вып. 1 2020 


18 


Ty коэффициент фильтрации скорости 
(значение зависит от признака наличия маневра 
объекта); 

V senTreS? V V een Tros параметры движения 
обобщенной трассы при предыдущем обновлении. 


Расчет доверительного интервала 
при фильтрации обобщенной трассы 


Координаты обобщенной трассы (например, 
VoenTrc) можно представить (см. формулы (3)) B ви- 
де суммы двух случайных величин Zope И 2.5! 


TgenTrc — TgenTrcE x (1 а) + Tmsg x a. 


Тогда дисперсия случайных величин будет вы- 
числяться по формуле 


= 2 2 
D genTre = DgenTrcE x (1 — а) + Рр ха. 


А доверительный интервал — 


Z 2 2 2 2 
tTaStgenTre = 4/tTaStzenTreE X (1 — 0)? + тавр x o. 


(4) 


Расчет оптимального значения 
коэффициента фильтрации 


Определим значение коэффициента фильтра- 
ции координат (о), обеспечивающего наименьшее 


значение trast,.,,,, рассчитанное по формуле (4). 


C. IO. САМСОН, А.В. НОВИКОВ, A. B. СОРОКИН 


/ аси. 


Для этого найдем производную (бав c.m.) 


d(trastsenTrc)/da = 
=d ( trast?, mg х (1 а)? + trast2.. х a?) /da = 


2 
trast hsg X 


y trast 


Приравняем нулю полученное выражение для про- 
изводной и найдем значение а, обеспечивающее 
данное равенство: 


а = trastz. my x (1— a) 


(=a)? 


2 ИИ 
msg X Q tr aSUsenTrcE 


trast?.. ха- trastéenTreE x (l-a) _ 


(1-а)? 


0. 
2 


9 
X а + trastsenTrcE 


«даф 
табби 


Таким образом, наименьшее значение 


trast обеспечивается при значении а: 


genTrc 


2 
genTrcE 


Е + trast 


trast 


2 
genTrc 


а = 


Е trast MN 


msg 

На рис. 3 приведены результаты расчета значе- 
ния коэффициента фильтрации координат (ось орди- 
нат) в зависимости от времени (ось абсцисс) в се- 
кундах в случае обнаружения маневра (manevr = 1) 
и в случае, когда маневр не обнаружен (тапеуг = 
= 0). Расчеты сделаны при следующих значени- 
ях параметров: trast g = 250, trasto = 0,2 · 9, 
trastyy =3 + 9. 
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Рис. 3. Результаты расчета значения коэффициента фильтрации координат 
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Расчет времени хранения трасс, 
потерянных источниками 
сопровождения, на основе 
доверительного интервала 


Обобщенная трасса, потерянная источниками 
сопровождения, должна храниться в системе обра- 
ботки информации до тех пор, пока значение экс- 
траполированного доверительного интервала по ко- 
ординатам обобщенной трассы не превысит порог, 
при котором отождествление с вновь поступившей 
трассой источника станет недостоверным. 

В качестве меры неопределенности информа- 
ции используется понятие информационной энтро- 
пии. Тогда, допуская, что оцениваемые величи- 
ны имеют непрерывные значения и распределе- 
ны по нормальному закону (характеризующемуся 
наибольшей энтропией), получаем из [4] следую- 
щие выражения: для энтропии экстраполированной 
обобщенной трассы 


1 In(216z., m6) 
Нрентев — 9 In(2) 


и энтропии входного сообщения о трассе источни- 
ка 
1 In(2702 се) 


s: In(2) 


msg ^ 9 


Количественная мера, которую может внести 
информация обобщенной трассы в процесс отож- 


дествления с вновь поступившим входным сообще- 
нием о трассе, будет определяться по формуле 


T= H msg = H рент = 
и Те) 1 In(2oz,. r6) = 


2 bhO 2 (9) 


1 1 ( О msg ) 
= n м 
In(2) О genTrcE 


Отношение дисперсий можно заменить отно- 
шением доверительных интервалов. Тогда 


1 trast msg 
d In(2) D (= ) | 


genTrcE 


Учитывая значение экстраполированного дове- 
рительного интервала (1), получим: 


1 
= n3 In (traste (ваё 
24 1/2 
SEU 2. (th 
к (1+2 ALT? (£ ини. (5 А 
(5) 


На рис. 4 приведена зависимость количествен- 
ной меры информации (— Г), которую может внести 
информация обобщенной трассы в процесс отож- 
дествления с вновь поступившим входным сообще- 
нием о трассе, от времени ее необновления (экстра- 
поляции), при значениях trust „g = 500, At = 10 с 
H trust = 0. 


0 
0 5 10 


15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 


Рис. 4. Зависимость количественной меры информации от времени ee необновления 
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Задавая критическое значение количественной 
меры I = Іо, можно получить значение порога по 
доверительному интервалу trusto, при достижении 
которого обобщенная трасса не вносит более досто- 
верной информации и сбрасывается с сопровожде- 


НИЯ. 
trast 
: In [ — 5 |. 
In(2) trastc 


Таким образом, для каждой обобщенной трас- 
сы, потерянной источниками сопровождения, мож- 
но рассчитать время ее хранения fq, после исте- 
чения которого данная трасса сбрасывается. 


1с — (6) 


На рис. 5 приведены результаты табличного 
расчета значения времени хранения (ось ординат), 
соответствующего достижению доверительного ин- 
тервала фтазёс = 20 км, в зависимости OT дове- 
рительного интервала по ускорению (ось абсцисс) 
и доверительного интервала по координатам обоб- 
щенной трассы на момент последнего ее обновле- 
ния (параметр графиков). 

На рис. 6 приведены результаты таблично- 
го расчета значения времени хранения, соответ- 
ствующего достижению доверительного интервала 
trasto = 40 км, в зависимости от доверительного 
интервала по ускорению и доверительного интерва- 
ла по координатам обобщенной трассы на момент 
последнего ее обновления. 


На рис. 5 и рис. 6 графики построены при зна- 


чениях параметров доверительного интервала, ука- 
занных в табл. 2. 


Таблица 2. Доверительный интервал по координатам 
обобщенной трассы как параметр зависимостей времени 
хранения от доверительного интервала по ускорению 


3 
4 
e 


250 
500 
750 
1000 
1250 


w 


Значение времени хранения вычисляется из 
уравнения 
[с == I. 


Учитывая выражения (5) и (6), получаем 


1 trast msg 
— ln — к 
In(2) trastc 
1 2. . Ёноа 

= (9) In СЯ (1 +2 а 


9 9 231/2 
+2. (=) + trast - thola 
At U 2 | 


Или 


2 
trastç = (is (1 +2. "iet +2. (в) ) + 


Учитывая, что tuod > lH fgg > At, упрощаем 
уравнение 


{тая 2 


2 2 
22 thog ~ trast = 0. 


Решая данное биквадратное уравнение, полу- 
чаем 


нов = 


где 


2 
D = (=) + trast%  trast2.. 
Ha рис. 7 приведены результаты расчета по 
формуле (7) значения времени хранения (ось ор- 
динат), соответствующие достижению доверитель- 
ным интервалам trasto = 20 км и trasto = 40 км, 
в зависимости от доверительного интервала по 
ускорению (ось абсцисс). Доверительный интервал 
по координатам обобщенной трассы на момент по- 
следнего ее обновления принят равным 250 м. 
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Рис. 5. Результаты табличного расчета значения времени хранения, соответствующего достижению доверительного 
интервала trasto = 20 км 
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Рис. 6. Результаты табличного расчета значения времени хранения, соответствующего достижению доверительного 
интервала trasto = 40 км 
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Рис. 7. Результаты расчета по формуле (7) значения времени хранения по достижению доверительного интервала 
(таѕіс = 20 км и trasto = 40 км в зависимости от доверительного интервала по ускорению 


Сравнение данных, приведенных на рис. 7 Список литературы 


с данными, приведенными на рис. 6 и 5, показыва- 
ет практическое совпадение результатов, получен- 
ных табличным методом и методом приближенных 
расчетов по формуле (7). 


Заключение 


Результаты расчетов параметров и коэффици- 
ентов программ третичной обработки информации, 
приведенные в работе, близки к эмпирическим зна- 
чениям данных параметров и коэффициентов, при- 
нятым в реализованных системах обработки ин- 
формации для частных ограниченных условий при- 
менения. 

Полученные методы позволяют гибко использо- 
вать программы третичной обработки информации, 
расширяя возможность их применения в условиях 
разрывного поля видимости источников информа- 
ции, их существенной разноточности и неопределен- 
ности гипотез движения объектов. Использование 
данных методов в специальном программном обес- 
печении перспективных центров обработки инфор- 
мации повысит достоверность формируемого описа- 
ния текущей объектовой обстановки. 
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Микрополосковый делитель мощности 
с тремя входами и тремя выходами 
для бортовой аппаратуры космических аппаратов 
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Аннотация. Представлены устройство, конструкция и экспериментальные электрические параметры микрополоскового делите- 
ля мощности с тремя входами и тремя выходами (3х 3), содержащее три направленных ответвителя (НО) Ланге, соединенных 
между собой тремя отрезками микрополосковой линии передачи. Приведены требуемые для достижения оптимальных парамет- 
ров значения коэффициентов передачи каждого НО и длин отрезка микрополосковой линии передачи. При подаче СВЧ-сигнала 
на любой из входов делителя мощности 3 x 3 теоретическое значение коэффициента передачи на каждый из трех его выхо- 
дов (без учета потерь в НО и в отрезках микрополосковой линии передачи) — не менее —4,78 дБ; при этом обеспечивается 
развязка между входом и двумя другими входами, а также между выходами не менее 20 дБ в 30-процентной полосе частот. 


Ключевые слова: микрополосковый делитель мощности, направленный ответвитель Ланге, коэффициент передачи, развязка 
между входами и между выходами 


Microstrip Power Divider 
with Three Inputs and Three Outputs 
ior Spacecrait Onboard Equipment 
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S. A. Yakhutin, otdeleniel /(Ospacecorp.ru 
Joint Stock Company "Russian Space Systems", Moscow, Russian Federation 


Abstract. The paper presents a device, design, and experimental electric parameters of a microstip power divider (3 x 3) with three 
inputs and three outputs consisting of three Lange directional couplers connected to each other by three sections of a microstrip trans- 
mission line. The article gives the required values of the transmission coefficients of each Lange directional coupler and lengths of a sec- 
tion of a microstrip transmission line to obtain optimum parameters. When sending a microwave signal to any of power divider (3 x 3) 
inputs, a theoretical value of the transmission coefficient to each of its three inputs (without including losses in the Lange directional 
coupler in the sections of a microstrip transmission line) is not less than —4.78 dB. At the same time decoupling between the input 
and two other inputs is provided, as well as between the outputs, which is not less than 20 dB in a 30% frequency band. 


Keywords: microstrip power divider, Lange directional coupler, transmission coefficient, inputs and outputs decoupling 
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Введение 


В бортовой аппаратуре космических аппа- 
ратов требуется соединять один из нескольких 
СВЧ-узлов или приборов с одним или нескольки- 
ми другими СВЧ-узлами или приборами. Таким 
образом, обеспечивается, например, резервирова- 
ние СВЧ-устройств. В качестве устройств пере- 
крестного резервирования применяются делители 
мощности, чаще всего с двумя или тремя входами 
и выходами (2 x 2, 3 x 2 или 3 x 3). Для борто- 
вой аппаратуры основными требованиями к дели- 
телям мощности являются: максимально возмож- 
ный коэффицент передачи, равномерность деления 
мощности между выходами, относительно большая 
полоса рабочих частот, развязка между входами 
и между выходами, компактность и удобство ин- 
теграции с другими СВЧ-узлами, для чего входы 
и выходы делителей мощности должны быть об- 
ращены к противоположным сторонам диэлектри- 
ческой подложки. Для ряда применений делителей 
мощности существенно их малое групповое время 
запаздывания в рабочем диапазоне частот. Извест- 
ные делители мощности размерностью 3 x 2 (2 x 3) 
и 3x3 [1,2] не удовлетворяют по крайней мере 
одному из перечисленных требований. 

Известны делители мощности с различным ко- 
личеством входов и выходов. Делители мощности 
2х2 реализуются на 3 дБ мостах различного типа 
и имеют теоретический (без учета потерь) коэффи- 
циент передачи —3 дБ [1,2]. 

Делители мощности 2 x З (3 x 2) в боль- 
шинстве строятся на основе каскадного включения 
3-децибельного моста и делителя мощности l x 2, 
подключаемого своим входом к одному из выхо- 
дов моста [3,4]. Недостаток такого делителя мощно- 
сти — неравмерность деления: на одном из выходов 
коэффициента передачи составляет порядка —3 ДБ, 
а на двух других выходах — не менее —6 дБ. 

Если необходимо иметь у делителя мощности 
три входа и три выхода, то рациональным реше- 
нием считается делитель мощности, изображенный 
на рис. 1 [5]. 

Делитель мощности содержит четыре моста 
Ланге, расположенных таким образом, что они 
в проекции образуют квадрат. Внутренние концы 
мостов Ланге соединены между собой, а внешние 


Вход 1 (Выход 1) Вход 2 (Выход 2) 


Выход 3 = Выход 1 
(Вход 3) X (Вход 1) 
X 
Выход 4 X R2 
(Вход 3) >= 
Вход 3 (Выход 3) RI 


Рис. 1. Схема делителя мощности 4 x 4 (3 x 3) [4] 


являются соответственно четырьмя входами и че- 
тырьмя выходами. Коэффициент передачи с любого 
входа на любой выход (без учета потерь в мо- 
стах Ланге) составляет —6 дБ. Если один из вхо- 
дов и один из выходов нагрузить на согласованную 
нагрузку в виде резистора (50 Ом), получаем де- 
литель мощности З х 3; при том же коэфициенте 
передачи, нагружая один из других входов или вы- 
ходов делителя, можно получить делители мощно- 
сти 2 x 3 (3 x 2). 

Недостаток такого делителя мощности — ма- 
лый коэффициент передачи и то, что входы и выхо- 
ды делителя мощности не могут быть расположе- 
ны на противоположных сторонах диэлектрической 
подложки. Для того чтобы обеспечить требуемое 
расположение, необходимо ввести на диэлектриче- 
ской подложке перекрестие входных и выходных 
отрезков микрополосковой линии, что, как пра- 
вило, приводит к уменьшению коэффициента пе- 
редачи. 

Предлагается новый микрополосковый дели- 
тель мощности 3 x 3 (2 x 3) [5], удовлетворяющий 
перечисленным выше требованиям. 


Устройство делителя мощности 3 x 3 


Схема предложенного делителя мощности 3x3 
представлена на рис. 2 [6]. 

Делитель мощности 3 x 3 выполнен на диэлек- 
трической подложке 1, на которой расположены 
три микрополосковых НО Ланге (НО1, HO2, НОЗ) 
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(Выход 8) | Вход 3 (Выход 2) li Вход 2 (Выход 1) 7 Вход 1 
1 (1) 4 (1) 
2 3 (1) 2 (1) 
HOI 
1 (3) 4 (3) 
3 (3) 2 (3) : 
HO3 
3 
1 (2) 4 (2) 
3 (2) -— 2 (2) 
(Вход 3) Выход 3 (Вход 2) { Выход 2 (Вход 1) | Выход 1 


Рис. 2. Схема делителя мощности 3 x 3 


Hà связанных линиях, однократно перекрещиваю- 
щихся, и при этом НОІ и НО2 имеют переход- 
ное ослабление З дБ, а НОЗ — 4,78 дБ. Каждый из 
двух концов связанных линий НОЗ, 4 (3) и 2 (3) 
соединен одним из двух одинаковых отрезков мик- 
рополосковой линии (МПЛ) 2 и З с одним из KOH- 
HOB отрезков связанных линий НОІ и НО? — 3 (1) 
и 1 (2); концы других отрезков связанных линий, 
расположенных с другой стороны НОТ и НО? — 
2 (1) ид (2), соединены между собой отрезком 
МПЛ 4, причем отрезки связанных линий HOI 
и НО2 соединены с отрезками МПЛ 2, 3, 4 так, что 
их свободные концы 4 (1), 1 (1) u 2 (2), 3 (2) o6pa- 
щены к противоположным краям диэлектрической 
подложки. Таким образом, концы 4 (1), 1 (1) и 1 (3) 
являются входами (выходами), а концы 2 (2), 3 (2), 
3 (3) — выходами (входами) делителя мощности 
3 x 3. Входной сигнал подается на одни из входов: 
Вход 1, или Вход 2, или Вход 3. 

С учетом величин модуля и фазы при распро- 
странении входного сигнала с одного из входов на 
выходы обеспечиваются одинаковые коэффициен- 
ты передачи. Для этого длина отрезка МПЛ 4 


должна быть больше суммарной длины отрезков 
МПЛ 2 и 3 на величину от 0,19 до 0,21 длины 
волны на средней частоте рабочего диапазона. Рас- 
смотрены пути прохождения входного сигнала на 
выходы делителя мощности 3 x 3. Коэффициент пе- 
редачи со Входа 2 или со Входа 3 на Выходы 1, 2, 
3 имеет то же самое значение. 

Коэффициент передачи на средней частоте 
диапазона с любого из входов на любой из выхо- 
дов равен —4,78 дБ без учета потерь в делителе 
мощности З X 3. 

Отрезки МПЛ 2 и З вводятся для того, чтобы 
оптимизировать развязки между входами делителя 
мощности. Их длина может быть выбрана в пре- 
делах от 0,125 до 0,25 длины волны на средней 
частоте рабочего дипазона. Приведенные парамет- 
ры делителя мощности 3 х 3 получены в результате 
моделирования в программе AWR Design Environ- 
ment. 

Подключение согласованной нагрузки (50 Ом) 
к одному из входов (выходов) делителя мощности 
3 x 3 превращает ero в делитель мощности 2 x 3 
(3 x 2). 
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Конструкция и результаты 
экспериментальных исследований 
делителя мощности З х З 


Макет делителя мощности 3X3, предназначен- 
ный для работы в С-диапазоне частот, выполнен на 
подложке из поликора толщиной 1 мм. Топология 
макета делителя мощности 3 х 3 представлена на 
рис. 3. 


Габариты подложки делителя мощности 3 x 3 
не превышают 24 x 30 мм. Делитель мощности З х 
х 3 может быть выполнен как в виде отдельного 
прибора, так и в виде функционального узла, инте- 
грированного в один из приборов бортовой аппара- 
туры, например в резервированный троированный 
усилитель мощности в качестве узла перекрестного 
резервирования. 

На рис. 4-8 приведены основные эксперимен- 
тальные характеристики делителя мощности З х 3 


Рис. 3. Топология макета делителя мощности 3 x 3 


Eo oc oU 


Коэффициент передачи, дБ 
MET Вход 1- us 
Выход 1 
=4 
—5 
xd Вход 1- 
Выход 2 
—6 
=r 
... Вход l- 
Выход 3 —8 
2,8 


3,3 
Частота, ГГц 


3,8 


Рис. 4. Зависимость коэффициента передачи делителя мощности 3 x 3 со Входа 1 на Выходы 1, 2, 3 от частоты 
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С-диапазона частот, полученные в нормальных yC- пового времени запаздывания от частоты. Как вид- 
ловиях: зависимости коэффициента передачи с од- но из рис. 4, вносимые потери делителя мощности 
ного из входов на каждый из трех выходов, развяз- 3 x 3 находятся в диапазоне от 4,9 до 5,5 дБ. 

ки между входами, коэффициента стоячей волны Развязка между входами (рис. 5) находится 
по напряжению (КСВН) входов и выходов и груп- в диапазоне от 19 до 32 дБ. 


Коэффициент передачи, дБ 
__ Вход 1- 
Вход 2 
—10 
—20 
ni Вход 1- 
Вход 3 
—30 
40 
__ Вход 2- 
Вход 3 —50 
2,8 3,3 3,8 
Частота, ГГц 


Рис. 5. Зависимость коэффициента передачи делителя мощности З x З со Входа 1 на Входы 2 и З со Входа 2 
на Вход 3 от частоты 


КСВН 
— Выход 1 
— Выход 2 1,5 
— Выход З 1 
28 3,3 3,8 
Частота, ГГц 


Рис. 6. Зависимость КСВН Входов 1, 2 и З делителя мощности З x З от частоты 
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КСВН входов и выходов (рис. би 7) не превы- Заключение 


шает 1,5 дБ. 
Групповое время запаздывания находится Представленный микрополосковый делитель 
в диапазоне от 0,5 до 0,58 нс в полосе частот от 2,8 мощности 3 x 3 имеет коэффициент передачи 
до 3,8 Гц (рис. 8). на 1,2 дБ выше, чем у ближайшего аналога [4]. 
КСВН 
2 
— Выход 1 
— Выход 2 1,5 
— Выход 3 1 
2,8 3,3 3,8 
Частота, ГГц 


Рис. 7. Зависимость KCBH Выходов 1, 2 и З делителя мощности З x З от частоты 


— Вход 1-Выход 1 ГВЗ, нс 
0,6 
— Вход 1-Выход 2 
0,58 
— Вход 1-Выход 3 
0,56 
— Вход 2-Выход 1 
4 
— Вход 2-Выход 2 o 
0,52 


— Вход 2-Выход 3 


— Вход 3-Выход 1 0,5 


Вход 3-Выход 2 | 0,48 
2,8 3,3 3,8 
—— Вход 3-Выход 3 Частота, ГГц 


Рис. 8. Зависимость группового времени запаздывания делителя мощности 3 X 3 с каждого из входов на каждый 
из выходов от частоты 
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Делитель мощности является компактным устрой- 
ством перекрестного резервирования с удобным 
расположением входов и выходов на противополож- 
ных сторонах делителя мощности 3 x 3. Он най- 
дет применение в бортовой аппаратуре командно- 
измерительной системы космического аппарата, бу- 
дучи установленным между выходами приемо-пе- 
редающего устройств и входами усилителей мощ- 
ности. 
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Построение многолучевых антенн 
абонентской радиолинии в $- и Ка-диапазонах частот 


А.А. Аджибеков, contact(Qspacecorp.ru 
АО «Российские космические системы», Москва, Российская Федерация 


Аннотация. Работа посвящена разработке концепции построения многолучевых антенно-фидерных устройств абонентской 
радиолинии в $ и Ка-диапазонах частот. Путем построения многолучевой антенны в виде цифровой антенной решетки обес- 
печивается формирование множества лучей и достигается плотное покрытие подспутниковой зоны в секторе углов +50°. 
Диаграммообразование осуществляется цифровыми методами, что обеспечивает возможность гибкого оперативного изменения 
формы и направления главного максимума ДН каждого луча без изменения архитектуры антенной решетки. 

В статье показано, что для удовлетворения заданным исходным данным, антенны целесообразно построить в виде циф- 
ровых антенных решеток: приемопередающей co 108 элементами в $-диапазоне частот, приемной и передающей co 108 sje- 
ментами в Ка-диапазоне частот. При этом каждая антенная решетка формирует 61 луч в секторе +50° с коэффициентом 
усиления 25,5 дБ и уровнем боковых лепестков менее —20 дБ. 


Ключевые слова: антенно-фидерное устройство, многолучевая антенна, диаграмообразующая схема, антенная решетка, диа- 
грамма направленности, коэффициент усиления 


Building ої Multibeam Antennas 
oi a Subscriber Radio Link т S- and Ka-Frequency Bands 


А.А. Adzhibekov, contact(Ospacecorp.ru 
Joint Stock Company "Russian Space Systems", Moscow, Russian Federation 


Abstract. The paper is devoted to the development of the concept for building multibeam antenna-feeder devices of a subscriber 
radio link in S- and Ka-Írequency bands. By building а multibeam antenna in the form of a digital antenna array, multiple beams 
are formed and tight coverage of a subsatellite zone is achieved in the angle sector +50°. Beam formation is performed by digital 
methods, which provides possibility of flexible operational change of shape and direction of the principal maximum of radiation 
pattern of each beam without change in antenna array architecture. 

The article shows that to meet the specified initial data, antennas are advisable to build in the form of digital antenna arrays: 
transceiver with 108 elements in the S-band, receiving and transmitting with 108 elements in the Ka-band. Each antenna array 
forms 61 beams in the sector 50° with the 25.5 dB gain and the side lobe level not less than —20 dB. 


Keywords: antenna-feeder device, multibeam antenna, beam-forming network, antenna array, radiation pattern, antenna gain 
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В настоящее время в наиболее развитых стра- 
нах мира активизировались работы по созданию 
глобальных спутниковых систем связи на осно- 
ве низкоорбитальных космических аппаратов (КА). 
Использование спутниковых группировок в низко- 
орбитальном сегменте для построения систем пе- 
редачи данных имеет ряд преимуществ по срав- 
нению с аналогичными системами на основе КА, 
размещенных на геостационарных орбитах, в том 
числе низкие задержки, возможность обеспечить 
устойчивую работу в области полюсов, менее жест- 
кие требования к энергетическому потенциалу ап- 
паратуры потребителя [1-3]. Тем не менее постро- 
ение таких систем требует создание антенно-фи- 
дерных устройств (АФУ) абонентской радиолинии, 
обеспечивающих формирование нескольких десят- 
ков (сотен) лучей, обеспечивающих плотное по- 
крытие в секторе углов до +50-60°. Многолучевые 
зеркальные антенны, широко применяемые при по- 
строении спутниковых систем связи в высокоорби- 
тальном и среднеорбитальном сегментах, затрудни- 
тельно применять на низких орбитах ввиду зависи- 
мости сектора сканирования от ширины диаграм- 
мы направленности (ДН) [4-6]. Согласно [6] при 
применении таких антенн может быть достигнут 
сектор сканирования до 10-20 ширин ДН, однако 
для обеспечения широкого сектора обзора вспомо- 
гательный рефлектор двухзеркальной антенны дол- 
жен быть увеличен практически до размера основ- 
ного зеркала, что приведет к увеличению массога- 
баритных параметров. Кроме того, в этом случае 
ширина луча по уровню —3 дБ на краю сектора 
обзора увеличивается практически в два раза. 

Целью работы являлась разработка концепции 
построения бортовых многолучевых антенн або- 
нентской радиолинии, работающих как на прием, 
так и на передачу, формирующих 61 луч в секто- 
ре +50° в $ и Ка-диапазонах частот. 


Синтез пассивной многолучевой 
антенной решетки 


Рассмотрим возможности построения борто- 
вой многолучевой антенны, удовлетворяющей при- 
веденным выше исходным данным. Для этого ис- 
следуем возможности ее построения в виде много- 


лучевой антенной решетки (АР) с пассивной диа- 
граммообразующей схемой (ДОС) без потерь [7,8]. 
Известно, что АР с пассивной ДОС без потерь мо- 
жет сформировать несколько лучей одновременно, 
если диаграммы направленности (ДН) этих лучей 
удовлетворяют условию ортогональности [7,8]: 


IBIUMZO 49 = à, (1) 


где F;(0), Р.(0) — нормированные ДН лучей; д2; — 
символ Кронекера. 

Таким образом, в многолучевой АР с пассив- 
ной ДОС без потерь ДН лучей не могут иметь про- 
извольную форму. Каждому ортогональному лучу 
при этом соответствует свое ортогональное ам- 
плитудно-фазовое распределение полей и, соот- 
ветственно, направление главного максимума ДН, 
а также уровень боковых лепестков (УБЛ). Среди 
ортогональных амплитудных распределений есть 
как равномерное распределение, так и распределе- 
ния, спадающие к краям. Для линейных и квад- 
ратных АР, к примеру, семейство ортогональных 
распределений имеет вид косинуса целой неот- 
рицательной степени. При этом наибольший уро- 
вень пересечения соседних лучей составляет около 
—4 дБ при УБЛ около —13 дБ, что не всегда яв- 
ляется приемлемым с точки зрения электромагнит- 
ной совместимости [9]. Использование косинусного 
распределения позволяет снизить УБЛ до величи- 
ны —23 дБ, однако при этом уровень пересечения 
соседних лучей снижается до величины —9 дБ [8]. 

Снижение УБЛ с сохранением (или минималь- 
ным снижением) уровня пересечения соседних лу- 
чей возможно при построении АР с контуром апер- 
туры, отличным от квадратного [10-13]. Для опре- 
деления вида контура апертуры и шага излучате- 
лей решетки, формирующих ортогональные лучи 
с заданным набором направлений, воспользуемся 
методикой, изложенной в работах [10-13]. Суть 
этой методики в том, что в пространстве направ- 
ляющих косинусов и, v (и = cos(0)sin(p), v = 
= sin(0) sin(q), где диф — полярный и азимуталь- 
ный углы сферической системы координат) опреде- 
ляется область, которую необходимо накрыть луча- 
ми. В этом пространстве выбирается сетка, в узлах 
которой будут расположены направления лучей, 
а также, исходя из их необходимого количества, 
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Рис. 1. Расположение лучей в пространстве направляющих косинусов (слева) и расположение излучателей (справа) 


Таблица. Характеристики бортовых антенн S- и Ка-диапазонов частот 


S-nuanason. Ка-диапазон 


Прием, передача 


передача прием 


61 61 61 
108 108 108 


Треугольная | Треугольная 
1116 x 966 111,2 х 96,6 | 74,4 x 64,4 


Количество лучей 


Количество элементов АР 


Расстояние между элементами, мм 


Тип сетки 


Поперечные габаритные размеры АР, мм 
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Рис. 2. Рассчитанная диаграмма направленности центрального луча антенных решеток S- и Ка-диапазонов частот 
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определяются угловые расстояния между соседни- 
ми лучами. С использованием следующих ниже со- 
отношений определяются: линейные размеры апер- 
туры в плоскости rj; параметры сетки, в узлах KO- 
торых располагаются излучатели антенной решет- 
ки; расстояния между ними: 


А 
А 
А 
d, = AS (4) 
À 


где D, — линейный размер апертуры по оси T; 


D, — линейный размер апертуры по оси у; 
d, — шаг решетки вдоль оси =; 
d, — шаг решетки вдоль оси у; 


Ôu д, — расстояния между соседними лучами 
по осям и и V в пространстве направляющих KOCH- 
нусов; 

A д, — параметры определяющие угловые 
направления дифракционных максимумов антен- 
ной решетки. 

Геометрия контура апертуры выбирается из сле- 
дующего правила: при многократной трансляции 
в плоскости ту этот контур может покрыть всю эту 
плоскость без пустот и пересечений. На рис. 1 про- 
иллюстрированы положения максимумов ДН АР 
в пространстве направляющих косинусов (круг на 
рисунке слева соответствует области углов 50°), 
формирующей 61 ортогональный луч; а также рас- 
положение излучателей в плоскости ху и контур 
апертуры АР. 

В таблице приведены рассчитанные по методи- 
ке [10-13] характеристики AP в Su Ка-диапазонах 
частот. 

Рассчитанная диаграмма направленности цен- 
трального луча антенной решетки представлена на 
рис. 2. 

Коэффициент усиления (КУ) антенной решет- 
ки составил 26,2 дБ, ширина ДН по уровню — 
—3 дБ — 8,9°. 

Как видно из рис. 2, УБЛ составляет —16- 
—17 дБ при уровне пересечения соседних лучей — 
—4,5 дБ [10-13]. При этом, согласно известной 


теореме из теории многолучевых антенн [7,8], 
снижение УБЛ в АР с пассивной ДОС без по- 
терь (при неизменной конфигурации АР) возможно 
только при соответствующем снижении уровня пе- 
ресечения соседних лучей. Таким образом, возни- 
кает противоречие: невозможно одновременно по- 
лучить плотное покрытие заданной области лучами 
при низком УБЛ. 


Построение облика многолучевой 
антенной решетки с цифровым 
диаграмообразованием 


Разрешить возникшее противоречие возможно 
путем использования цифровых методов диаграм- 
мообразования [14, 15]. Сформулированная для 
многолучевых антенн с пассивными ДОС без по- 
терь теорема к АР с цифровым диаграммообразова- 
нием неприменима [8, 14, 15] (ввиду того, что циф- 
ровые ДОС пассивными не являются), так что при 
использовании цифровых методов диаграммообра- 
зования возможно формирование лучей практиче- 
ски с любым сочетанием уровней пересечения со- 
седних лучей и УБЛ. Функциональные схемы при- 
емного и передающего модулей АР с цифровым 
диаграммообразованием представлена на рис. 3. 
Приемный модуль состоит из малошумящего уси- 
лителя МШУ, преобразователя частоты ПЧ, квад- 
ратурного демодулятора КДМ и аналого-цифровых 
преобразователей АЦП. Передающий модуль со- 
стоит из усилителя мощности УМ, квадратурного 
модулятора КМ и цифро-аналогового преобразова- 
теля ЦАП. 

Рассмотрим возможности формирования ДН 
с пониженным УБЛ в АР с параметрами, приве- 
денными в таблице 1. Для этого аналитически по- 
строим ДН АР с помощью следующего соотноше- 
ния [7]: 


FO, p) = (0, p) У Aix 


x ei (аз) собе) у, sin(8)sin(g))-i, (6) 


где 9 — полярный угол сферической системы Koop- 
динат; 

ф — азимутальный угол сферической системы 
координат; 
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Рис. 3. Функциональные схемы приемного (слева) и передающего (справа) модулей антенной решетки с цифровым 
диаграммообразованием [14] 
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Рис. 4. Рассчитанная диаграммы направленности центрального луча антенной решетки с распределением XaHceHa 


700,0) — ДН слабонаправленного излуча- 


теля; 


А; — амплитуда возбуждения 1-го излучателя; 


0; — фаза возбуждения 1-го излучателя; 


т, у — координаты 1-го излучателя; 


k — волновое число; 


1 — мнимая единица. 


В качестве амплитудного распределения вы- 


берем однопараметрическое распределение Хансе- 
на [8]. На рис. 4 представлено сечение рассчитан- 
ной с помощью (6) ДН АР. Как видно из рисунка, 
КУ составляет 25,6 дБ, при УБЛ —28 дБ и ширине 
ДН по уровню —3 дБ 9-10°. При этом могут быть 
обеспечены направления лучей, определенные ра- 
нее при синтезе АР, представленые на рис. 1, что 
позволяет сохранить уровень пересечения лучей, 
полученный для АР с пассивной ДОС. 


Выводы 


В статье дано обоснование способа формирова- 
ния облика бортовых антенн абонентской радиоли- 
нии в S и Ка-диапазонах частот. Показано, что для 
удовлетворения заданным исходным данным эти 
антенны целесообразно построить в виде цифровых 
антенных решеток: приемопередающей со 108 эле- 
ментами в 5-диапазоне частот; приемной и nepe- 
дающей со 108 элементами в Ка-диапазоне частот. 
При этом каждая АР формирует 61 луч в секторе 
X50? с коэффициентом усиления 25,5 дБ и УБЛ 
менее —20 дБ. 
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Аннотация. В рамках разработки опытного образца приемо-передающего модуля (ППМ) активной фазированной антенной 
решетки (АФАР) на основе комбинирования технологии многослойных печатных СВЧ-плат и управляющей системы на кристал- 
ле был разработан монолитный малошумящий усилитель (МШУ) для диапазона частот 30-36 ГГц, созданный по технологии 
РНЕМТ GaAs с топологической нормой 0,15 мкм. Разработка выполнена с помощью современных систем автоматического про- 
ектирования и моделирования. Проведены исследования по оптимизации топологии усилителя для ППМ АФАР. Результаты рас- 
четов были использованы для проектирования схемы и топологии усилителя в средах Microwave Office и Advanced Design System. 

В статье показано, что на частоте 33 ГГц коэффициент усиления МШУ достигает 15 дБ при коэффициенте шума не бо- 
лее 3 дБ. 


Ключевые слова: монолитная интегральная схема, малошумящий усилитель, рНЕМТ, активная фазированная решетка 


Low-Noise Amplifier ior the Transceiver Module 
oi a Ka-band Active Phased Array 


A. V. Buyankin, buyankin.av(Qspacecorp.ru 
Joint Stock Company "Russian Space Systems", Moscow, Russian Federation 


V. E. Poimalin, poimalin.ve(Qspacecorp.ru 
Joint Stock Company "Russian Space Systems", Moscow, Russian Federation 


А.А. Nelin, Сапа. Sci. (Mil.), innovation(Qspacecorp.ru 
Joint Stock Company "Russian Space Systems", Moscow, Russian Federation 


L. E. Ragulina, Cand. Sci. (Phys.-Math.), ragulina.le(Qmipt.ru 
Moscow Institute of Physics and Technology, Dolgoprudny, Moscow region, Russian Federation 


М. V. Ryzhakov, mryzhakov(Qapplmech.mipt.ru 
Moscow Institute of Physics and Technology, Dolgoprudny, Moscow region, Russian Federation 


Abstract. А monolithic low-noise amplifier (ГМА) for the 30-36 GHz frequency band with the implementation of pHEMT GaAs 
technology with а 0.15 um layout rule was created as part of the process of designing a prototype transceiver module (TM) 
for an active phased array antenna (APAA) based on the combination of multilayer microwave (MW) printed circuit boards and 
control system-on-chip technologies. The development was carried out with the help of present-day systems ої automated design 
and simulation. Research on optimizing the amplifier layout for the APAA TM was carried out. Results of calculations were used 
to design the circuit and layout of the amplifier in the Microwave ОШсе and Advanced Design System environments. 

The article demonstrates that LNA gain reaches 15 dB with a noise ligure of no more than З dB at a frequency of 33 GHz. 
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Введение 


В настоящее время наблюдается стремитель- 
ное развитие космических систем радиолокации 
и связи, что влечет за собой существенное уве- 
личение количества космических аппаратов (КА) 
и увеличение объемов передаваемой информации 
по каналам связи. При этом частотный ресурс 
в разрешенных диапазонах P, L, S, C, X весь- 
ма ограничен и B большинстве случаев занят суще- 
ствующими системами. В связи с этим перед разра- 
ботчиками перспективных космических систем свя- 
зи стоит задача освоения миллиметровых диапазо- 
нов К- и Ка- для связи как между КА, так и меж- 
ду КА и наземными станциями [1]. При этом одной 
из современных тенденций является широкое при- 
менение активных фазированных антенных решеток 
(АФАР) в космических радиолиниях. 

Актуальность работы подтверждается прово- 
димыми в настоящее время комплексными про- 
ектами по созданию высокотехнологичных произ- 
водств, в том числе и приемо-передающих модулей 
с активной фазированной антенной решеткой на ос- 
нове современных технологий. Проект от 31.05.2018 
№074-11-2018-014, выполняемый в рамках согла- 
шения с Минобрнауки России и кооперации с вуза- 
ми (Постановление Правительства Российской Фе- 
дерации от 9 апреля 2010 г. №218), — тому под- 
тверждение. Решаемая комплексная задача при этом 
является сложной как с научной, так и с производ- 
ственной и технологической точек зрения. 

Одна из главных задач при создании приемо- 
передающих модулей Ка-диапазона частот — уве- 
личение чувствительности приемника [2]. Для ре- 
шения данной задачи требуется создание монолит- 
ной интегральной схемы (МИС) МШУ с достаточ- 
но малым коэффициентом шума. Наиболее подхо- 
дящей для создания МШУ Ка-диапазона частот 
является технология pHEMT GaAs с топологиче- 
ской нормой 0,15 мкм [3]. В связи с необходи- 
мостью интеграции МШУ непосредственно в ан- 
тенну к блоку МШУ могут быть предъявлены бо- 
лее жесткие требования по внешним воздействую- 
щим факторам (ВВФ), в особенности по спецфак- 
торам. Поэтому технология создания МИС должна 
быть в наименьшей степени подвержена воздей- 
ствию космических излучений и др. [4]. 


Таблица. Требуемые параметры МШУ 


Норма 
Наименование Буквенное параметра 
параметра, обозначение 
единица измерения | параметра He 
более 
36 
Коэффициент 
усиления, дБ 
Коэффициент 31 
шума, дБ у 
Максимальная 
выходная мощность, 
дБмВт 
КСВН входа/выхода 2,0 
6 | Ток потребления, мА " 80 


Требуемые параметры 
в таблице. 

На основе процесса pHEMT 0,15 мкм разра- 
ботан 2-каскадный МШУ с напряжением питания 
5 В и напряжением смещения 0,7 В (рис. 1). 

Принципиальная схема МШУ при реализации 
на микрополосковых линиях представлена на рис. 2. 

На основе принципиальной схемы сгенериро- 
вана топология усилителя (рис. 3). Размер МИС 
МШУ 2500 x 1000 мкм. Расчетный ток потребле- 
ния — менее 80 мА. 

В процессе разработки МШУ проведено моде- 
лирование выходной мощности. В результате до- 
стигнуто значение выходной мощности МШУ на 
уровне 15 дБм. На рис. 4 представлен гармониче- 
ский сигнал усилителя при максимальной мощно- 
сти на выходе. 

Моделирование топологии было произведено 
методом косимуляции (рис. 5), смысл которого со- 
стоит в том, что электромагнитное моделирование 
производится без транзисторов. Модели транзисто- 
ров при получениии характеристик МШУ подгру- 
жаются из схемотехнического редактора. 

Результаты моделирования максимального ко- 
эффициента усиления МИС МШУ в схемотехни- 
ческом и топологическом редакторах представлены 
на рис. 6. 


МШУ представлены 
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Рис. 2. Принципиальная схема МШУ при реализации Hà микрополосковых линиях 


Рис. 3. Топология МИС МШУ 


РАКЕТНО-КОСМИЧЕСКОЕ ПРИБОРОСТРОЕНИЕ И ИНФОРМАЦИОННЫЕ СИСТЕМЫ т. 7 вып. 1 2020 


МАЛОШУМЯЩИЙ УСИЛИТЕЛЬ ПРИЕМО-ПЕРЕДАЮЩЕГО МОДУЛЯ 39 


Спектр на выходе 
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Рис. 4. Значение максимальной выходной мощности МИС МШУ и ближайших гармоник 


Рис. 5. Режим косимуляции МИС МШУ 


Результаты топологического моделирования 
соответствуют предварительным оценкам. Усиле- 
ние, полученное в топологии, в среднем на 2 дБ 
меньше, чем результаты схемотехнического проек- 
тирования. На рис. 7 даны результаты моделирова- 


ния частотной зависимости коэффициентов отра- 
жения по входу и выходу. 


Результаты топологического моделирования 
в части коэффициента отражения по входу ока- 
зались лучше, чем результаты схемотехнического 
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4 Результат топологического 


а моделирования 
p Результат схемотехнического 
r^ 10 | моделирования 
Ul 4 |21 m2 
- | Частота = 32,00 TTu| |Hacrora = 33,00 ГГц 
5 1/5(2,1) = 15,969 дБ 5 (2,1) = 15,394 дБ 


35 
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Рис. 6. Сравнение результатов моделирования коэффи- 
циента усиления (521) МИС МШУ 
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Рис. 7. дБ МШУ 


моделирования. Результаты топологического моде- 
лирования в части коэффициента отражения по вы- 
ходу частично не соответствуют результатам схе- 
мотехнического моделирования, однако не выхо- 
дят за границы требований в частотном диапазоне 
30-34 ГГц. 

На рис. 8 приведены результаты моделирова- 
ния коэффициентов шума МИС МШУ. 


T T T T 
Ки (топологическое моделирование) 


K „ (схемотехническое моделирование) 


Пе 


30 31 32 33 34 35 36 37 
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Рис. 8. Сравнение результатов моделирования коэффи- 
циентов шума МИС МШУ 


По результатам топологического моделировния 
значение коэффициента шума не превышает 3 дБ, 
а в диапазоне частот 30-33 ГГц коэффициент шума 
составляет 2,8 дБ. 


Выводы 


Исходя из результатов проведенного модели- 
рования можно прийти к выводу, что прогноз 
по наиболее оптимальному использованию тех- 
нологии РНЕМТ 0,15 мкм для создания МШУ 
Ка-диапазона частот был подтвержден. Разрабо- 
танная топология МИС МШУ может быть исполь- 
зована с необходимыми конструктивными измене- 
ниями для создания приемного тракта кор-чипа 
ППМ АФАР Ка-диапазона частот. 

Статья подготовлена в рамках соглашения 
АО «Российские космические системы» с Минобр- 
науки России от 31.05.2018 М 074-11-2018-014. 
Соглашение заключено на основании Постанов- 
ления Правительства Российской Федерации от 
9 апреля 2010 e. №218 «О мерах государствен- 
ной поддержки развития кооперации российских 
образовательных организаций высшего образо- 
вания, государственных учреждений и органи- 
заций» для реализации комплексных проектов 
по созданию высокотехнологичного производ- 
ства в рамках подпрограммы «Институциональ- 
ное развитие научно-исследовательского сек- 
тора» государственной программы Российской 
Федерации «Развитие науки и технологий на 
2013-2020 годы». 


РАКЕТНО-КОСМИЧЕСКОЕ ПРИБОРОСТРОЕНИЕ И ИНФОРМАЦИОННЫЕ СИСТЕМЫ т. 7 вып. 1 2020 


МАЛОШУМЯЩИЙ УСИЛИТЕЛЬ ПРИЕМО-ПЕРЕДАЮЩЕГО МОДУЛЯ 41 


Список литературы 3. Wang Gang, Chen Wei, Liu Jiarui, Mo Jiongjiong, 


Chen Hua, Wang Zhiyu, Yu Faxin et al. Design 
1. Гнатюк Д.Л., Федоров IO. B., Галиев T.B., Гали- 


ев P.P., Щербакова М.Ю. Монолитные интеграль- 
ные схемы малошумящих усилителей КВЧ-диапа- 
зона на GaAs рНЕМТ-гетероструктурах // Доклады 
ТУСУРа, 2010, №2 (22), ч. 1. С. 49-55. 


. Yang Hou, Ruming Wen, lingyun Li et al. Analysis and 
optimum design ої impedance matching for Ka-band 
cryogenic low noise amplifiers [C] // Asia Pacific 
Microwave Conference. Singapore, 2009. P. 1593-1596. 


of a broadband Ka-band MMIC LNA using deep nega- 
tive feedback loop // IEICE Electronics Express, 2018, 
vol. 15, iss. 10. 


. Wei-Tsung Li, Jeng-Han Tsai, Hong-Yuan Yang, Wei- 


Hung Chou, Shyh-Buu Gea, Hsin-Chia Lu, Tian- 
Wei Huang et al. Parasitic-Insensitive Linearization 
Methods for 60-GHz 90-nm CMOS LNAs // IEEE 
Transactions on Microwave Theory and Techniques, 
Aug. 2012, vol. 60, no. 8. P. 2512-2523. 


РАКЕТНО-КОСМИЧЕСКОЕ ПРИБОРОСТРОЕНИЕ И ИНФОРМАЦИОННЫЕ СИСТЕМЫ т. 7 вып. 1 2020 


РАКЕТНО-КОСМИЧЕСКОЕ ПРИБОРОСТРОЕНИЕ И ИНФОРМАЦИОННЫЕ СИСТЕМЫ 
2020, том 7, выпуск 1, с. 42-50 


РАДИОТЕХНИКА И КОСМИЧЕСКАЯ СВЯЗЬ 


УДК 621.317.08 РОГ 10.30894/issn2409-0239.2020.7.1.42.50 


Приемный канал СВЧ-гиперспектрометра с двумя 
приемниками для аэрокосмического зондирования Земли 


А. В. Убайчин, к. т. n., с.н. с., доцент, ubaychin(Qtpu.ru 
Национальный исследовательский Томский политехнический университет, г. Томск, Российская Федерация 


Аннотация. Приведена разработанная структурная схема приемного канала СВЧ-гиперспектрометра на основе нового 
СВЧ-радиометра с двумя приемниками. Описан алгоритм функционирования приемного канала СВЧ-гиперспектрометра. По- 
казанный вариант использования двух полупроводниковых генераторов шума для построения передаточной характеристики 
СВЧ-радиометра позволяет уменьшить зависимость результатов измерений от дестабилизирующих факторов радиометрических 
измерений в виде дрейфа собственных шумов и коэффициента передачи радиометрического приемника. 

Предложен способ формирования опорного генератора низкотемпературного шума из собственных шумов малошумящего 
усилителя. Способ основывается на реализации низкотемпературного шума, обусловленного шумовой волной на входе мало- 
шумящего усилителя, работающего в согласованном режиме, при постоянных шумовых параметрах. Данный способ позволяет 
осуществить перенос направленного ответвителя из антенного канала СВЧ-радиометра в канал опорного генератора шума, что 
обусловливает повышение флуктуационной чувствительности за счет снижения потерь на входе системы. 

Представлена математическая модель оценки погрешности измерений как функции точности термостатирования, темпе- 
ратурных коэффициентов и величины потерь в элементах и узлах входной высокочастотной части. Показана возможность 
повышения разрешающей способности по частоте в два раза при заданной флуктуационной чувствительности, длительности 
одного измерения и уровне шумов за счет использования двух работающих параллельно радиометрических приемников. На ос- 
новании разработанной математической модели приведены результаты численного эксперимента по оценке точности измерений 
при различных технических параметрах входной высокочастотной части. 

Использование предложенного СВЧ-радиометра позволяет уменьшить зависимость результатов измерений от точности 
термостатирования по сравнению с существующими СВЧ-радиометрами, использующими пассивные генераторы шума. 


Ключевые слова: дистанционное зондирование, научное приборостроение, температурная стабильность, СВЧ-гиперспектраль- 
ные исследования, погрешность измерений 


Receiving Channel of UHF-Hyperspectrometer 


with Two Receivers ior Earth Remote Sensing 


A. V. Ubaychin, Сапа. Sci. (Engineering), senior researcher, associate professor, ubagchin(Qtpu.ru 
Tomsk Polytechnic University, Tomsk, Russian Federation 


Abstract. The developed block diagram of the receiving channel of a microwave hyperspectrometer based on а new microwave 
radiometer with two receivers is presented. Ап algorithm for the functioning of the receiving channel of the microwave hyperspec- 
trometer is described. The presented method of using two semiconductor noise generators for building the transfer characteristic 
of a microwave radiometer makes it possible to reduce the dependence of the measurement results on the destabilizing factors of 
radiometric measurements in the form of a drift of selí-noise and the transmission coefficient of the radiometric receiver. 

А method ої forming a reference generator ої low-temperature noise from the self-noise of a low-noise amplifier is proposed. 
The method is based on the implementation of low-temperature noise caused by a noise wave at the input of the low-noise amplifier 
operating in a coordinated mode with constant noise parameters. This method allows the transfer of the directional coupler from 
the antenna channel of the microwave radiometer to the channel of the reference noise generator, which leads to an increase 
in fluctuation sensitivity by reducing losses at the input of the system. 

А mathematical model for estimating the measurement error as a function of thermostating accuracy, temperature coefficients, 
and losses in the elements and nodes of the high-frequency input is presented. The possibility of increasing the frequency resolution 
by a factor of two at a given fluctuation sensitivity, duration of one measurement and noise level due to the use of two radiometric 
receivers operating in parallel is presented. Based on the developed mathematical model, the results of a numerical experiment 
to assess the accuracy of measurements at various technical parameters ої the high-frequency input are presented. 

Using the proposed microwave radiometer makes it possible to reduce the dependence of the measurement results on the 
accuracy of temperature control compared to existing microwave radiometers using passive noise generators. 


Keywords: remote sensing, scientific instrumentation, temperature stability, microwave hyperspectral studies, measurement error 
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Введение 


Неуклонный рост потребности к повышению 
качества и количества информации, получаемой при 
дистанционном зондировании, приводит к необходи- 
мости поиска новых методов и подходов к его прове- 
дению. Одним из наиболее перспективных является 
гиперспектральный метод исследований, позволяю- 
щий получить информацию об объекте путем изме- 
рения его излучения в различных участках спек- 
тра. Повышение объема данных о спектральном 
составе излучения приводит к улучшению эффек- 
тивности решения обратных задач характеризации 
исследуемого объекта. 

Наибольшее распространение гиперспектраль- 
ные методы получили в оптическом и ИК-диапа- 
зонах. Это связано с технологической простотой 
реализации большого числа сенсоров, работающих 
в различных спектральных диапазонах [1]. 

Прогресс развития и миниатюризации СВЧ-те- 
хники позволил реализовать гиперспектральные ис- 
следования в СВЧ-диапазоне. СВЧ-гиперспектроме- 
тры состоят из нескольких СВЧ-радиометров. Один 
или несколько СВЧ-радиометров, работающих 
на близких частотах, образуют приемный канал 
СВЧ-гиперспектрометра. Современные СВЧ-гипер- 
спектрометры обладают десятками и сотнями при- 
емных каналов и перекрывают диапазон частот от 
единиц до сотен гигагерц [2]. 

Повышение информативности при СВЧ-гипер- 
спектральных исследованиях обусловливает требо- 
вание по увеличению разрешающей способности по 
частоте. Это требование приводит к необходимости 
уменьшения рабочей полосы частот в применяемых 
СВЧ-радиометрах. Уменьшение рабочей полосы ча- 
стот при заданной шумовой температуре радиомет- 
рического приемника приводит к ухудшению флук- 
туационной чувствительности и увеличению дли- 
тельности измерений. 

Бортовое базирование СВЧ-гиперспектроме- 
тров определяет повышенные требования по устой- 
чивости к изменяющимся факторам внешней сре- 
ды. Современные СВЧ-гиперспектрометры, как 
правило, выполняются на классических СВЧ-ра- 
диометрах. Такие СВЧ-радиометры требуют допол- 
нительной стабилизации радиометрического при- 
емника, что в случаях аэрокосмического базиро- 


вания является фактором, ограничивающим их эф- 
фективное применение. 

В связи с упомянутыми выше особенностя- 
ми реализации СВЧ-гиперспектральных измерений 
актуальной является задача разработки нового под- 
хода к реализации приемных каналов на осно- 
ве СВЧ-радиометров нового типа, позволяющих ре- 
ализовать совокупность повышенного спектрального 
разрешения и устойчивости к изменяющимся факто- 
рам внешней среды при заданном уровне шумов ра- 
диометрического приемника, флуктуационной чув- 
ствительности и длительности одного измерения. 

В статье описан приемный канал СВЧ-ги- 
перспектрометра, выполненный на основе нового 
СВЧ-радиометра, обладающего повышенными раз- 
решающей способностью по частоте и устойчиво- 
стью к изменяющимся факторам внешней среды. 


Теория 


В соответствии с описанными выше требова- 
ниями предложен приемный канал СВЧ-гиперспек- 
трометра, основанный на новом типе СВЧ-радио- 
метров. Разработанный СВЧ-радиометр обладает 
двумя радиометрическим приемниками и использу- 
ет активные генераторы шума в качестве опорных 
точек передаточной характеристики. Структурная 
схема разработанного СВЧ-радиометра представле- 
на на рис. 1. 

Структурная схема разработанного СВЧ-ра- 
диометра состоит из приемной антенны А с шу- 
мовой температурой Т,, термостата ТС с под- 
держиваемой температурой Тст и размещенными 
в нем генератором шума ГШ с шумовой темпера- 
турой Ттит, высокочастотным ключом КЛ, nanpas- 
ленным ответвителем НО с коэффициентом свя- 
зи C, опорным малошумящим усилителем ОМШУ, 
согласованной нагрузки СН, высокочастотным пе- 
реключателем ПК, управляемого гетеродина Гет, 
цифрового блока управления ЦБУ, двухканально- 
го аналого-цифрового преобразователя АЦП, двух 
идентичных радиометрических приемников РПІ 
и PII2 с шумовой температурой Ти, коэффици- 
ентом усиления С, шириной рабочей полосы ча- 
стот Af, состоящих из малошумящего усилите- 
ля МШУ, полосно-пропускающего фильтра ППФ, 
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Рис. 1. Структурная схема разработанного СВЧ-радиометра 


смесителя СМ, усилителя промежуточной частоты 
УПЧ, квадратичного детектора КД, усилителя по- 
стоянного тока УПТ. 

Одной из ключевых особенностей предложен- 
ного СВЧ-радиометра является использование двух 
активных опорных полупроводниковых генераторов 
шума для построения передаточной характеристики. 
Сигнал первого опорного генератора шума форми- 
руется из собственных шумов опорного малошумя- 
щего усилителя. Для реализации этого вход опор- 
ного малошумящего усилителя подключен к входу 
радиометрического приемника. Таким образом, вход 
ОМШУ согласован с входом радиометрического 
приемника. Это обеспечивает отсутствие интерфе- 
ренции и постоянные шумовые параметры ОМШУ. 
При этом на вход радиометрического приемника по- 
ступает низкотемпературный шум, обусловленный 
шумовой волной на входе ОМШУ. Уровень форми- 
руемого на входе ОМШУ шума пропорционален его 
собственным шумам. Этот факт позволяет использо- 
вать ОМШУ в качестве опорного генератора шума 
с шумовой температурой Tog, [3]. 

Второй генератор шума реализован на лавин- 
но-пролетном диоде в режиме обратного пробоя 
и формирует в антенном канале шумовой сиг- 
нал Topo- 

Функционирование предложенного СВЧ-ра- 
диометра осуществляется следующим образом. 
В радиометрических приемниках происходят по- 
следовательные операции усиления, полосовой 
фильтрации, переноса по частоте, усиления по про- 
межуточной частоте, квадратичного детектирова- 
ния и усиления детектированного сигнала. 


ЦБУ формирует сигналы амплитудно-импульс- 
ной модуляции ѓ,имо, lagi И ѓлимә: На рис. 2 
представлены временные диаграммы, поясняющие 
принцип работы предложенного СВЧ-радиометра. 
Ha рис. 2, а изображены сигналы CAMMI им! 
и Флимо Соответственно. На рис. 2,6 и рис. 2, в 


10 


I 
їлимі 


Рис. 2. Временные диаграммы работы СВЧ-радиометра 
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изображены сигналы на выходах РПІ и PII2 co- 
ответственно. Под действием низкого уровня сиг- 
нала [ay Посредством ПК осуществляется nog- 
ключение антенны на вход РП? и входа ОМШУ 
на вход РПІ, под действием высокого уровня сиг- 
нала tapymo антенна подключается на вход РПІ, 
а вход ОМШУ подключается на вход РП2. При вы- 
соком уровне taym] Или ѓ,умә Кл переходит в pe- 
жим низкого ослабления, что приводит к передаче 
сигнала ГШ к дополнительному каналу НО. При 
низком уровне їум Или {имә Кл переходит B pe- 
жим высокого ослабления. Последовательное пере- 
ключение на входах радиометрических приемников 
сигналов антенны, ОМШУ и комбинации ОМШУ 
и ГШ приводит к появлению на их выходе импульс- 
но-модулированной последовательности. 

Величина шумов первого опорного генератора 
шума Тору на входе ПК в соответствии со струк- 
турной схемой (рис. 1) равна 


(1) 


Топ z Tom 


где Tom — величина собственных шумов ОМШУ 
на его входе. Величина шумов второго опорного 
генератора шума Топо на входе ПК равна 


Tono = СТрш, (2) 
где C — коэффициент связи направленного ответ- 
вителя. 

Обозначим уровни сигналов на выходе пер- 
вого радиометрического приемника, соответствую- 
щие подключению на его вход антенны ОМШУ 
и комбинации OMINY и ГШ через Од, Ов и Uc 
соответственно. Уровни Ид, Ив пропорциональны 
шумовым температурам опорных генераторов шума 
и остаются неизменными (теоретически) в процес- 
се эксплуатации. 

Уровень Ис определяется величиной шумо- 
вой температуры антенны. Для его измерения ис- 
пользуются различные методики [4]. В представ- 
ленной реализации используется непосредственное 
измерение уровней U4, Ив и Ос путем оцифровки. 
По сравнению с использованием компаратора и ура- 
внивания вольт-секундных площадей [4] этот спо- 
соб обладает потенциалом по увеличению динами- 
ки измерений при скачкообразных изменениях шу- 
мовой температуры антенны [5]. 


В соответствии с описанным выше алгоритмом 
переключения источников сигнала на входе РПІ 
уровни U4, Up и Ос c учетом (1), (2) равны 

Ол = [Tom + Tu]: G Af- k= 

= [Tom + Tu]: С: Af - К, 
Ug = [Tom + Tm + Trig: C]: G- Af- k — 
= [Tomi + Tono + Tul: G: Af- k, 


(3) 


(4) 
(5) 


где Тош — низкотемпературный шум, пропорцио- 
нальный собственным шумам OMINY, Ти — wy- 
мовая температура РПІ, Т, - шумовая темпера- 
тура антенны, С — полный коэффициент переда- 
чи РПІ, Af — рабочая полоса частот РПІ, k — 
постоянная Больцмана, С — коэффициент связи 
направленного ответвителя. 

Последовательная оцифровка в аналого-циф- 
ровом преобразователе позволяет измерить уровни 
сигналов U4, Св и Ос. Результаты измерений ne- 
редаются в ЦБУ, где осуществляется расчет шу- 
мовой температуры антенны. Расчет шумовой тем- 
пературы антенны происходит исходя из следую- 
mero. Для описанных выше уровней сигналов U4, 
Up и Ос справедливо равенство 


Tono — ТА 
ТА > Tom 


_ Ов - Uc 


0-0. (6) 


Выражая из (6) Th, получим 


pe Tone : [Uc — UA] + Тот : [Up — Uc] (7) 


Выражение (7) связывает измеренные уровни 
сигналов U,, Ив и Uc CO значениями шумовых 
температур опорных источников. Важным след- 
ствием формулы (7) является отсутствие в ней пара- 
метров приемника G и Ту. Этот факт обусловливает 
теоретическое отсутствие влияния на результаты из- 
мерений изменений коэффициента передачи и дрей- 
фа собственных шумов радиометрического приемни- 
ка [4]. Численные значения параметров Топ! и Tono 
определяются при калибровке СВЧ-радиометра по 
двум эталонным источникам шумового излучения 
путем установления соответствия измеренного уров- 
ня сигналов U 4, Upg заданным шумовым температу- 
рам на входе СВЧ-радиометра. Второй радиометри- 
ческий приемник функционирует аналогично РПІ. 
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В аспектах применения предложенного СВЧ-ра- 
диометра в качестве приемного канала СВЧ-ги- 
перспектрометра преимуществом является возмож- 
ность уменьшения рабочей полосы частот в два 
раза по сравнению с классическими СВЧ-радио- 
метрами с одним радиометрическим приемником. 
Этот эффект обусловлен следующим. В работе [6] 
показано, что чувствительность СВЧ-радиометров 
с двумя приемниками в корень квадратный из двух 
выше чувствительности СВЧ-радиометра с одним 
приемником. Аналитически это следует из увели- 
чения числа измерений в два раза на заданном 
временном интервале за счет параллельной рабо- 
ты двух радиометрических приемников. Увеличен- 
ное в два раза число измерений приводит к умень- 
шению дисперсии измеряемой величины в корень 
квадратный из двух. Минимальная чувствитель- 
ность СВЧ-радиометров, использующих два опор- 
ных источника шума и один радиометрический при- 
емник, при верхней границе динамического диапа- 
зона не превышающей Topo, описывается выраже- 
нием [4] 


"ED | 4(Toni + Tin)? + Tón 
"un dfinpT 8 в 


где арр — рабочая полоса частот СВЧ-радиоме- 
тра с одним приемником, т — постоянная времени 
интегратора, характеризующая длительность одно- 
го измерения. 

Учитывая (8), выражение для оценки рабочей 
полосы частот принимает вид 


1 4(Ton + Tin)?  Tóno l 


d = 9 
inp АТ т 8 (9) 


В соответствии с теорией ошибок и исследо- 
ваниями, проведенными в [6], чувствительность 
СВЧ-радиометров с двумя радиометрическими при- 
емниками увеличивается в корень квадратный из 
двух. В этом случае справедливо равенство 


АТ 
АТ, = AL. 10 
T ETC (10) 
где АТ, — флуктуационная чувствительность 


СВЧ-радиометра с двумя приемниками. 
Выражая из (10) AT, и подставляя в (9), за- 
пишем выражение для оценки рабочей полосы ча- 


стот предложенного СВЧ-радиометра с двумя при- 
емниками 


Пра = 59 = дт 16 


Анализ совокупности выражений (8)-(11) поз- 
воляет сделать вывод, что в предложенном СВЧ-ра- 
диометре в случае сохранения длительности изме- 
рений на заданном уровне и флуктуационной чув- 
ствительности на уровне СВЧ-радиометра с одним 
радиометрическим приемником рабочая полоса ча- 
стот сокращается в два раза. 

В процессе эксплуатации СВЧ-гиперспектро- 
метра важным с точки зрения сохранения ста- 
бильности измерений на заданном метрологиче- 
ском уровне является точность поддержания уров- 
ня шумов опорных генераторов шума. В работе [7] 
показано, что использование в качестве опорного 
генератора шума согласованной нагрузки позволя- 
ет достигнуть минимальной точности измерений, 
равной величине ошибки поддержания температу- 
ры термостата. Очевидно, что для улучшения точ- 
ности измерений необходимо применение генерато- 
ров шума, обеспечивающих меньшей зависимостью 
от изменяющейся температуры по сравнению с пас- 
сивными генераторами шума. 

В разработанном СВЧ-радиометре, по срав- 
нению с [8], температурная стабильность выше 
при одинаковой точности поддержания температу- 
ры термостата за счет использования активных ге- 
нераторов шума. Для подтверждения этого прове- 
ден анализ погрешностей радиометрических изме- 
рений, результаты которого приведены ниже. 

Основными факторами, приводящими к по- 
грешности измерений в разработанном СВЧ-радио- 
метре, являются: погрешность поддержания тем- 
пературы термостата дТст, температурная неста- 
бильность шумов ОМШУ дТоц, температурная 
нестабильность шумов генератора шума 9Т-ит, Be- 
личина потерь в линии передачи между ОМШУ 
и HO Іомшу, величина потерь в линии передачи 
между HO и ПК Г. 

Для проведения анализа приняты следующие 
условия. Полные потери в линии передачи меж- 
ду ОМШУ и ПК равны сумме lomy и Ги. Be- 
личиной потерь в линии передачи между антен- 
ной и ПК можно пренебречь, так как вклад этих 
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потерь не поддается контролю и свойствен всем 
СВЧ-радиометрам. Потери в ПК и линии передачи 
между ПК и входом радиометрического приемни- 
ка не влияют на результаты измерений ввиду их 
взаимного вычитания при использовании выраже- 
ния (9). 

Обозначим через 007, и дИв погрешности из- 
мерения уровней сигналов U 4, Ив опорных генера- 
торов шума. Учитывая потери в линиях передачи 
и погрешности поддержания температуры термо- 
стата (формулы (3)-(5), (7)), получим выражения 
для оценки величин ÔU 4 и дОь: 


óU 4 = [6Ттош · |9 Тст| - (Y — [lom + trul) + 
+ |6Тст| - (lom + irm) + Tig]: C- Af - К, (12) 


óUg = [бТтош · |óTer| - (1 — [lom + trul) + 
+ ôTrrm · [6 Тст| : lrum  9Ter - Пош + lru) + Тш] х 
x G-Af-k, (13) 


где ОТтонь» 9Ттги — модули максимального откло- 
нения шумовой температуры ОМШУ и ГШ в пре- 
делах, обусловленных погрешностью поддержания 
температуры термостата. 

Запишем выражение (7) с учетом погрешно- 
стей измерения уровней сигналов U ,, Up: 


ТА — 
_ Tom : [Uc - (UA  9U4)] — Tom [Че - (Usg + ôUp)] 
(Ug --6Ug) — (Ол +804) 


(14) 


Выражение (14) с учетом (12) и (13) позволяет 
аналитически оценить погрешность измерений при 
известных параметрах ЭТст, Тош» 1ш, lomy 


и п: 


Результаты моделирования 


Для оценки абсолютной точности измерений 
разработанного СВЧ-радиометра проведено числен- 
ное моделирование на основе заданных парамет- 
ров входной высокочастотной части. В качестве 
ОМШУ выбран малошумящий усилитель НМС517 
с номинальной шумовой температурой 170 К и тем- 
пературным коэффициентом собственных шумов 
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порядка 0,77 К/°С, для моделирования ГШ исполь- 
зован лавинно-пролетный диод NWA402 c Temne- 
ратурным коэффициентом собственных шумов по- 
рядка 2100 К/°С при шумовой температуре Ти = 
= 9,17. 10? К. 

Анализ температурной нестабильности уров- 
ня шумов ГШ и схемы его подключения поз- 
воляет сделать вывод о следующем. Использова- 
ние направленного ответвителя, обеспечивающего 
направленный перенос в основной канал сигнала 
ГШ, равного верхнему динамическому диапазону 
СВЧ-радиометра, позволяет понизить значение теп- 
ловой нестабильности пропорционально коэффици- 
енту связи С. Для рассмотренного примера тепло- 
вая нестабильность сигнала ГШ при верхней гра- 
нице динамического диапазона, равной Тмдкс = 
= 300 К, в основном канале направленного OTBeT- 
вителя составляет 


T 
ӧТугщ = С. Ттешном =. ЭТтешном = 
ТГШНОМ 
a ор. 103 = 0,69 К/С. 
9,17. 10 


Для проведения численного эксперимента ис- 
пользованы различные точности термостатирова- 
ния входной высокочастотной части и величины 
потерь Іомшу/Ігш: Результаты численного экспе- 
римента на разработанной модели представлены 
в табл. 1. 

Для оценки влияния температурной неста- 
бильности опорных генераторов шума на погреш- 
ность измерений проведено исследование разрабо- 
танной модели при следующих параметрах: точ- 
ность термостатирования дТст = +1 К, величина 
потерь lomy = 0,1, [и = 0,1, шумовая темпера- 
тура антенны 7, = 300 К. Результаты исследова- 
ния приведены в табл. 2. 

Для оценки распределения погрешности из- 
мерений по динамическому диапазону проведено 
моделирование изменения погрешности измерений 
как функции шумовой температуры антенны при 
различных величинах температурной нестабильно- 
сти генераторов шума и заданных параметрах точ- 
ности термостатирования ОТст = +1 К, величине 
потерь юмшу = 0,1, Іп = 0,1. Результаты модели- 
рования представлены в табл. 3. 
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Таблица 1. Модуль максимальной погрешности измерений как функции точности термостатирования и потерь 
во входной высокочастотной части 


Погрешность измерений, К 


Величина потерь, омшу/ lru 


0,5/0,5 | 0,33/0,33 | 0,25/0,25 | 0,2/0,2 | 0,125/0,125 | 0,1/0,1 


0,672 0,576 0,528 0,500 0,457 0,442 


Loan | 


Таблица 2. Модуль максимальной погрешности измерений как функция изменяющейся температурной стабиль- 
ности Топ! и Tono 


Погрешность измерений, К 


Топ 


> 
am 


1,1 


1,4 


1,5 


ое 
o 
ре) 


— 


19 


1,21 £1,3 
0,20 | 0,21 


0,22 


0,23 


0,28 


osz [ossosa [oss [oes osr [oss oes [oro [or [ors 


0,55 


H 
— 
[99] 
© 


| |е [е |е |е |е 
Е е [iens = 


— 


-1,5 


1,61 


0,78 | 0,79 | 0,80 | 0,81 | 0,82 | 0,83 | 06,84 | 0,85 | 0,86 | 0,87 | 0,92 


1,62 | 1,63 | 1,64 | 1,65 | 1,66 | 1,67 | 1,72 


Обсуждение результатов 


Применение в разработанном СВЧ-радиометре 
малошумящего усилителя, подключенного входом 
к входу радиометрического приемника, позволяет 
использовать его в качестве опорного генератора 
низкотемпературного шума. 

Вследствие сравнительно низкого темпера- 
турного коэффициента в примененных опор- 


ных генераторах шума, выполненных на осно- 
ве малошумящего усилителя и лавинно-пролет- 
ного диода, разработанный СВЧ-радиометр обла- 
дает повышенной устойчивостью к изменяющим- 
ся факторам внешней среды при заданной точ- 
ности термостатирования входной высокочастот- 
ной части по сравнению с СВЧ-радиометрами, 
использующими пассивные опорные генераторы 
шума. 
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Таблица 3. Модуль максимальной погрешности измерений как функция изменяющейся шумовой температуры 
антенны и температурной стабильности дТон и Tono 


Погрешность измерений, К 
Температурная стабильность опорных генераторов шума 
0,5 | +0,6 | +0,7 | +0,8 | +0,9 1,11+1,2[+1,3 ; 1,5| +2 
ifo [ous озо [ов [29 [039 [oar [040 [044 [ола [051 [oss олз 
3| 60 | 0,24 | 0,28 | 06,33 | 0,38 | 0,43 | 0,47 | 0,52 | 0,57 | 0,62 | 0,66 | 0,71 10,95 
Fe [iso оз» [nas Гола о 095 127 
210 | 0,37 | 0,44 | 0,52 | 0,59 | 0,67 | 0,74 | 0,81 | 0,89 | 0,96 | 1,04 | 1,11 | 1,48 


В отличие от существующих СВЧ-радиоме- 
тров [4], в представленном решении минимизиро- 
ваны потери в антенном канале. Это достигается 
за счет переноса направленного ответвителя в ка- 
нал второго опорного генератора шума. Уменьше- 
ние потерь в антенном канале повышает чувстви- 
тельность в предложенном СВЧ-радиометре. 


Использование двух параллельно работающих 
радиометрических приемников позволяет на задан- 
ном временном интервале накапливать в два ра- 
за больше результатов измерений. В соответствии 
с приведенными аналитическими соотношениями 
это уменьшает рабочую полосу частот в два раза 
при заданных длительности измерения и флуктуа- 
ЦИОННОЙ чувствительности. 


Разработанная математическая модель оцен- 
ки абсолютной точности измерений позволяет на 
этапе проектирования СВЧ-радиометра оценить 
уровень погрешности измерений при известных 
параметрах точности термостатирования, потерь 
во входной высокочастотной части и температур- 
ных коэффициентах опорных генераторов шума. 

Результаты численного моделирования показы- 
вают, что при температурных коэффициентах пер- 
вого и второго опорных генераторов шума 0,77 
и 069 К/°С соответственно и потерях во входной 


высокочастотной части юмшу, Пи» Не превышаю- 
щих суммарно уровень 0,2, величина максималь- 
ной погрешности измерений не превышает точно- 
сти термостатирования. 

Результаты исследований влияния изменяюще- 
гося уровня температурных коэффициентов опор- 
ных генераторов шума показывают, что для раз- 
работанного СВЧ-радиометра наибольшее влияние 
оказывает температурный коэффициент ОМШУ. 

Исследование распределения погрешности из- 
мерений по динамическому диапазону, составля- 
ющему 0-300 К, показало, что величина измене- 
ния погрешности измерений не превышает двух раз 
при равных величинах температурных коэффици- 
ентов опорных генераторов шума, величине потерь 
во входной высокочастотной части 0,1 и величине 
погрешности поддержания температуры термостата 
не более +1 К. 


Заключение 


Разработана схема приемного канала СВЧ-ги- 
перспектрометра, использующая два радиометри- 
ческих приемника и два опорных полупроводнико- 
вых генератора шума. 
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В качестве опорных генераторов шума приме- 
нены активные полупроводниковые генераторы шу- 
ма — лавинно-пролетный диод в режиме обратного 
пробоя и малошумящий усилитель, подключенный 
входом к входу радиометрического приемника. 

Подключение малошумящего усилителя вхо- 
дом к входу радиометрического приемника позволя- 
ет формировать на входе радиометрического прием- 
ника низкотемпературный шум, пропорциональный 
собственным шумам малошумящего усилителя. 

Использование двух опорных генераторов шу- 
ма позволяет минимизировать влияние дестабили- 
зирующих факторов радиометрических измерений, 
а именно дрейфа собственных шумов и коэффици- 
ента передачи радиометрического приемника. Это 
позволяет упростить калибровку в гиперспектраль- 
ном режиме и повысить устойчивость к изменяю- 
щимся факторам внешней среды, в первую очередь 
температуры. 

Перенос направленного ответвителя из антен- 
ного канала уменьшает потери на входе системы 
и приводит к повышению чувствительности. 

Использование двух приемников позволяет уве- 
личить разрешающую способность по частоте при 
проведении СВЧ-гиперспектральных исследований 
в два раза за счет пропорционального увеличения 
времени наблюдения за сигналом антенны по срав- 
нению с приемными каналами, выполненными на 
основе классических СВЧ-радиометров с одним ра- 
диометрическим приемником. 

Разработанная математическая модель оценки 
погрешности радиометрических измерений позво- 
лит на этапе проектирования оценить точность из- 
мерений в зависимости от технических параметров 
входной высокочастотной части. 

Результаты численного эксперимента на разра- 
ботанной модели показали, что влияние точности 
термостатирования на результаты измерений сни- 
жено по сравнению СВЧ-радиометрами, использу- 
ющими пассивные опорные генераторы шума. 

В разработанном СВЧ-радиометре основной 
вклад в погрешность измерений обусловлен тем- 
пературной нестабильностью уровня собственных 
шумов опорного малошумящего усилителя. 

Проведенные в статье результаты являются ос- 
новой для последующих практических эксперимен- 


тальных исследований технических характеристик 
разрабатываемого двухприемникового СВЧ-гипер- 
спектрометра для дистанционного зондирования 
природных сред. 


Исследование выполнено за счет гранта Рос- 
сийского научного фонда (проект Л№ 18-79-00045). 
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Особенности формирования 
информационно-вычислительной системы 
для решения измерительных задач 
в условиях различной степени структуризации 


В. К. Ларин, к. m. n., c. н. с., contact(Qspacecorp.ru 
АО «Российские космические системы», Москва, Российская Федерация 


Аннотация. В статье рассмотрены вопросы формирования информационно-вычислительной системы (ИВС) для решения из- 
мерительных задач в условиях различной степени структуризации. Решение задачи основано на принципах структурного 
анализа, центральным из которых является разделение сложной системы на простые составляющие с определением функ- 
циональных связей между ними. Вид решения может быть представлен как в графической, так и в аналитической формах. 
В работе сформулирована последовательность этапов решения поставленной задачи. 

Приводится общая структурная схема ИВС для решения технических задач. В соответствии с указанными принципами 
решение разбивается на методическую и технологическую части. Даны структурные схемы этих частей. Выделен фрагмент 
решения измерительной задачи, отражающий специфику перехода решения от предыдущего блока к последующему в зависимо- 
сти от выполнения «условий перехода». Показана зависимость «условий перехода» от значений информационных параметров. 
На примере задачи определения параметров движения КА рассмотрено применение сформированной ИВС для ее решения. 
Итогом исследований является структурная схема универсальной ИВС для решения технических задач различного типа. 


Ключевые слова: информационно-вычислительная система, измерительная задача, структура, условия перехода, информаци- 
онные параметры 
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ior Solving Measuring Problems 
in Conditions ої Varying Degrees ої Structuring 
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Abstract. The article considers the creation of an information computer system (ICS) for solving measuring problems in conditions 
of different degrees of structuring. The solution is based on the principles of structural analysis. The key principle is the division 
of a complex system into simple constituents and the delinition of functional relationships between them. The solution type can be 
given in both graphical and analytical form. The paper formulates a sequence of stages for solving the problem. 

А general block diagram of the ICS Гог solving technical problems is given. In accordance with the aforementioned principles, 
the solution is divided into a methodological and a technological part. Block diagrams of these parts are presented. А fragment 
of the solution to the measuring problem, which reflects the particular nature of solution transition from the previous block 
to the next depending on the fulfillment of the *transition conditions", is highlighted. The dependence of the "transition conditions" 
on the values of information parameters is shown. The problem of determining the motion parameters of a spacecralt is considered 
as an example ої implementing the developed ICS for its solution. 


Keywords: information computer system, measuring task, structure, transition conditions, information parameters 


52 В.К. ЛАРИН 


Введение 


Объектом исследования является некоторая 
измерительная задача, в которой определены вход- 
ные данные, вид результата и подход к решению 
данной задачи. 

Измерительная задача — это задача обра- 
ботки измерений с целью определения физических 
свойств (или характеристик) измеряемого объек- 
qa [1]. 

Примером может служить задача определения 
параметров движения КА, где объектом является 
орбита КА. 

В основу исследований положены принципы 
структурного анализа [2]. 

Решение любой инженерной задачи можно раз- 
делить на две части: методическую, состоящую из 
постановки задачи и алгоритма, и технологическую, 
включающую программу, информационно-вычисли- 
тельную систему (ИВС) и процедуру реализации 
программы на вычислительных средствах ИВС. 

Согласно ГОСТ Р 53622-2009 термин инфор- 
мационно-вычислительная система означает со- 
вокупность данных (или баз данных), систем управ- 
ления базами данных и прикладных программ, 
функционирующих на вычислительных средствах 
как единое целое для решения определенных за- 
дач [3]. 

Решение задачи зависит от степени ее струк- 
туризации. Согласно классификации, предложен- 
ной в [4], проблемы (к ним относятся и измери- 
тельные задачи) можно разделить по степени струк- 
туризации на три типа: структурированные (C3) — 
имеющие численное решение; неструктурирован- 
ные (HC3) — не имеющие решения; и слабострук- 
турированные (ССЗ), имеющие на фоне формаль- 


ного (численного) решения качественное решение, 
а в пределе — не имеющие решения. 


Решение задач различной степени 
структуризации 


На основе вышесказанного можно предпо- 
ложить, что слабоструктурированные задачи за- 
нимают пространство между структурированными 
и неструктурированными задачами. Очевидно, что 
степень структуризации задач на границе масси- 
ва структурированных задач будет максимальной 
(примем равной 1) и минимальной — на границе 
массива неструктурированных задач (примем рав- 
ной 0). 

С большой долей вероятности можно предполо- 
жить, что задачи со значением степени структури- 
зации, близким к 1, должны иметь количественное 
или качественное решение, задачи со значением сте- 
пени структуризации, близким к 0, должны иметь 
качественное решение, а в пределе — не иметь его. 

Для достаточно большого диапазона, в кото- 
рый попадают ССЗ, и неопределенной корректно- 
сти их решения, в [5] было предложено разделить 
массив слабоструктурированных задач на две ча- 
сти: задачи, примыкающие к массиву структуриро- 
ванных задач, названные квазиструктурированными 
задачами (КСЗ); и собственно слабоструктурирован- 
ные задачи. 

С учетом вышесказанного, а также приведен- 
ных выше значений степени структуризации (o60- 
значим через К) массив КСЗ приближенно характе- 
ризуется диапазоном К от 1,0 до 0,5, массив ССЗ — 
соответственно диапазоном К от 0,5 до 0. В табл. 1 
приведены в общем виде данные взаимосвязи вида 
решения задачи и степени структуризации. 


Таблица 1. Соответствие вида решения задачи и степени структуризации 


Аналитические методы 


Численное решение 


Аналитические методы 


Численное решение с отклонением 
точности результата 


Эвристические методы 


Качественное решение 
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Практически все измерительные задачи изна- 
чально имеют численное решение. При реализа- 
ции задачи на ИВС возможно получение результата 
с различной степенью точности. Это объясняется 
ошибками, допущенными при разработке алгорит- 
ма и программы, а также при формировании отдель- 
ных частей ИВС. Примем, что решение, получен- 
ное с ошибкой больше заданного предела, является 
непригодным для дальнейшего рассмотрения. 

Таким образом, цель настоящих исследований 
заключается в формировании ИВС, которая позво- 
ляла бы получить оптимальный результат для из- 
мерительных задач квазиструктурированного типа 
в виде численного решения, так как качественное 
решение для измерительных задач не имеет прак- 
тического значения. 

В статье рассмотрена технология построения 
ИВС для решения данной проблемы. 

Решение поставленной задачи заключается 
в последовательном выполнении следующих этапов. 

Этап 1. Определение структуры ИВС. На пер- 
вом этапе определяется структура ИВС, исполь- 
зуемая для решения информационных задач об- 
щего типа. Состав такой ИВС определен в ГОСТ 
Р 53622-2009. 

Этап 2. Определение методической и техно- 
логической частей решения измерительных задач. 
На этом этапе определяется структура методической 
и технологической частей решения измерительных 
задач, проводится анализ структур с целью выделе- 
ния блоков, несущих наибольшую нагрузку некор- 
ректности в элементах, и определяется влияние 
изменения характеристик элементов некорректных 
блоков на точность решения задачи в заданной 
предметной области. 

Этап 3. Формирование основной структуры 
ИВС. На третьем этапе формируется основная 
структура ИВС с блокировкой решений, получен- 
ных с ошибками, превышающими предельные зна- 
чения. 


Этап 1. Определение структуры ИВС 


В соответствии с требованиям ГОСТ Р 53622- 
2009 ИВС включает следующие виды обеспече- 
ния: информационное, программное, алгоритмиче- 
ское и техническое. С учетом сказанного струк- 


Специальное 
IIO 


Результат 


Рис. 1. Структурная схема ИВС 


турная схема ИВС общего типа для решения тех- 
нических задач в технологическом варианте будет 
иметь следующий вид (рис. 1). 

На рис. 1 введены сокращения: ИО — инфор- 
мационное обеспечение, ПО — программное обес- 
печение, ТО — техническое обеспечение, БД — 6a- 
за данных, СУБД — система управления БД. 

Алгоритмическое обеспечение не включено 
в схему, так как оно отнесено к методической части 
и рассматривается отдельно. 


Этап 2. Определение методической 
и технологической частей решения 
измерительных задач 


Специфика измерительных задач (ИЗ) харак- 
теризуется следующими свойствами: 

— большим объемом входных данных (измере- 
ний), ошибки которых имеют в основном случай- 
ный характер; 

— используемыми методами обработки — филь- 
трация измерений по критериям с последующей об- 
работкой статистическими методами (метод наи- 
меньших квадратов, метод Монте-Карло); 

— строгим соблюдением правил формирования 
и решения условных и нормальных уравнений по- 
правок к измерениям (соответствие числа уравне- 
ний числу неизвестных); 

— заданием оптимальных значений пределов 
внутренней сходимости в процедурах решения 
с использованием метода последовательных при- 
ближений (согласование быстродействия и точно- 
сти решения). 
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Указанные свойства измерительных задач зна- 
чительно усложняют поиск решения. Во-первых, 
необходимостью формализации свойств для их уче- 
та и, во-вторых, определением возможного диа- 
пазона их значений. Для решения этих вопросов 
целесообразно использовать методы структурного 
анализа, а также подход, используемый в [6], где 
рекомендуется предварительно разделить процеду- 
ру решения на методическую и технологическую 
части, представив каждую часть в виде струк- 
турной схемы. Структурные схемы методической 
и технологической частей решения представлены 
на рис. 2 и 3 соответственно. 


Сам по себе термин «слабоструктурированная 
проблема (задача)» говорит о «слабых» структур- 
ных связях частей системы или об их отсутствии. 
Аналитически связи можно описать с помощью 
условий перехода между соседними блоками про- 
граммы. По аналогии со структурными связями 
условия могут отсутствовать или быть не точно 
сформулированными. В таком случае либо может 
произойти сбой решения, либо получен отличаю- 
щийся по точности результат. 

Возможны следующие варианты 
некотором шаге программы: 

— остановка расчета вследствие невыполнения 
условия перехода к очередному блоку программы 


«сбоя» на 


Постановка (команда «Останов»); 
задачи 
а Входные данные | - получение неверных данных в очередном 
блоке программы вследствие влияния критического 
Метод решения | 
значения информационного параметра (ИП) и пере- 
Решение | ход к следующему блоку программы на очередном 
шаге за счет изменения условий перехода. 
Алгоритм Структурная схема фрагмента «сбойного» ва- 
Исходные данные | 
рианта решения при возможном нарушении усло- 
| Фильтрация | вий перехода между соседними блоками програм- 
измерений 
P мы представлен на рис. 4. 
Статистические На рис. 4: U(k,) — номинальное условие пе- 
методы 
рехода от блока n к блоку n + 1; 1, — ИП блока m; 
Методы U(k,,1,) — измененное условие перехода от бло- 
математического 
анализа ка п к блоку n + 1; № — номинальное значение 
Результат 
ешения 
E 1. Решение 
Рис. 2. Структура методической части решения | 
2. Блок п 
"m | Ens kn 
Системная 
часть 3. U(k = 4. 1 
Программа Us.) Her ii 
Функциональная 
часть Да 
Y 
Т. Блок n4-1 
Информационное EN CYB UL р k 5. U (knla) 
ИВС обеспечение Д, Д yep В 
р Y 
Архивы | | Ta No. N «UI « N, 
d Her 
Техническое ЛВС, сервер | 8. Результат - 
обеспечение Останов 


ЛВС — локальная вычислительная сеть 


Рис. 3. Структура технологической части решения 


Рис. 4. Схема фрагмента перехода расчетов между со- 
седними блоками программы m, п + 1 
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условия перехода; N,—N, — диапазон значений 
условия перехода для конкретного блока; Ul = 
= U (kn, In) (в блоке 6); Fa, №, — функциональ- 
ное содержание блоков п, n + 1. 

В алгоритмической форме этот фрагмент мо- 
жет быть представлен в следующем виде: 

1. Е, — расчет блока п. 

2U (k) =N > Fnr A Ip 

3. U (kn, In) = Ul — коррекция условия перехо- 
да (УП) как функции информационного параметра. 

4. N SUIL < №, > т >-› Останов, 
где (7, — номинал условия перехода блока п; 

I, — информационный параметр блока m; 

k, — критерий расчета Fp; 

К: — критерий расчета F5 

ô — поправка к УП. 

Для анализа условий перехода представим про- 
грамму в виде набора замкнутых математических 
процедур, где переход между ними осуществляется 
по выполнению условий. 

Вид условия зависит от содержания очередной 
процедуры (блока). 

Для измерительных задач возможны следую- 
щие условия перехода: 

— полученные данные на очередном шаге долж- 
ны быть больше, меньше или равны некоторому зна- 
чению; 

— полученные данные должны находиться в за- 
данном диапазоне значений; 

— единицы измерений полученных данных 
должны соответствовать единицам измерений ве- 
личин, используемых на следующем шаге реше- 
ния (проверка на идентичность единиц измерения 
должна проходить на каждом шаге программы). 

В качестве примера рассмотрим задачу опреде- 
ления параметров движения космического аппара- 
та (ОПД КА), условиями перехода в которой могут 
служить следующие соотношения рассчитываемых 
параметров: 

— при формировании двойных разностей из- 
мерений фазы длина базовой линии должна нахо- 
диться в заданном диапазоне; 

— при фильтрации измерений (кода или фазы) 
значение сигнал/шум каждого измерения не дол- 
жен быть меньше заданной величины; 

— при фильтрации измерений (кода или фазы) 
по углу места, при котором производилось измере- 


ние, его значение не должно быть меньше задан- 
ного; 

— при формировании сеанса количество изме- 
рений не должно быть меньше заданного; 

— при формировании математической модели 
движения (ММД) КА учитываемые возмущающие 
ускорения должны обеспечивать заданную точ- 
ность параметров орбиты КА; 

- при формировании математической моде- 
ли измерений (ММИ) учитываемые составляю- 
щие должны соответствовать заданному набору по- 
правок; 

— при оценке решения «краевой задачи» по дан- 
ным определения местоположения отличия полу- 
ченных значений координат от расчетных не долж- 
ны превышать заданной величины. 

В соответствии c [7] ИП зависят от значе- 
ний функций и процедур, составляющих основу ре- 
шения, и могут совпадать с условиями перехода 
и влиять на точность результата. Ниже представ- 
лен перечень основных ИП ОПД КА. 

Как и в предыдущем случае, ИП целесообраз- 
но разделить по их влиянию на две части: мето- 
дическую (постановка и алгоритм) и технологиче- 
скую (информационное обеспечение и программа). 

ИП элементов методической части решения 
и их критические значения 

Постановка задачи 

1. Цель решения (неоднозначно сформулиро- 
вана). 

Алгоритм 

1. Количество измерений в сеансе (меньше 
нормы). 

2. Учет возмущающих сил в модели движения 
(недостаточно). 

3. Учет поправок в модели измерений (недо- 
статочно). 

4. Критерии фильтрации измерений (не соот- 
ветствуют нормам). 

5. Диапазон базовых линий (задан неверно). 

6. Величина разрядки измерений (задана не- 
верно). 

ИП элементов технологической части реше- 
ния и их критические значения 

Информационное обеспечение 

1. Координаты измерительных станций (зада- 
ны неверно). 
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Таблица 2. Перечень ИП и их значения 


Наименование Символ 


Соотношение 


TEK x Примечание 


и номинального 
значений ИП 


Методическая часть 


Возмущающие силы в ММД Недостаток учета возмущающих сил 
Поправки в модель измерений ôI, = АГ, Недостаток учета поправок в ММИ 
Критерии фильтрации Е; + Fy Критерий не соответствует норме 
Диапазон базовых линий AL, + № Базовые линии вне диапазона 
Величина разрядки измерений At, X Aty Разрядка He соответствует норме 


Технологическая часть 


Измерительные станции 


Таблицы 


Данные в ячейках имеют ошибки 


Настройки 


Настройки не равны заданному 
значению 


Выбор исходных данных из архива 


2. Формирование таблиц (с ошибками). 

3. Выбор данных из таблиц (откат транзак- 
ЦИЙ). 

Программа 

1. Настройки (с ошибками). 


2. Выбор исходных данных из архивов (с ошиб- 
ками). 


З. Циклы (с ошибками). 
4. Условия (с ошибками). 


5. Выбор и использование стандартных под- 
программ (с ошибками). 


6. Формирование результата (с ошибками). 


ИД с ошибками 


а Га EID 


Несоответствие норме выбранной 
Стандартная программа Sp Æ Sp, 
а a программы 


Формирование результата DEMO | Формирование неверного результата | неверного результата 


В табл. 2 приведен перечень ИП и их значе- 
ний, приводящих решение к одному из сбойных ва- 
риантов. 

Как показывает опыт работы, два случая могут 
привести к невозможности получить решение: 

— неопределенность цели решения, что соот- 
ветствует ИП в постановке задачи; 

— ошибки при настройке программы, что соот- 
ветствует ИП в программе. 

Все остальные случаи приводят к неточному 
результату. 

Результаты отбора ИП второго случая, в соот- 
ветствии с [8], представлены в табл. 3. 
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Таблица 3. Результаты отбора ИП 


№ 


Название ИП 
п/п 


Критическое условие ИП 


Примечание 


3 | Исходные данные из архивов 


Архив 


2 | Координаты измерительных станций В.Г, = Вр, [и 


Тр # І, 


Базовые линии вне диапазона 
Координаты ИС заданы неверно 


ИД заданы неверно 


Возмущающие силы в ММД 


БД, СУБД 


5 | Выбор данных из таблиц БД 


Программа 


Ошибки в ячейках 


Откат транзакций 


Недостаток учета возмущающих сил 


6 
7 | Поправки в модель измерений 
8 


9 | Величина разрядки измерений At, X Aty 


Недостаток учета поправок 
Критерий не соответствует норме 


Разрядка не соответствуют норме 


10 | Циклы Параметры цикла заданы неверно 
11 | Логические условия If, Z If, Условие задано неверно 
12 | Стандартные подпрограммы Sp = Sp, Несоответствие норме выбранной 


стандартной подпрограммы 


Этап 3. Формирование основной 
структуры ИВС 


С учетом схемы фрагмента программы (рис. 4) 
и данных табл. 3 можно сформировать ИВС для из- 
мерительной задачи определения параметров движе- 
ния КА, учитывающую степень ее структуризации. 
Структурная схема ИВС представлена на рис. 5. 

ИПА 1, 2, 3 — информационные параметры 
блока «Архив»; ИПД 1, 2 — информационные па- 
раметры блока БД, СУБД; блок 1,2,...,n—1,n — 
блоки программы; фрагмент 1,2,...,n—1,n — nog- 
программа перехода между блоками, индекс соот- 
ветствует номеру блока (рис. 4); о, — среднеквад- 
ратическая ошибка результата; № — предельное 
значение ошибки результата. 

ИПА и ИПД могут оказать влияние на реше- 
ние задачи при выполнении команд считывания ин- 
формации. 

Безусловный переход от блока через фраг- 
мент Å к очередному блоку обеспечивается скор- 
ректированным «условием перехода» за счет ИП. 


Таким образом, на рис. 5 представлена ИВС, 
сформированная с учетом возможной некоррект- 
ности отдельных элементов составных частей, вы- 
раженных в формате ИП. Конкретное содержание 
блоков и фрагментов ИВС должно соответствовать 
функциональному назначению задачи. 

Из материалов статьи следует, что точность 
решения задачи зависит от величины отклонения 
условий перехода от нормы, что аналогично степе- 
ни структуризации составных частей ИВС. 


Заключение 


1. Разработана технология формирования 
ИВС для решения измерительных задач с различ- 
ной степенью структуризации. 

2. В основу решения проблемы положены прин- 
ципы структурного анализа, позволяющие значи- 
тельно упростить процесс решения подобных задач. 

3. В качестве примера рассмотрена задача опре- 
деления параметров движения КА, где показано 
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Рис. 5. Структурная схема ИВС 


влияние «условий перехода» на точность решения 
задачи. 

4. Предложенная структура ИВС имеет уни- 
версальный характер и может использоваться для 
решения измерительных задач различного типа. 
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Введение 


Координационная группа по метеорологическим 
спутникам CGMS была образована в 1972 году. 
Росгидромет присоединился к CGMS в 1973 году. 
В настоящее время членами CGMS также явля- 
ются Госкорпорация по космической деятельности 
«Роскосмос», Китайская метеорологическая адми- 
нистрация (СМА), Китайское национальное кос- 
мическое управление (CNSA), Европейская органи- 
зация по эксплуатации метеорологических спутни- 
ков (EUMETSAT), Европейское космическое агент- 
ство (ЕЅА), Французское космическое агентство 
(CNES), Индийская организация космических ис- 
следований (ISRO), Индийский метеорологический 
департамент (IMD), Японское агентство аэрокос- 
мических исследований (JAXA), Японское метео- 
рологическое агентство (JMA), Корейская метео- 
рологическая администрация (КМА), Националь- 
ное управление по аэронавтике и исследованию 
космического пространства США (NASA), Наци- 
ональное управление океанических и атмосферных 
исследований США (NOAA), Всемирная метеоро- 
логическая организация (WMO), Межправитель- 
ственная океанографическая комиссия ЮНЕСКО 
(IOS-UNESCO) 

Кроме того, B состав CGMS со статусом на- 
блюдателей входят Канадское космическое агент- 
ство (CSA), Глобальная система наблюдений за кли- 
матом (GCOS), Корейский институт аэрокосмиче- 
ских разработок (КАВТ), Корейский институт оке- 
анических исследований и развития (KORDI), Го- 
сударственное океанологическое управление Китая 
(SOA). 

Международная рабочая группа по Глобаль- 
ной системе интеркалибровки спутниковых данных 
GSICS была образована совместно ВМО и CGMS 
в 2005 году с целью мониторинга и улучшения каче- 
ства наблюдений с оперативных метеорологических 
спутников Глобальной системы наблюдений ВМО. 

Росгидромет и Роскосмос имеют своих пред- 
ставителей в CGMS и GSICS, используют участие 
в работе этих международных мероприятий для вы- 
работки общих подходов к формированию данных 
космической съемки, более эффективного проведе- 
ния работ по калибровке/валидации информацион- 
ных данных отечественных метеоспутников и совер- 


шенствования аппаратуры с целью повышения точ- 
ности и стабильности радиометрических измерений, 
а также точности их геопривязки. 


Основные принципы организации 
работ метеорологических агентств, 
космических агентств и операторов 
КА ДЗЗ по калибровке/валидации 
спутниковых данных 


Спутниковые данные для мониторинга климата 
и других задач гидрометеорологии требуют точной 
калибровки, согласованной для различных однотип- 
ных приборов и обеспечивающей в течение длитель- 
ного времени штатной эксплуатации КА хорошее 
соответствие с данными референсного (эталонного) 
прибора. Работы в рамках Международной системы 
калибровки и валидации спутниковых приборов, ко- 
ординируемых CGMS и GSICS, имеют важное 3Ha- 
чение как для исследования климатических процес- 
сов, так и для совершенствования организационных, 
аппаратурных, технологических и методических ас- 
пектов всех видов калибровок. 


Основная идея интеркалибровки заключается 
в том, что два КА проводят съемку одной и той же 
территории в одно и то же время, с одинаковыми 
пространственными характеристиками и геометри- 
ей обзора. В результате интеркалибровки нужно 
количественно определить относительную погреш- 
ность измерений между проверяемым и эталонным 
приборами для конкретной тандемной съемки, най- 
ти причины несоответствия данных, скорректиро- 
вать или устранить их. 


Как показано в работе [1] на примере аппа- 
ратуры МСУ-МР, интеркалибровка двух приборов 
разных космических аппаратов требует точного со- 
гласования в 8-мерном пространстве, включающем 
время, широту, долготу, высоту, солнечный зенит- 
ный угол, угол обзора по азимуту, зенитный угол 
наблюдения и спектральный диапазон. Значитель- 
ные погрешности в любой из перечисленных со- 
ставляющих могут привести к случайным и систе- 
матическим ошибкам, что впоследствии приведет 
к неточностям измерения. 


РАКЕТНО-КОСМИЧЕСКОЕ ПРИБОРОСТРОЕНИЕ И ИНФОРМАЦИОННЫЕ СИСТЕМЫ т. 7 вып. 1 2020 


62 IO. M. ГЕКТИН, C. M. 3OPHH, Д.О. ТРОФИМОВ, И. A. БАРСУКОВ, К. И. ЖУКОВСКАЯ 


Выбор эталонного прибора также вызыва- 
ет сложности — ключевой его характеристикой 
должна являться радиометрическая стабильность, 
но нужно учитывать и другие важные факторы 
(орбиты спутников должны позволять проводить 
съемку одних и тех же территорий, приборы долж- 
ны иметь схожие спектральные характеристики 
и др.). Кроме того, эталонный прибор выбирается 
путем проверки качества и признания его таковым 
странами — членами международного сообщества 
в рамках CGMS и GSICS с учетом малых погреш- 
ностей радиометрических измерений, долговремен- 
ной стабильности измерительных характеристик, 
а также условий и точности наземной радиомет- 
рической калибровки и наличия высокостабильной 
системы бортовой калибровки. 

Например, как показано в [1,2], для аппа- 
ратуры МСУ-МР и МСУ-ГС видимого диапазо- 
на в качестве опорных радиометров использова- 
лись соответственно AVHRR (КА MetOp, Espo- 
союз) и VIIRS (KA Suomi NPP, США). 

А для интеркалибровки спутниковой ИК-аппа- 
ратуры в настоящее время в рамках GSICS в ка- 
честве основного эталонного прибора используется 
интерферометр IASI (КА MetOp, Евросоюз). Кроме 
того, для интеркалибровки могут использоваться 
данные инфракрасного (ИК) зондировщика AIRS 
(KA Aqua, США). 

Помимо IASI и AIRS в качестве промежуточ- 
ных эталонов для ИК-диапазона спектра использу- 
ется ряд приборов, стабильность и точность радио- 
метрических характеристик которых была подтвер- 
ждена в ходе различных интеркалибровок целевой 
аппаратуры (ЦА): 

- ИК-зондировщик | CrIS/SNPP (работы по 
проверке качества абсолютной калибровки завер- 
шились в 2013 г.); 

- сканирующий радиометр MODIS/EOS Terra, 
Aqua; 

- сканирующий радиометр AVHRR/NOAA; 

– сканирующий радиометр VIIRS/SNPP; 

- сканирующий радиометр SEVIRI/Meteosat 10. 

Данные измерений перечисленных приборов 
(регионального и глобального покрытия) имеют- 
ся в открытом доступе или могут быть получены 
из спутниковых архивов, поддерживаемых МОАА, 
NASA, EUMETSAT и используются для интерка- 


либровки данных космической съемки ЦА, уста- 
новленной на КА серий «Mereop-M», «Электро-Л», 
«Канопус-В-ИК» [1-3]. 

Опыт разработки и эксплуатации целевой аппа- 
ратуры микроволнового диапазона (например, мик- 
роволновых сканеров/зондировщиков температур- 
но-влажностного зондирования атмосферы МТВЗА 
и МТВЗА-ГЯ для KA «Mereop-3M» №1 и КА ce- 
рии «Mereop-M») подтверждает необходимость со- 
здания процедур интеркалибровки и внешней ка- 
либровки спутниковых данных. В настоящее вре- 
мя микроволновый сканер/зондировщик МТВЗА-ГЯ 
является одним из основных видов целевой аппара- 
туры на метеорологических спутниках серии «Ме- 
Teop-M». По информационным характеристикам 
и области применения МТВЗА-ГЯ соответствует 
зарубежным аналогам — зондировщикам AMSU-A 
и AMSU-B (КА NOAA, США и MetOp, Espo- 
союз), ATMS (KA Suomi NPP, США) [4], ска- 
нерам  AMSR-E (KA Aqua США), AMSR2 
(KA GCOM-WI, Япония) [4] и сканеру/зондиров- 
щику SSMIS (KA DMSP Е16-Е19, США). 


Для внешней калибровки микроволновых ска- 
неров/зондировщиков выполняются сопоставления 
измеренной антенной температуры (отсчетов) и ра- 
счетной яркостной температуры над протяженны- 
ми однородными «холодными» и «горячими» участ- 
ками океана и земной поверхности [4, 5]. В каче- 
стве «горячей» области используются широколист- 
венные леса бассейна Амазонки, при этом выбира- 
ются измерения Hà витках, когда над участком со 
сплошным лесным покровом отсутствует мощная 
облачность и осадки. «Холодные» области выбира- 
ются над отдельными районами Антарктиды и над 
районами океана с температурой поверхности воды 
от 0 до +12 °С, слабым ветром, отсутствием обла- 
ков и низкими значениями паросодержания атмо- 
сферы [4]. 

Для интеркалибровки пассивных измерений ми- 
кроволновых приборов типа сканера/зондировщика 
МТВЗА-ГЯ в качестве зарубежных референсных 
приборов можно использовать прежде всего аппара- 
туру с коническим режимом сканирования — ска- 
неры/зондировщики SSMIS (КА DMSP Е16-Е19, 
США), AMSR-E/EOS (KA Aqua, США), AMSR2 
(KA GCOM-WI, Япония). При отсутствии доступа 
к информации указанных приборов для интерка- 
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либровки данных в отдельных каналах МТВЗА-ГЯ 
с определенными ограничениями можно использо- 
вать измерения микроволновых зондировщиков 
АМЅ0-А и AMSU-B (КА NOAA, США и MetOp, 
Евросоюз) и АТМ$ (KA Suomi NPP и JWSS, 
США) с поперечным режимом сканирования. 


Краткий обзор состояния 
зарубежных спутниковых 
группировок и целевой аппаратуры 
гидрометеорологического 
назначения по результатам 
докладов на заседаниях ССМ$-47 

и GSICS-EP-20 (г. Сочи, Россия, 
16-24 мая 2019 г.) 


Лидирующее положение в вопросах координа- 
ции и организации международного сотрудничества 
в вопросах глобальных спутниковых систем наблю- 
дения за метеорологическими (и смежными с ними) 
параметрами в целях анализа и прогноза состояния 


окружающей среды, в том числе опасных явлений 
погоды и их заблаговременного предупреждения, за- 
нимают США, Евросоюз, Япония и Китай. 
Соответственно важен анализ состояния спут- 
никовых группировок и ЦА гидрометеорологическо- 
го назначения в первую очередь указанных стран. 


США (NASA, NOAA) 


На рис. 1, 2 представлен состав действующей 
и перспективной группировок гидрометеорологиче- 
ских спутников США из доклада NASA [6]. 

По состоянию на май 2019 г. США имеют 
23 действующих гидрометеорологических спутни- 
ка (NASA), без учета геостационарной космиче- 
ской системы GOES (NOAA). С момента прошлого 
заседания CGMS-46 было запущено 10 новых MHC- 
сий, включающих: ICESat-2, GRACE-FO, OSO-3, 
ECOSTRESS, GEDI и пять ультракомпактных КА 
(типа CubeSat): RainCube, Cube ВВТ, CSIM-FD, 
TEMPTEST-D, MinXSS-2. B период 2019-2025 rr. 
NASA планирует запустить еще 15 миссий, 
в их числе TEMPO, Landsat-9, 10, Sentinel-6A/B, 
SWOT, PACE, CLARREO-PF и др. (рис. 1). Таким 


NASA: Спутниковые системы ДЗЗ 


Миссии: от настоящего времени до 2023 г. 


О 


Шпродленные 


Аппаратура на МКС 
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BOE INOAA] 
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Рис. 1. Группировка гидрометеорологических спутников ДЗЗ США [6] 
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Обзор текущих и разрабатываемых спуниковых систем ДЗЗ Л 


Система 
Landsat-7 1999 
Terra 1999 
Aqua 
SORCE 
Aura 
CALIPSO 
CloudSat 
Jason-2 
Suomi-NPP 2011 
Landsat-8 
TCTE 
GPM Core 
OCO-2 


Текущее количество KA — c 
(исключая CubeSa 
(на апреля 2019 г.) 


просуществуют дольше 


DSCOVR 
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CYGNSS 
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LIS-ISS 
TSIS-1-ISS 
GRACE-FO 
ECOSTRESS-ISS 
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(исключая CubeSat) 


К 2020 г. будут запущены 3 дополнительных КА и 2 прибора 


Стандартные спутниковые системы NASA рассчитаны на 3-5 лет 
функционирования, но этот срок, как правило, выше 


GEDL-ISS 


(не показаны в таблице) 
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Планируемые КА и приборы c запуском после 2020 г. 
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Рис. 2. Состав действующей и перспективной группировок спутников ДЗЗ США [6] 


образом, суммарный состав группировки гидроме- 
теорологических спутников США превысит 30 КА. 

При этом на данный момент в составе груп- 
пировки спутников ДЗЗ США есть значительное 
число КА со сроком службы 15-20 лет, которые 
продолжают успешно функционировать (рис. 2). 

Также увеличивается количественный и ка- 
чественный состав аппаратуры ДЗЗ США (LIS, 
SAGE Ш, Т515-1, OCO-3, ECOSTRESS, GEDI, 
CLARREO-PF, EMIT) на международной космиче- 
ской станции (МКС) (рис. 3) [7]. 

Необходимо отметить, что в 2017-2019 гг. бы- 
ли запущены и введены в эксплуатацию ультраком- 
пактные КА (типа CubeSat), такие как MiRaTA, 
HARP, RAVAN, IceCube, RainCube, TEMPEST-D, 
CubeRRT, CIRiS, которые подтвердили перспектив- 
ность развития данного направления ДЗЗ в США. 

Кроме того, в NASA и NOAA явно обозначе- 
ны тенденции по коммерциализации всех звеньев 
создания и функционирования КС от разработки 
техники до распространения данных ДЗЗ, форми- 
рования прогнозов погоды и климата. 

По состоянию на сентябрь 2018 г. были заклю- 
чены контракты в рамках пилотных коммерческих 
программ по созданию многоспутниковых систем 
на основе малых КА: 


- DigitalGlobe, 5 спутников наблюдения Земли 
с разрешением до 30 см (СеоЕуе-1, WorldView-1, 
WorldView-2, WorldView-3, WorldView-4); 

— Planet (ранее Planet Labs), три спутниковые 
системы (SkySat, Dove, RapidEye) с группировкой 
более 150 спутников наблюдения Земли со сред- 
ним и высоким разрешением и высокой частотой 
обновления данных; 

— SPIRE, более 60 спутников с приборами pa- 
диозатменного зондирования, определения коорди- 
нат летательных аппаратов и судов. 


Евросоюз (ESA, EUMETSAT) 


Спутниковая группировка Евросоюза по состо- 
янию на май 2019 г. содержит 15 действующих 
спутников, 25 находятся в разработке, 12 готовят- 
ся к запуску. За последнее время можно отме- 
тить успешный запуск и ввод в эксплуатацию КА 
Sentinel-5P (октябрь 2017 г.), Sentinel-3B (апрель 
2018 r.), Aeolus (август 2018 r.), MetOp-C (ноябрь 
2018 r.). 

Ha puc. 4 представлен план-график развития 
спутниковой группировки ЕЗА до 2040 г. [8]. Сле- 
дует отметить представленные согласованные планы 
NASA и ЕЗА (программа Copernicus) по развитию 
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Рис. 3. Состав аппаратуры ДЗЗ США на МКС [7] 


Планируемые спутниковые системы EUMETSAT 
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Рис. 4. Спутниковая группировка ДЗЗ Евросоюза [8] 


научной составляющей программ наблюдения Зем- 
ли из космоса, которые предполагают создание но- 
вых приборов, нацеленных на изучение качества 
атмосферного воздуха, состава атмосферы, опреде- 


ление содержания углекислого газа, озона, метана 
и других парниковых газов, а также на изучение 
интенсивности процессов фотосинтеза и роста рас- 
тений. 
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Перечисленные задачи требуют создания прин- 
ципиально новых приборов, таких как спектро- 
метры сверхвысокого спектрального разрешения 
FLEX (Fluorence Explorer), OCO (Orbiting Carbon 
Observatory), GeoCarb (Geostationary Carbon Cycle 
Observatory), a также ИК-радиометры и спектро- 
метры с чувствительностью в дальнем ИК- и мик- 
роволновом диапазонах, такие как CORSAIR (Cali- 
brated Observations of Radiance Spectra in the Far- 
Infrared), CLARREO (Climate Absolut Radiance and 
Refractivity Observatory), ECOSTRESS (ECOsys- 
tem Spaceborne Thermal Radiometer Experiment on 
Space Station), ICESat-2 (Ice, Cloud and land Eleva- 
tion Satellite 2), GEDI (Global Ecosystem Dynamics 
Investigation), а также ряд приборов в рамках про- 
граммы MeTop-SG — MWI, MWS, ICI и др. 


Япония (JAXA, JMA) 


Основные КА японской программы наблюде- 
ния Земли из космоса представлены на рис. 5 [9]. 
Японская аппаратура гидрометеорологического 
назначения, помимо приборов видимого, ближнего 
ИК- и ближнего УФ-диапазонов (SGLI), включает 
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также ИК-радиометры (IRS, VNR, AHI), ИК-фурье 
спектрометры (TANSO-FTS), приборы наблюдения 
за парниковыми газами (КА ШиК, 7 приборов) 
и микроволновые радиометры (GMI, АМ$В). Необ- 
ходимо отметить тесное сотрудничество японских 
агентств ЈАХА и ЈМА с соответствующими агент- 
ствами США и Евросоюза. На заседаниях ра- 
бочих групп GSICS и CGMS часто обсуждают- 
ся вопросы интеркалибровки ИК-радиометров AHI 
(Advanced Himawari Imager), которые изготовлены 
с участием США (аналоги ABI), установлены на 
КА Himavari-8, 9 и хорошо себя зарекомендова- 
ли при длительной эксплуатации, подтвердив сли- 
чения с данными, полученными с приборов IASI 
и CrIS/SNPP (США). На заседаниях WGII co сто- 
роны представителя США была выражена особая 
заинтересованность в данных микроволнового ска- 
нера с коническим сканированием и обеспечении 
запуска AMSR-3 в 2022-2023 гг. 


Китай (СМА, CNSA) 


На рис. 6 представлены история и стратегия 
запусков гидрометеорологических КА китайской 
программы наблюдения Земли из космоса [10]. 
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Рис. 5. Японская аппаратура гидрометеорологического назначения [9] 
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Координационная группа по метеорологическим спутникам — CGMS 


4.Будущие программы 


Государственная программа метерологических спутников Фэнъюнь (Fengyun) (2011-2020) 
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Рис. 6. История и стратегия запусков гидрометеорологических КА по программе КНР Фэнъюнь (Fengyun) [10] 


На данный момент на орбите работают 7 КА 
КНР (4 на геостационарной и 3 на низкой около- 
земной орбитах). Геостационарный КА РУ-2Н за- 
пущен в июне 2018 г. и принят в эксплуатацию 
в январе 2019. 

В составе ЦА китайской группировки пред- 
ставлен весь спектр приборов мониторинга атмо- 
сферы, океана, Солнца и Луны, МСУ, микровол- 
новые сканеры-зондировщики, ИК-фурье-спектро- 
метры (ИКФС), регистраторы молний (lightning 
таррег) и т. д. 

С 2019 по 2023 гг. планируется запустить 6 но- 
вых КА (FY-3E, FY-4B, ЕУ-4Е, FY-RM, FY-4C, 
FY-3G) (см. рис. 6) [10]. 

В докладе СМА на GSICS отмечалось, что 
КНР планирует масштабные проекты по изучению 
Луны, которые помогут получить также и необхо- 
димые данные для радиометрической калибровки 
аппаратуры ДЗЗ по лунному диску. 

В ряде докладов специалистов КНР сообщается 
о проекте разработки национального радиометриче- 


ского эталонного стенда космического базирования 
(Space-based Radiometric Benchmark) с целью прове- 
дения радиометрических калибровок ЦА в абсолют- 
ных единицах системы СИ (рис. 7) [10]. 

Основные исполнители проекта: 

— Академия оптоэлектроники (AOE), Китай- 
ская академия наук (CAS); 

— Национальный спутниковый метеорологиче- 
ский центр (NSMC), СМА; 

— Шанхайский институт технической физики 
(SITP), CAS; 

— Хэфэйский институт естественных наук 
(HIPS), CAS. 

Этап b (2018-2022 гг.) предполагает финанси- 
рование 300 млн юаней (около 3 млрд руб.) [10], 
что говорит о значимости данного направления 
в планах развития систем ДЗЗ. Данная програм- 
ма предусматривает такие исследования в об- 
ласти метрологии, как измерения характеристик 
АЧТ, коэффициентов излучения различных мате- 
риалов, создание бортовых АЧТ с применением 
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Проекты бортовых радиометрических эталонов в Китае 
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Рис. 7. Радиометрический эталонный стенд космического базирования КНР [10] 


калибровочных ячеек на фазовых переходах спла- 
вов галлия, воды и органических веществ, приме- 
нение лазерной техники. 

Данный подход позволит снизить неопределен- 
ность калибровки по температуре полостного бор- 
тового AYT до 10 мК для диапазона температур 
от 270 до 350 К (коэффициент излучения не менее 
0,997 и погрешность воспроизведения температуры 
не более 0,15 К). 


Заключение 


Таким образом, можно выделить следующие 
основные направления развития КС ДЗЗ гидроме- 
теорологического назначения. 

1. В части точных радиометрических изме- 
рений для прогноза погоды, мониторинга со- 
стояния климата, акваторий морей и океанов 
большое внимание уделяется приборам ИК-диа- 
пазона, таким как многозональные сканиру- 
ющие устройства (МСУ), ИК-фурье-спектроме- 
тры (ИКФС), а также микроволновые пассивные 


РАКЕТНО-КОСМИЧЕСКОЕ ПРИБОРОСТРОЕНИЕ И 


сканирующие радиометры. Основные требования 
к этим приборам состоят в долговременной стабиль- 
ности измерительных характеристик и минималь- 
но достижимых инструментальных погрешностях. 
Так, например, представленные результаты интер- 
калибровки Himawari-8/AHI (Япония), проведенной 
в 2019 г., при сличении с данными эталонных при- 
боров MetOp-A/IASI (США) и SNPP/CIIS (США) 
показали стандартное отклонение измеренных тем- 
ператур в разных ИК каналах от 0,005 до 0,039 К, 
а значения систематической погрешности от 0,01 
до 0,24 К за год. Эти данные подтверждают очень 
высокую стабильность и точность измерений ради- 
ационной температуры зарубежной ИК-аппаратуры. 

2. Признано, что алгоритмы GSICS и эталон- 
ные приборы GSICS (IASI, CrIS, L1 Microwave 
SSMIS, AMSR-E, AMSR-2, MSU/AMSU FCDR) 
могут быть использованы для мониторинга и pe- 
калибровки соответствующих российских приборов 
МСУ-ГС, MCY-MP, ИКФС-2, МТВЗА-ГЯ. 

З. В рамках GSICS в период до марта 
2019 г. проводились работы по интеркалибровке 
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следующей ИК аппаратуры: Meteosat-8-11/SEVIRI 
(EUMETSAT), MTSAT-2 Imager, Himawairi-8-9/ 
АНІ (JMA); GOES-11, -12, -15 Imager, GOES-16, 
-17/АВІ (NASA/NOAA); FY-2C,D,E,F,G/VISSR, 
FY-4A/AGRII/GIIRS (СМА); COMS, KOMP- 
ЅАТ-2А (КМА); INSAT-3D/Imager, INSAT-3D/ 
Sounder/Imager, INSAT-3R (ISRO); «Mereop-M»/ 
ИКФС-2, «Электро-Л»/МСУ-ГС, | «Mereop-M»/ 
МСУ-МР (Росгидромет, Роскосмос). 

4. Отмечена перспективность проведения ра- 
бот по использованию лунной калибровки аппа- 
ратуры, разработке соответствующих алгоритмов 
и расчетных моделей. 

5. Космическим агентством КНР запланиро- 
ван масштабный проект разработки национально- 
го радиометрического эталонного стенда космиче- 
ского базирования (Space-based Radiometric Bench- 
mark) с целью проведения радиометрических калиб- 
ровок ЦА от видимого до ИК-диапазона в абсолют- 
ных единицах системы СИ с высокой точностью. 

6. Современные вызовы, связанные с большим 
количеством успешных стартапов в области мно- 
госпутниковых систем на базе малых КА, суще- 
ственным образом трансформировали рынок ДЗЗ. 
Эксперты предполагают, что к 2023 г. на орбитах 
появится более 1100 новых спутников ДЗЗ. Толь- 
ко четыре компании (Planet (ранее Planet Labs), 
Maxar, Spire и BlackSky Globe) планируют за 
этот период запустить более 970 спутников [11]. 
При этом необходимо отметить, что многоспутни- 
ковые системы ДЗЗ на основе малых КА (типа 
CubeSat) на данном этапе предназначены для высо- 
кооперативного наблюдения и не решают задач точ- 
ных радиометрических измерений характеристик 
атмосферы, океана и поверхности Земли в различ- 
ных спектральных диапазонах. Однако, например, 
NASA выполняет большое количество проектов на 
основе КА типа CubeSat в рамках программы Earth 
Science Technology Program Elements (ESTO), 
которые предполагают создание ультракомпакт- 
ной ЦА, включая ИК-гиперспектрометры (HyTI), 
спектрометры солнечного излучения (СТМ), ра- 
диолокационные радиометры (SNoOPI, RainCube, 
TEMPEST-D), гиперспектральные трассовые спек- 
трометры (TACOS) и т.д. 

7. Все большее внимание уделяется спутни- 
ковым системам ДЗЗ арктической зоны, которые 


предназначены для обеспечения исследований, кон- 
троля и прогнозирования состояния атмосферы, 
океана, арктических морских льдов, прибрежной 
зоны, высоты снежного покрова, аномальных яв- 
лений в приполярных районах. 

8. Следует отметить, что в России большую 
часть работы по распространению данных ЦА гид- 
рометеорологического назначения, оценке качества 
полученных данных, проведению интеркалибровки 
и представлению научно-технических и организа- 
ционных результатов в рамках GSICS и CGMS 
проводит Росгидромет, в частности ФГБУ «НИЦ 
«Планета». Однако, эта работа невозможна без 
участия предприятий — разработчиков бортовой 
аппаратуры. Например, АО «Российские космиче- 
ские системы» является разработчиком и изгото- 
вителем такой бортовой аппаратуры, как многозо- 
нальные сканирующие устройства МСУ-МР (КА 
«Mereop-M»), МСУ-ГС (КА «Электро-Л», «Арк- 
тика-М»), МСУ-ИК-СРМ (КА «Канопус-В-ИК»), 
МСУ-МР-МП (КА «Mereop-MII»), инфракрасный 
радиометр высокого разрешения (РИВР) для мно- 
гоцелевого лабораторного модуля «Наука» МКС, 
сканеров/зондировщиков МТВЗА-ГЯ (КА «Мете- 
op-M») и МТВЗА-ГЯ-МП (KA «Mereop-MII»). 

Все вышеперечисленные приборы имеют OCO- 
бенности декодирования и обработки поступающей 
информации при штатной эксплуатации, особенно- 
сти измерительных трактов спектральных каналов, 
особенности наземной и бортовой калибровки, что 
требует полетного сопровождения данной целевой 
аппаратуры специалистами разработчика при оцен- 
ке качества данных космической съемки для ис- 
ключения систематических ошибок при проведении 
интеркалибровок. 

В докладах, представленных специалиста- 
ми АО «Российские космические системы» на 
GSICS-EP-20 и CGMS-47, авторы-разработчики 
сканеров/зондировщиков МТВЗА-ГЯ (КА серии 
«Метеор») [5], многозональных сканирующих 
устройств МСУ-ГС (КА «Электро-Л») [12] и МСУ- 
ИК-СРМ (КА «Канопус-В-ИК») [13] описывают 
принципы работы, основные технические характе- 
ристики и особенности бортовой калибровки аппа- 
ратуры с учетом оценки полученных данных кос- 
мической съемки, в том числе на основе интерка- 
либровки. 
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В частности, B докладе [12] отмечено, что 
оценка метрологических характеристик радиометра 
МСУ-ГС в составе КА «Электро-Л» №2 в сравне- 
нии с радиометрами УПВ$ (KA Suomi МРР, схема 
GEO-LEO) и SEVIRI (KA Meteosat-10, схема GEO- 
GEO), показала хорошее соответствие получаемых 
этими приборами радиометрических данных, а так- 
же возможность построения на основе данных кос- 
мических съемок МСУ-ГС тематических продуктов, 
соответствующих мировым стандартам качества. 

Результаты интеркалибровки радиометра МСУ- 
ИК-СРМ (KA «Канопус-В-ИК») [13] с исполь- 
зованием в качестве опорных приборов IASI 
(КА MetOp-A,B) и ИКФС-2 (KA «Метеор-М» №2) 
показали, что относительная разность радиометри- 
ческих данных в канале 2 МСУ-ИК-СРМ по срав- 
нению с данными [ASI (MetOp-B) и ИКФС-2 для 
спектральной плотности энергетической яркости на 
эффективной длине волны не превышает 5%. Это 
свидетельствует о хорошем качестве абсолютной 
калибровки МСУ-ИК-СРМ. 

Таким образом, интеркалибровка является хо- 
рошим инструментом для отслеживания метро- 
логических характеристик космических радиомет- 
ров при условии минимизации методических оши- 
бок, связанных с согласованием времени, широ- 
ты, долготы, высоты, солнечного зенитного угла, 
угла обзора по азимуту, зенитного угла наблюде- 
ния и спектрального диапазона при съемке одного 
объекта двумя различными приборами. 
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Аннотация. В статье рассматриваются особенности обработки данных сверхвысокого разрешения аппаратуры «Геотон-Л1» 
космического аппарата «Ресурс-П» в обеспечение формирования бесшовных сплошных покрытий (БСП) на регионы Россий- 
ской Федерации. Приведены методика формирования БСП и требования к опорной информации. Обозначены проблемные 
вопросы отбора данных и создания продуктов стандартной обработки данных «Геотон-Л1». Приведена технология формиро- 
вания бесшовных сплошных покрытий, обеспечивающая точность геодезической привязки на опорных точках не хуже 5 м. 
На примере обработки данных Белгородской области (более 27 тыс. км”) дано распределение временных затрат по техно- 
логическим процессам: 2/3 всего времени занимает обработка данных с участием оператора, при этом 1/2 всего времени — 
стандартная обработка. В рамках выполнения работ по заказу ГК «Роскосмос» с использованием рассматриваемой технологии 
сформированы БСП за летние периоды 2013-2018 гг. на 10 регионов РФ: Москва и Московская область, Санкт-Петербург 
и Ленинградская область, Севастополь и Республика Крым, Белгородская, Брянская, Воронежская, Калининградская, Курская, 
Самарская, Орловская области. 
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Abstract. The article discusses the ultra-high resolution data processing of the Geoton-Ll equipment installed on the Resurs-P 
spacecraft to ensure the generation of seamless continuous coverage (SCC) for regions ої Russia. The method for achieving SCC 
and the requirements for supporting information are given. The problematic issues of data selection and the creation of standard data 
processing Geoton-L1 products are indicated. The technology of achieving seamless continuous coverage, which ensures the accuracy 
of geodetic reference at ground control points of no worse than 5 m is described. Taking the processing of the Belgorod Region 
data (more than 27,000 square kilometers) as an example, the time distribution by technological processes is presented: 2/3 of total 
time is data processing with the operator involvement and 1/2 of total time is standard level processing. As part оѓ the efforts 
undertaken at the request of the Roscosmos State Corporation for Space Activities, using the above technology, the SCCs were 
generated for 10 regions of the Russian Federation during summer seasons ої 2013-2018: Moscow and Moscow Region, Saint 
Petersburg and Leningrad Region, Sevastopol and the Republic of Crimea, Belgorod Region, Bryansk Region, Voronezh Region, 
Kaliningrad Region, Kursk Region, Samara Region, and Orel Region. 


Keywords: Earth remote sensing, Resurs-P spacecraft, very-high resolution, seamless continuous coverage, mosaic, basic product 
bank 
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Введение 


Основные тенденции развития мировой груп- 
пировки космических аппаратов дистанционного 
зондирования Земли высокого (1-10 м) и сверхвы- 
сокого (лучше 1 м) пространственного разрешения 
[1-3] предусматривают: 1) улучшение простран- 
ственного разрешения (до 0,1-0,3 м) для новых 
спутников; 2) повышение частоты съемки за счет со- 
здания многочисленных группировок малых спутни- 
ков. Кроме того, наметился устойчивый тренд опе- 
ративного предоставления космической информации 
и информационных продуктов дистанционного зон- 
дирования Земли (ДЗЗ) на основе веб-сервисов (на- 
пример, https://eos.com/ landviewer). При этом cy- 
щественное увеличение объемов космической ин- 
формации ДЗЗ высокого разрешения определило 
еще один тренд — регламентированное формиро- 
вание и обновление мозаик или бесшовных сплош- 
ных покрытий для отдельных регионов и всего ми- 
ра (например, https://www.planet.com). 

Отечественная группировка в период 2013- 
2018 гг. насчитывала до трех действующих кос- 
мических аппаратов (КА) ДЗЗ типа «Ресурс-П», 
оборудованных сенсором «Геотон-Л1» сверхвысо- 
кого разрешения [4]. На начало 2020 г. по це- 
левому назначению функционировал только один 
КА ДЗЗ типа «Ресурс-П». Вместе с тем вос- 
требованность в данных и информационных про- 
дуктах ДЗЗ неуклонно возрастает. Более того, 
немалый интерес проявляют потребители к бес- 
шовным сплошным покрытиям больших террито- 
рий на основе данных ДЗЗ, технологии форми- 
рования которых созданы АО «Российские кос- 
мические системы» в рамках выполнения опыт- 
но-конструкторских работ Федеральной космиче- 
ской программы России на 2016-2025 гг. Примеры 
применения этих технологий для автоматическо- 
го формирования бесшовных сплошных покры- 
тий (БСП) на основе данных среднего разре- 
шения КА «Ресурс-П» и «Метеор-М» рассмот- 
рены в работах [5, 6]. Предоставление инфор- 
мационных продуктов и БСП потребителям осу- 
ществляется оператором КС ДЗЗ с использовани- 
ем геоинформационных сервисов «Банк базовых 
продуктов» (https://bbp.ntsomz.ru) [7,8] и «Геопор- 
тал Роскосмоса» (https://gptl.ru) [9,10]. В рамках 


данной статьи рассматриваются особенности обра- 
ботки данных сверхвысокого разрешения аппарату- 
ры «Геотон-Л1» космического аппарата «Ресурс-П» 
при формировании бесшовных сплошных покрытий 
на регионы Российской Федерации. 


Методика формирования 
бесшовного сплошного покрытия 


Общая методика формирования бесшовного 
сплошного покрытия на основе набора сцен дан- 
ных аппаратуры «Геотон-Л1» КА серии «Ресурс-П» 
включает следующие основные этапы. 

1. Стандартная обработка набора сцен до уров- 
ня 4À в соответствии с утвержденным Госкорпора- 
цией «Роскосмос» «Перечнем стандартных продук- 
тов, изготавливаемых по материалам ДЗЗ с КА ти- 
па “Ресурс-П”» и поставляемых потребителям [11]: 


e геометрическая и радиометрическая коррекции; 

® трансформирование на плоскость; 

e расчет параметров ВРС; 

e синтез спектральных ЁСВ-каналов и паншар- 
пенинг на основе данных панхроматического 
спектрального канала. 

2. Фотограмметрическая обработка набора сцен: 


ө поиск И задание связующих и опорных точек; 


e блочное уравнивание набора сцен на основе 
множества связующих и опорных точек; 


e ортотрансформирование набора сцен. 
3. Построение БСП: 


e построение и редактирование линий порезов; 
e расчет параметров яркостного выравнивания; 


e формирование единого результирующего изоб- 
ражения БСП. 


Определим требования к условиям съемки 
и опорной информации в обеспечение формирова- 


ния БСП. Во-первых, с учетом пространственно- 
го разрешения (0,7 м) аппаратуры «Геотон-Л1» КА 
«Ресурс-П» необходимая точность планового поло- 
жения опорных точек должна составлять не менее 
половины пикселя — 0,35 м. Во-вторых, на рис. 1 
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Рис. 1. Зависимость точности задания высоты наблюда- 

емой точки местности на снимке «Геотон-Л1» КА серии 

«Ресурс-П» от угла крена КА, обеспечивающая смеще- 
ние образа точки на снимке в пределах 0,5 пикселя 


приведена зависимость ошибки задания геодези- 
ческой высоты наблюдаемой точки местности от 
углов крена съемочной аппаратуры «Геотон-Л1» 
КА «Ресурс-П», при которых наблюдается смеще- 


ние образа точки местности на снимке в преде- 
лах половины пикселя. Соответственно для «на- 
дирного» снимка точность задания высоты долж- 
на быть лучше 10 м, при углах крена менее 30° — 
не хуже 1 м. Таким образом, для формирования 
БСП с точностью (среднеквадратическим отклоне- 
нием (CKO)) до 3-4 м необходимо использовать: 
1) снимки с углами визирования/крена не более 
10°; 2) опорные точки с точностью 0,35 м в плане 
и по высоте не хуже 2 м; 3) цифровую модель ре- 
льефа (ЦМР) с максимальной ошибкой не более 
12 м (3 СКО). 

Далее следует отметить проблему выбора ис- 
ходных данных аппаратуры «Геотон-Л1» KA «Pe- 
сурс-П», связанную с отсутствием планового 
сплошного покрытия территории РФ. На рис. 2 
приведено покрытие Ленинградской области марш- 
рутами (сценами) съемки за летние периоды 
2013-2018 гг. с уровнем облачности менее 50%. 


Рис. 2. Набор маршрутов/сцен съемки аппаратуры «Геотон-Л1» КА серии «Ресурс-П» территории Ленинградской 
области за летние периоды 2013-2018 гг. 
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Фактически отснято 2/3 региона, однако ввиду 
отсутствия перекрытия между сценами (с уров- 
нем облачности менее 30%) представляется воз- 
можным формирование сплошного покрытия толь- 
ко для 1/3 региона (граница обозначена желтым 
цветом на рис. 2). Для увеличения площади бес- 
шовного покрытия необходимо использовать сце- 
ны, отснятые в условиях визирования более 20-30? 
по крену и/или с уровнем облачности 50-100%, 
что не всегда целесообразно, поскольку от этого 
зависит качество результирующего изображения. 


Проблемные вопросы стандартной 
обработки данных «Teoron-JI1» КА 
серии «Ресурс-П» в обеспечение 
формирования БСП 


На рис. 3 представлена схема расположе- 
ния ПЗС-матриц в фокальной плоскости сенсора 
«Геотон-Л1» КА серии «Ресурс-П». Как видно из 
схемы, спектральные каналы Ве4 (0,61-0,68 мкм), 
Green (0,52-0,60 мкм) находятся в рамках одной 
группы ПЗС-матриц, канал Вше (0,45-0,52 мкм) — 
в рамках другой группы ПЗС-матриц. Каждая из 
этих групп равноудалена от группы ПЗС-матриц 
панхроматического канала (0,58-0,80 мкм), нахо- 
дящейся в плоскости оптического центра. 

Приведенная оптическая схема расположения 
различных групп ПЗС-матриц обусловливает про- 


блему ВСВ-синтеза спектральных каналов, а имен- 
но разновременную съемку одной и той же точ- 
ки местности разными группами ПЗС-матриц. Как 
следствие, попиксельное сведение спектральных 
каналов возможно за счет использования высоты 
наблюдаемой точки местности (на рис. 4 приве- 
дена траектория a,—Pj-b,). При этом ошибки зада- 
ния высоты наблюдаемой точки местности приводят 
к ошибкам синтеза спектральных каналов (на рис. 4 
приведена траектория a,—P5-b,). Ha рис. 5 при- 
веден фрагмент изображения, на котором демон- 
стрируются ошибки синтеза спектральных кана- 
лов, связанных с неточностью задания высоты на- 
блюдаемого высотного объекта. 

Решение проблем синтеза спектральных ка- 
налов возможно в случае реализации следующих 
конструктивных и технологических задач. 

1. Применение светоразделительных оптиче- 
ских систем, обеспечивающих одномоментную 
съемку в различных спектральных каналах (дан- 
ный прием реализован для аппаратуры ГСА КА 
«Ресурс-П»). 

2. Минимизация расстояния между группами 
ПЗС-матриц таким образом, чтобы временной ин- 
тервал задержки наблюдения одной и той же точки 
местности был согласован с пространственным раз- 
решением аппаратуры (данное решение реализовано 
для спектральных каналов Red и Green в рамках од- 
ной группы ПЗС-матриц аппаратуры «Геотон-Л1»). 

3. Согласование пространственного разреше- 
ния аппаратуры с цифровой моделью рельефа 
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Рис. З. Схема расположения ПЗС-матриц B фокальной плоскости сенсора «Геотон-ЛІ» КА серии «Ресурс-П». 
Справа приведены параметры светофильтров ПЗС 
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Рис. 4. Схема синтеза спектральных каналов при ошибках задания высоты наблюдаемого объекта 


Рис. 5. Ошибки синтеза спектральных каналов данных «Геотон-Л1» КА серии «Ресурс-П» для высотных объектов 


(ЦМР), используемой при синтезе спектральных 
каналов (данное требование характерно для ап- 
паратуры среднего пространственного разрешения 
КА «Ресурс-П» — КШМСА [9]). 

4. Синтез изображений различных спектраль- 
ных каналов на основе нелинейных деформа- 
ций изображений (данное решение реализовано 
в НКПОР KA «Ресурс-П» в виде специального 
программного обеспечения, требующего участия 
оператора для контроля триангуляционной сети, 
см. рис. 6). 


Технология формирования БСП 


Общая технология формирования БСП, по дан- 
ным «Геотон-Л1» КА серии «Ресурс-П», приведе- 
на на рис. 7. Технология включает этап формиро- 
вания продуктов стандартной обработки, ориентиро- 
ванный на программные решения НКПОР, в рамках 
которых выполняется создание продуктов уровня 4А 
по «Перечню стандартных продуктов, изготавлива- 
емых по материалам ДЗЗ с КА типа “Ресурс-П”» 
[11] B автоматизированном режиме (с участием 
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Рис. 6. Пример синтеза спектральных каналов данных «Геотон-Л1» КА серии «Ресурс-П» с использованием KOP- 
реляционного сопоставления и построения триангуляционной сети 
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Рис. 7. Технология формирования БСП по данным аппаратуры «Геотон-Л1» КА типа «Ресурс-П» 


оператора). Далее для каждого маршрута выполня- а также фрагментирование (см. рис. 8) маршрута 
ется автоматическое уточнение геопривязки (на ос- на сцены. Обработанные маршруты/сцены аккуму- 
нове опорных данных и цифровой модели рельефа), лируются в обменном буфере, а также загружаются 
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Рис. 8. Примеры фрагментированных маршрутов аппаратуры «Геотон-Л1» КА серии «Ресурс-П» на территорию 
Белгородской области 


в геоинформационный сервис «Банк базовых про- 
дуктов», где производится их каталогизация и ар- 
хивирование. 

В рамках данной работы для автоматическо- 
го измерения плановых координат опорных точек 
использовалось покрытие, предоставляемое Planet 
Labs (пространственное разрешение — З м), точ- 
ность геодезической привязки которого не хуже 
10м [12]. Доступ к этому покрытию реализован 
в рамках библиотеки GDAL (https://gdal.org). B ка- 
честве ЦМР применялась SRTM (https://doi.org/ 
10.5066/Е7РВ7ТЕТ). 

Далее выполнялось блочное уравнивание для 
заданного набора маршрутов/сцен и последующее 


их ортотрансформирование. Например, 
ном уравнивании данных «Геотон-Л1» КА серии 
«Ресурс-П» на территорию Республики Крым (пе- 
репад высот более 1500 м) СКО на опорных точ- 
ках составило не более 5 м [13]. При этом для pe- 
гионов с перепадом высот в несколько сотен мет- 
ров, например Белгородской области, СКО привяз- 
ки отдельных сцен составляло около 1,5 пикселей 
опоры или не более 5 м (см. рис. 9). Соответствен- 
но для таких «равнинных» регионов формирова- 
ние БСП осуществлялось без блочного уравнива- 
ния и для сведения контуров ортотрасформирован- 
ных маршрутов/сцен, при необходимости применя- 
лись методы локальной деформации, в частности 


при блоч- 
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Рис. 9. Распределение а) систематических и б) случайных ошибок (м) геопривязки сцен «Геотон-Л1» КА серии 
«Ресурс-П», используемых при формировании БСП на территорию Белгородской области 


Рис. 10. Пример БСП территории Белгородской области, сформированного по данным аппаратуры «Геотон-Л1» 
КА серии «Ресурс-П» за летние периоды 2013-2018 гг. 
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Таблица. Временные оценки формирования БСП Бел- 
городской области 


Технологический этап 


Время, ч | Объем, ТБ 


Стандартная обработка 
Фрагментирование 


Уточнение привязки 


Формирование результирую- 
щего изображения 


* — этап с участием оператора. ** — данные со сжа- 


тием. 


реализованные в фотограмметрическом пакете 
Photomod. Результирующее изображение БСП Бел- 
городской области приведено на рис. 10. 

В таблице показаны временные оценки техно- 
логических этапов формирования покрытия (вклю- 
чая накладные расходы, связанные с копированием 
данных) на примере формирования БСП Белгород- 
ской области на основе обработки более 30 марш- 
рутов данных «Геотон-Л1» КА серии «Ресурс-П». 
Соответствующая таблице диаграмма распределе- 
ния временных затрат показана на рис. 11. Анализ 
диаграммы показывает, что работа оператора со- 
ставляет 2/3 всего времени формирования конеч- 
ного продукта, при этом более половины време- 
ни уходит на стандартную обработку данных «Гео- 
тон-Л1» КА «Ресурс-П». 


m Стандартная обработка 
B Фрагментирование 

и Уточнение привязки 

B Ортотрансформирование 
E Редактирование порезов 


B Формирование 
результирующего 
изображения 


Рис. 11. Диаграмма распределения времени между тех- 
нологическими этапами формирования БСП Белгород- 
ской области 


Следует отметить, что этапы автоматической 
обработки данных выполнялись на двух вычис- 
лительных узлах, организованных в соответствии 
с моделью линейного конвейера. Узел фрагменти- 
рования в том числе являлся буфером стандарт- 
ных продуктов НКПОР. Для уточнения геопривяз- 
ки и ортотрансформирования применялся вычисли- 
тельный узел, обеспечивающий также копирование 
ортотрансформированных сцен на вычислительный 
узел оператора формирования БСП. С учетом меж- 
сетевого взаимодействия, а также накладных рас- 
ходов одновременного чтения/записи обменных бу- 
феров обеспечивались номинальные показатели ав- 
томатической обработки данных — 15-25 Мб/с. 
Для увеличения скорости автоматической обработ- 
ки планируется дальнейшая оптимизация техноло- 
гической схемы. 


Заключение 


В рамках выполнения контрактных работ по 
заказу Госкорпорации «Роскосмос» сформированы 
бесшовные сплошные покрытия по данным аппара- 
туры «Геотон-Л1» КА серии «Ресурс-П» за летние 
периоды 2013-2018 гг. на следующие регионы РФ 
(см. рис. 10, 12): Москва и Московская область, 
Санкт-Петербург и Ленинградская область, Сева- 
стополь и Республика Крым, Белгородская, Брян- 
ская, Воронежская, Калининградская, Курская, Са- 
марская, Орловская области. 

В общей сложности на основе данных стандарт- 
ных продуктов ДЗЗ объемом около 10 ТБ (c уче- 
том сжатия) сформировано более 330 тыс. км? 
(более 3,5 ТБ) бесшовных сплошных покрытий ре- 
гионов Российской Федерации сверхвысокого про- 
странственного разрешения. 
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Аннотация. В сетях подвижной связи стандарта GSM на 51М-карте абонента хранятся постоянные ключ аутентификации 
и алгоритм аутентификации. Для аутентификации абонента используются постоянные ключ аутентификации, алгоритм аутен- 
тификации и алгоритм шифрования, которые задействуются для аутентификации абонента и шифрования передаваемой инфор- 
ации. В случае компрометации ключа аутентификации алгоритм аутентификации и алгоритм шифрования злоумышленника 
смогут входить в связь с базовой станцией под чужим идентификатором телефона IMSI (International Mobile Subscriber Iden- 
tity). Представлена технология усиления криптостойкости алгоритмов аутентификации и шифрования информации, применяе- 
ых в сетях подвижной связи GSM, в основе которой вместо постоянных ключа аутентификации, алгоритма аутентификации 
алгоритма шифрования лежит применение динамично изменяющихся матриц ключей и алгоритмов. Предлагаемое техниче- 
кое решение может иметь широкое применение в системах радиосвязи, сетях подвижной связи и персональной подвижной 
путниковой связи, а также позволит избирательно подходить к дополнительному усилению криптостойкости передаваемой 
нформации для отдельных категорий абонентов. Положительным эффектом также является то, что повышение криптостойко- 
сти сети связи осуществляется без значительного увеличения вычислительной нагрузки на процессор абонентского терминала, 
что обеспечивает экономию заряда аккумуляторной батареи абонентского терминала. 
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Abstract. Іп GSM mobile networks, а subscriber's SIM card stores a permanent authentication key and an authentication algo- 
rithm. For authentication of the subscriber is used a permanent authentication key, an authentication algorithm, as well as an 
encryption algorithm, that used to encrypt the transmitted information. If the authentication key is compromised, an attacker will 
be able to communicate with the base station using someone else's IMSI (International Mobile Subscriber Identity). The technology 
lor enhancing the cryptographic strength of authentication algorithms and information encryption used in GSM mobile networks 
is presented, which, instead of using a permanent authentication key, authentication algorithm and encryption algorithm, is based 
on the use of dynamically changing key matrices and algorithms. The proposed technical solution can be widely used in radio 
communication systems, mobile networks, and personal mobile satellite communications. It will also make it possible to differ- 
entiate additional strengthening of the cryptographic stability of the transmitted information for certain categories of subscribers. 
An additional positive effect is that the increase in the cryptographic strength of the communication network is carried out without 
a significant increase in the computational load on the processor ої the subscriber terminal, which saves the battery power of the 
subscriber terminal. 
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Введение 


По оценке аналитиков, уровень потерь опера- 
торов мобильной связи от разного рода мошенниче- 
ства и вредительства составляет от 2-6% до 25% 
от общего объема трафика. Причем атаки мошенни- 
ков направлены как против операторов, так и про- 
тив абонентов [1]. Решение проблемы обеспечения 
безопасности в российских сетях связи осложняет- 
ся широким использованием технических средств 
зарубежного производства, что создает возмож- 
ность реализации не декларируемых поставщиками 
возможностей [1]. 

Для передачи информации цифровой телефо- 
нии и коротких сообщений (SMS) с использованием 
персональных абонентских терминалов наибольшее 
распространение получил стандарт GSM. В сетях 
подвижной связи (СПС) GSM для аутентификации 
абонента используются постоянные ключ аутенти- 
фикации К; и алгоритм аутентификации A;, а для 
шифрования информации применяются постоян- 
ный алгоритм шифрования Ag, которые привязаны 
к SIM-xapre абонента [2]. Злоумышленник, KOTO- 
рый единожды вскрыл ключ аутентификации K;, 
алгоритм аутентификации А; и алгоритм шифро- 
вания Ág, может постоянно расшифровывать HH- 
формацию в абонентском канале связи с конкрет- 
ным пользователем СПС GSM. Отдельные катего- 
рии абонентов, например представители крупного 
бизнеса [3,4] или госслужащие [5], нуждаются 
в повышении криптостойкости услуг связи СПС 
GSM [6] и сетях подвижной персональной спут- 
никовой связи [7]. 


Алгоритмы аутентификации 
и шифрования в СПС GSM 


В современных сетях подвижной (мобильной) 
связи применяется алгоритм шифрования, основан- 
ный на следующей технологии [2]. 

Для аутентификации информации абонента 
в сетях подвижной связи стандарта GSM каждый 
абонентский терминал (АТ) содержит записанные 
в 51М-карте свои индивидуальный ключ аутенти- 
фикации К; и алгоритм аутентификации A;. baso- 
вая станция (БС) сети подвижной связи передает 


на АТ случайное число RAND, а АТ на основании 
числа RAND, ключа аутентификации K; и алго- 
ритма аутентификации А, вычисляет значение OT- 
вета SRES: ЗВЕ$ = A;(K;; RAND). Далее AT or- 
правляет ответ SRES на базовую станцию, KOTO- 
рая использует ключ К, и алгоритм аутентифика- 
ции А, вызываемого абонента, хранящиеся в базах 
данных зарегистрированных абонентов мобильного 
центра коммутации, вычисляет число SRES и срав- 
нивает вычисленное число SRES с полученным от 
абонентского терминала. Если вычисленное и полу- 
ченное от АТ число SRES совпадают, то БС всту- 
пает в связь с АТ. 

То есть в данном случае используется аутенти- 
фикация с открытым ключом, характеризующаяся 
следующими недостатками: 

— ключ К; и алгоритм аутентификации А; no- 
стоянные, следствием чего является низкая крип- 
тостойкость; 

— случайное число RAND передается по от- 
крытому каналу связи и при компрометации ключа 
и алгоритма аутентификации злоумышленник мо- 
жет осуществить ложное подключение к базовой 
станции (т.е. создать ложный АТ, осуществив под- 
лог), пока на АТ, ключи и алгоритмы аутентифи- 
кации которого скомпрометированы, не будет сме- 
нена SIM-xapra; 

— отсутствует возможность увеличения крип- 
тостойкости аутентификации для отдельных ка- 
тегорий АТ без применения дополнительных про- 
граммных средств. 

Для шифрования информации абонентской 
радиолинии в сети подвижной связи стандар- 
та GSM БС передает Ha АТ случайное число 
RAND. Абонентский терминал на основании чис- 
ла RAND, ключа К; и алгоритма шифрования Аз 
вычисляет ключ шифрования на сеанс связи К’: 
К. = AQ(K;RAND). Передаваемая в абонент- 
ской радиолинии (от базовой станции и от або- 
нентского терминала) информация зашифровыва- 
ется с применением ключа шифрования на се- 
анс связи К. и алгоритма шифрования Ag: Е, = 
= Аз(Ко; Ie), где I, — открытая информация в ce- 
ансе связи. 

То есть используют шифрование с открытым 
ключом, характеризующееся следующими недо- 
статками: 
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— алгоритм вычисления А; ключа аутентифи- 
кации /;, из которого формируется ключ шифро- 
вания информации в сеансе связи А’, постоянный, 
поэтому злоумышленник имеет возможность ском- 
прометировать ключ аутентификации J€;, а затем 
научится формировать ключи шифрования инфор- 
мации в сеансе связи Ke; 


- алгоритм шифрования Ág, используемый 
для формирования ключа шифрования информации 
в сеансе связи А, из ключа аутентификации К, 
и применяемый для шифрования информации, по- 
стоянный, поэтому данный способ шифрования об- 
ладает низкой криптостойкостью; 


— случайное число RAND передается по от- 
крытому каналу связи и при компрометации ключа 
аутентификации K;, алгоритма аутентификации A; 
и алгоритма шифрования А, злоумышленник мо- 
жет прослушивать трафик АТ, пока не будет сме- 
нена 51М-карта; 


- для отдельных категорий АТ отсутствует 
возможность увеличения криптостойкости шифро- 
вания без применения дополнительных програм- 
мных средств. 


Метод аутентификации абонента 
с применением динамично 
изменяющихся матриц ключей 

и алгоритмов 


В предлагаемом методе аутентификации [8] 
данные аутентификации и шифрования передава- 
емой информации задействуются при установлении 
каждого соединения и передаче данных и органи- 
зованы в памяти каждого составляющего систе- 
му устройства как матрица ключей аутентифика- 
ции Mg; матрица порядков смены позиций Мск; 
ключей аутентификации К; в матрице ключей 
аутентификации Mg; матрица алгоритмов аутен- 
тификации Mj; матрица порядков смены пози- 
ций Мод; алгоритмов аутентификации А, в матри- 
це M 4;; матрица М дз алгоритмов шифрования Аз; 
матрица порядков смены позиций Мү. дз алгорит- 
мов шифрования Ag в матрице Мл алгоритмов 
шифрования. 


При работе предложенной системы передачи 
данных периодические изменения позиций ключей 
в матрице ключей, а также изменение алгорит- 
мов в матрице алгоритмов являются альтернати- 
вой хранения в постоянном запоминающем устрой- 
стве большого количества ключей и алгоритмов. 
Работа системы основана на двух последователь- 
ных процессах: аутентификации информации або- 
нентской радиолинии и шифрования абонентской 
радиолинии. Порядок смены ключей аутентифика- 
ции А,, формулы математических операций и Mexa- 
низм смены алгоритмов аутентификации А, и алго- 
ритмов шифрования А; известны каждому центру 
коммутации системы — наземной станции сети по- 
движной связи и/или спутнику-ретранслятору сети 
персональной спутниковой связи. Процесс аутен- 
тификации удостоверяет, что абонент имеет право 
доступа к услугам связи и предшествует процедуре 
установления соединения. После установления со- 
единения начинает работу процедура шифрования 
информации в абонентской радиолинии. 

Пример простейшей матрицы ключей аутенти- 
фикации — 


Fig, Kio R 


где В — это принятое от БС число RAND. 
Пример простейшей матрицы смены позиций 
ключей аутентификации — 


Скар Criz: Сказ 
Мск: = Срд C icio C ki23 
Сказ Cki Сказ 


Пример простейшей матрицы алгоритмов ау- 
тентификации — 


Аз 1 Ai» Ais 
Ma = Ар А99 А әз 
Азу iso ias 


Пример простейшей матрицы смены позиций 
алгоритмов аутентификации — 


С лиг Cai» Caig 
Moai = |CAi Cai» Caiz 
Ci Cai Caisa 
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Пример простейшей матрицы алгоритмов шиф- 
рования — 


Agi 1 Аз1о Аз 
М дз = Agi Ago» gos 
831 gap Аззз 


Пример простейшей матрицы смены позиций 


алгоритмов == 
Саву C дв? С Ag13 
Moae = С A821 C A829 C Ag23 
С A831 C A832 C A833 


Аутентификация AT выполняется следующим 
образом. 

Шаг 1. Матрицы My. ключей аутентифика- 
ции, матрицы порядков смены позиций Мск; клю- 
чей аутентификации, матрицы алгоритмов аутен- 
тификации Мд.,, матрицы порядков смены пози- 
ций Mg, алгоритмов аутентификации, матрицы 
алгоритмов шифрования M 4g и матрицы порядков 
смены позиций M 4g алгоритмов шифрования за- 
писываются в постоянное энергонезависимое запо- 
минающее устройство (ПЗУ) каждого АТ (или на 
SIM-kapry), в защищенные базы данных центров 
коммутации и являются уникальными для каж- 
дого АТ. 

Шаг 2. Принятое АТ от центра коммутации 
случайное число RAND записывается в одну из по- 
зиций матрицы Мк; ключей аутентификации для 
проведения процедуры вычисления числа SRES. 
Позиция матрицы Мк; ключей аутентификации, 
в которую записывается принятое случайное число 
RAND, определяется порядком смены позиций C ge; 
ключа аутентификации. 

Шаг 3. Значения ключей аутентификации К, 
и принятое случайное число ВАМО, записанных 
в матрице Mg; ключей аутентификации, периоди- 
чески перемещаются по позициям матрицы клю- 
чей аутентификации по определенному алгоритму 
в соответствии с действующим на данный момент 
порядком смены позиций Ск; ключей аутентифи- 
кации. Порядковые номера порядков смены пози- 
ций Ср; ключа аутентификации записаны в MAT- 
рице Мск, порядков смены позиций ключа аутен- 
тификации, и каждому номеру порядков смены по- 
зиций Ск; ключа аутентификации соответствует 
определенная последовательность смены позиций. 


После прохождения полного цикла смены по- 
зиций ключами аутентификации K; и принятым 
случайным числом RAND ключи аутентифика- 
ции А, и принятое случайное число RAND Bos- 
вращаются на исходные позиции и изменяется по- 
рядок смены их позиций Cg; После прохождения 
полного цикла изменения порядков смены пози- 
ций Скд; ключами аутентификации порядки сме- 
ны позиций Ск д; возвращаются на исходные пози- 
ции в матрице M opgi порядков смены позиций и по- 
рядок смены позиций Ск д; изменяется по опреде- 
ленному закону. 

Шаг 4. Номера алгоритмов А, вычисления чис- 
ла SRES (алгоритмов аутентификации) записыва- 
ются в матрицу аутентификации алгоритмов M ,;, 
и каждому порядковому номеру соответствует своя 
математическая формула. 

Шаг 5. Номера алгоритмов аутентификации Aj, 
записанных в матрице M4; алгоритмов аутенти- 
фикации, перемещаются по позициям матрицы ал- 
горитмов аутентификации по определенному пра- 
вилу в соответствии с действующим на данный 
момент порядком смены позиций C; алгорит- 
мов аутентификации. Вначале номера алгоритмов 
аутентификации А; изменяются по одному прави- 
лу Сан действующему в настоящий момент, за- 
тем порядок смены позиций алгоритма аутенти- 
фикации C; изменяется. Номера порядков смены 
позиций C4; алгоритма аутентификации записаны 
в матрице Мод; порядков смены позиций алгорит- 
ма аутентификации, и каждому номеру порядков 
смены позиций C 4; алгоритма аутентификации CO- 
ответствует определенная последовательность сме- 
ны позиций. 

После прохождения полного цикла смены по- 
зиций номерами алгоритмов аутентификации А, 
номера алгоритмов аутентификации возвращаются 
на исходные позиции и изменяется порядок смены 
их позиций C ,;. После прохождения полного цик- 
ла изменения порядков смены позиций C ,; номера- 
ми алгоритмов аутентификации А, порядки смены 
позиций С’; возвращаются на исходные позиции 
в матрице Meg; порядков смены позиций и поря- 
док смены позиций C ,; изменяется по определен- 
ному закону. 

Смена позиций ключами аутентификации Кд; 
и номеров алгоритмов аутентификации А; осущест- 
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вляется через строго определенные интервалы 
времени (от 1 суток до полугода), известные АТ 
и центрам коммутации. Все алгоритмы вычисле- 
ний A; ответа SRES содержат однонаправленные 
функции. 

Шаг 6. В AT на основании матрицы M g; клю- 
чей аутентификации, принятого случайного числа 
RAND и матрицы алгоритмов аутентификации Мл; 
производятся математические операции для вычис- 
ления числа SRES. 

Шаг 7. Вычисленное число SRES передается 
в центр коммутации, где также вычисляется число 
SRES. Вычисленное число SRES сравнивается с по- 
лученным от АТ числом SRES. Абонент считается 
прошедшим аутентификацию при совпадении чисел 
SRES, вычисленных АТ и центром коммутации. 

Криптостойкость процесса аутентификации 
определяется: 

— длиной ключа аутентификации Ly; (для уве- 
личения криптостойкости длина ключа шифрова- 
ния Lg; должна стремится к максимально возмож- 
ному значению); 

— количеством применяемых ключей аутенти- 
фикации Мк; и количеством перестановок клю- 
чей Мск; (для увеличения криптостойкости KOJIH- 
чество применяемых ключей аутентификации и ко- 
личество перестановок ключей должны стремиться 
к максимально возможным значениям); 

— сложностью и количеством алгоритмов ау- 
тентификации IN 4; (для увеличения криптостойко- 
сти сложность и количество алгоритмов аутенти- 
фикации должны стремиться к максимально воз- 
можным значениям); 

— количеством перестановок алгоритмов аутен- 
тификации No 4; (для увеличения криптостойкости 
количество перестановок алгоритмов аутентифика- 
ции должно стремиться к максимально возможным 
значениям); 

— периодом действия ключей аутентифика- 
ции Тк; и алгоритмов аутентификации T; (для 
увеличения криптостойкости период действия клю- 
чей аутентификации и алгоритмов аутентификации 
должны стремиться к минимально возможным зна- 
чениям). 

Алгоритм работы АТ и центра коммутации 
СПС или СР сети ПСС при выполнении процедуры 
аутентификации изображен на рис. 1 и 2. 


Метод шифрования информации 
с применением динамично 
изменяющихся матриц ключей 

и алгоритмов 


В предлагаемом методе [8] шифрование ин- 
формации выполняется следующим образом. 

Шаг 1. Ключ шифрования K, на сеанс связи 
определяется из вычисленного на этапе аутенти- 
фикации ключа аутентификации К; с применени- 
ем действующего в настоящий момент алгоритма 
шифрования As. 

Шаг 2. Номера алгоритмов шифрования А; 
записываются в матрицу алгоритмов шифрова- 
ния Mg. Каждому порядковому номеру соответ- 
ствует заданная математическая формула. Номера 
алгоритмов шифрования Аз, записанных B матри- 
це алгоритмов шифрования M ,„, перемещаются по 
позициям данной матрицы по определенному пра- 
вилу в соответствии с действующим на данный 
момент порядком смены позиций C 4g алгоритмов 
шифрования. 

Вначале номера алгоритмов шифрования А; 
изменяются по одному правилу C 4g, действующе- 
му в настоящий момент, затем порядок смены по- 
зиций алгоритма шифрования C 4, изменяется. Ho- 
мера порядков смены позиций С’; алгоритма шиф- 
рования записаны в матрице Mo 4s порядков смены 
позиций алгоритма шифрования, и каждому номе- 
ру порядков смены позиций C 4g алгоритма шифро- 
вания соответствует определенная последователь- 
ность смены позиций. После прохождения полного 
цикла смены позиций номерами алгоритмов шиф- 
рования А, номера алгоритмов шифрования воз- 
вращаются на исходные позиции и изменяется по- 
рядок смены их позиций C дз. 

После прохождения полного цикла изменения 
порядков смены позиций C 4g номерами алгоритмов 
шифрования А; порядки смены позиций C ,g BO3- 
вращаются на исходные позиции в матрице Mo дв 
порядков смены позиций и порядок смены пози- 
ций C 4g изменяется по определенному закону. 

Смена позиций номеров алгоритмов шифрова- 
ния Ag осуществляется через строго определенные 
интервалы времени (от 1 суток до полугода), из- 
вестные АТ и центрам коммутации. Все алгоритмы 
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Начало работы процедуры 
аутентификации АТ 


| 


Получение абонентом 
случайного числа ВАМО 


О ЗИ 
| | 


Определение текущего состояния Определение текущего алгоритма 
матрицы ключей аутентификации аутентификации 


Время действия 
текущего состояния матрицы 
ключей аутентификации 
истекло? 


Время действия 
текущего алгоритма 
аутентификации 
истекло? 


Определение текущего порядка Определение текущего порядка 
смены позиций матрицы ключей смены алгоритмов 
аутентификации аутентификации 


Текущий порядок 
смены позиций алгоритмов 
аутентификации прошел 
полный цикл? 


Текущий порядок 
смены позиций ключей 
аутентификации прошел 
полный цикл? 


Выбор очередного порядка смены 


» » Вы чередного порядка смены 
позиций матрицы ключей бор очередного порядка сме 


алгоритмов аутентификации 


аутентификации 
Изменение позиций матрицы Выбор очередного алгоритма 
ключей аутентификации по правилу аутентификации по правилу 
очередного порядка смены позиций порядка смены позиций 


| l 
m 


Процедура вычисления 
числа SRES 


| 


Отправка числа SRES на БС 


Рис. 1. Алгоритм работы АТ при выполнении процедуры аутентификации 
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| 


| 


Числа SRES, 
вычисленное и полученное 
от абонента, 
совпадают? 


Нет 


Y 
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Аутентификация абонента 
не выполнена 


Y 
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Рис. 2. Алгоритм работы центра коммутации сети подвижной связи или спутника-ретранслятора сети подвижной 
персональной спутниковой связи при выполнении процедуры аутентификации 
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вычислений Ag ключа шифрования К. содержат 
однонаправленные функции. 

Шаг 3. Информация зашифровывается с помо- 
щью ключа шифрования К’ на сеанс связи и алго- 
ритма шифрования Аз. 

Криптостойкость процесса шифрования ин- 
формации определяется: 

- длиной ключа аутентификации Ly; (для 
увеличения криптостойкости длина ключа шифро- 
вания Lg; должна стремиться к максимально воз- 
можному значению); 

— количеством применяемых ключей аутенти- 
фикации Мк; и количеством перестановок клю- 
чей Nog, (для увеличения криптостойкости KOJIH- 
чество применяемых ключей аутентификации и ко- 
личество перестановок ключей должны стремиться 
к максимально возможным значениям); 

— сложностью и количеством алгоритмов аутен- 
тификации №; и алгоритмов шифрования Ns 
(для увеличения криптостойкости сложность и ко- 
личество алгоритмов аутентификации и алгорит- 
мов шифрования должны стремиться к максимально 
возможным значениям); 

— количеством перестановок алгоритмов аутен- 
тификации No; и алгоритмов шифрования No дв 
(для увеличения криптостойкости количество пе- 
рестановок алгоритмов аутентификации должно 
стремиться к максимально возможным значениям); 

— периодом действия ключей аутентифика- 
ции Тк,, алгоритмов аутентификации T; и an- 
горитмов шифрования T's (для увеличения крип- 
тостойкости период действия ключей аутентифи- 
кации, алгоритмов аутентификации и алгоритмов 
шифрования должен стремиться к минимально воз- 
можным значениям). 

Алгоритм процедуры шифрования в сети по- 
движной связи или спутника-ретранслятора в сети 
подвижной персональной спутниковой связи изоб- 
ражен на рис. 3. 

Дешифрование информации выполняется сле- 
дующим образом. 

Шаг 1. В центре коммутации так же, как 
и в АТ, на основании вычисленного на этапе аутен- 
тификации ключа аутентификации К; и матрицы 
алгоритмов шифрования M agi производятся мате- 
матические операции для вычисления ключа шиф- 
рования на сеанс связи К... 


Начало работы процедуры 
шифрования информации 


| 


Определение текущего алгоритма 
шифрования 


Время действия 
текущего алгоритма 
шифрования 
истекло? 


Нет 


Определение текущего порядка 
смены позиций алгоритмов 
шифрования 


Текущий порядок 
смены позиций алгоритмов 
шифрования прошел 
полный цикл? 


Выбор очередного порядка смены 
позиций алгоритмов шифрования 


3 


Выбор очередного алгоритма 
шифрования по правилу 
очередного порядка смены позиций 


| 


Процедура вычисления ключа 
шифрования на сеанс связи 
на основании ключа аутентификации, 
принятого числа RAND и алгоритма 
шифрования: К, = А, (К; RAND) 


| 


Процедура шифрования 
информации на основании 
ключа шифрования и алгоритма 
шифрования: E, = А,(К,; [.) 


Рис. 3. Алгоритм процедуры шифрования в сети подвиж- 
ной связи или спутника-ретранслятора в сети подвиж- 


ной персональной спутниковой связи 
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Шаг 2. Затем в центре коммутации на основа- 
нии вычисленного ключа шифрования на сеанс свя- 
зи К, и действующего в настоящий момент алго- 
ритма шифрования А; выполняется дешифрование 
информации (функция, обратная шифрованию). 


Методика управления степенью 
криптостойкости передаваемой 
информации 


Для управления степенью криптостойкости 
передаваемой информации и, соответственно, отка- 
за от усложнения системы в ряде случаев абонен- 
тов классифицируют по восьми категориям от 0-й 
до 7-й [8]: 

1) абоненты с минимальными требованиями 
к криптостойкости — 0-й, 1-й, 2-й категорий — 
«массовый» пользователь (потенциальный наруши- 
тель — хакер, работающий единолично на обычных 
бытовых компьютерах); 

2) абоненты с повышенными требованиями 
к криптостойкости — 3-й, 4-й категорий — высоко- 
поставленные представители бизнеса и госслужа- 
щие (потенциальный нарушитель — организован- 
ная группа хакеров, работающая в интересах пре- 
ступных организаций или конкурентной разведки 
в небольших центрах обработки данных); 

3) абоненты с высокими требованиями к крип- 
тостойкости — 5-й, 6-й, 7-й категорий — высокопо- 
ставленные госслужащие и представители силовых 
структур и ведомств (потенциальный нарушитель — 
инженеры, имеющие в распоряжении крупные цен- 
тры обработки данных). 

В зависимости от категории абонентов приме- 
няют: 

— разные по размеру матрицы ключей аутен- 
тификации, матрицы порядков смены позиций клю- 
чей аутентификации, матрицы алгоритмов аутенти- 
фикации и матрицы порядков смены позиций алго- 
ритмов аутентификации; 

— разные по размеру матрицы алгоритмов 
шифрования и матрицы порядков смены позиций 
алгоритмов шифрования; 

— разные временные сроки действия опреде- 
ленной комбинации ключей аутентификации; 


— разные временные сроки действия опреде- 
ленного алгоритма аутентификации и определен- 
ного алгоритма шифрования; 

— разные по сложности алгоритмы аутентифи- 
кации и алгоритмы шифрования; 

— разную длительность интервала времени 
действия определенной комбинации ключей аутен- 
тификации, алгоритмов аутентификации и алгорит- 
мов шифрования, которая должна быть меньше пе- 
риода, необходимого для компрометации ключей 
и алгоритмов. 

Например, для категории абонентов с макси- 
мальными требованиями к криптостойкости инфор- 
мации ключи можно записывать не в одну, а в две 
матрицы и производить операцию умножения од- 
ной матрицы на другую. 


Заключение 


В рассмотренных методах аутентификации 
и шифрования информации за счет одновременно- 
го обеспечения на всех устройствах большого ко- 
личества вариантов числа аутентификации и клю- 
ча шифрования информации абонентского канала 
связи обеспечивается повышение криптостойкости 
системы передачи данных с минимальной допол- 
нительной нагрузкой на вычислительные мощно- 
сти АТ. 

Таким образом, разработана и обоснована ме- 
тодика аутентификации абонентов и шифрования 
информации, обладающая следующими преимуще- 
ствами: 

— простотой реализации, не требующей высо- 
кой производительности вычислительных средств; 

— универсальностью применения в различных 
системах радиосвязи и спутниковой связи, в том 
числе в СПС и в СППСС военного и двойного на- 
значения; 

— гибким подходом к изменению степени крип- 
тостойкости для различных категорий трафика 
и различных категорий абонентов. 
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Аннотация. Исследованы процессы гальваномагнитного транспорта в тонких пленках мягкого ферромагнетного сплава 
Мі, 6 ео 18 Соо ов, применяемого в анизотропных магниторезистивных датчиках (АМР) индукции магнитного поля космическо- 
го применения. Описаны результаты измерения зависимостей магнитофизических и магнитоэлектрических параметров тонкой 
пленки Nio 6 ео 18 Сооов 6-зондовым методом. Обнаружена величина Ap/p > 5,0%. Проведен анализ процесса токопереноса 
в сплаве в рамках модели трехмерной топологии поверхностей Ферми Лифшица-Каганова. Сделан вывод, что наблюдаемое по- 
вышение Ap/p — следствие эффекта разрыва дырочной зоны проводимости из-за присутствия атомов кобальта в кубической 
подрешетке железа. Результаты могут быть применены для повышения чувствительности АМР-сенсоров и для разработки 
новых материалов спинтроники. 
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Abstract. Galvanomagnetic transfer processes in thin films of the soft magnetic alloy № в Ее 18 Содоо used in space-grade 
anisotropic magnetoresistive (AMR) magnetic field's sensors are studied. The results of measuring the magnetophysical and 
magnetoelectric parameters of a thin film of Nigz;FegsCogos alloy using a 6-probe method are described. The value Ap/p > 5.0% 
was detected. The analysis of the current transfer process in the alloy is carried out according to the model the three-dimensional 
topology of the Fermi surfaces proposed by Lifshitz-Kaganov. It is concluded that the increase in Ар/р is a consequence ої the 
rupture of the hole conduction band due to the presence of cobalt atoms in the cubic iron sublattice. Obtained results can be used 
Гог upgrade high-sensitive AMR-sensors and for new spintronic materials development. 
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Введение 


Происходящее в настоящий момент одновре- 
менное развитие наноспутниковых технологий и пе- 
реход к высокоинформативным методам многоплат- 
форменных исследований геомагнитного поля ост- 
ро поднимает вопрос о миниатюризации магни- 
тометрической аппаратуры космических аппаратов 
(ММА КА) [1]. В основу принципов действия чув- 
ствительных элементов (ЧЭ) ММА КА могут быть 
положены более десятка разных физических эф- 
фектов, но наилучшим потенциалом миниатюриза- 
ции при высокой чувствительности, малом энерго- 
потреблении и высокой стойкости к внешним фак- 
торам обладают анизотропные магниторезистивные 
сенсоры (АМР-сенсоры) [1,2]. АМР-сенсоры пред- 
ставляют собой специальные резисторы, сделанные 
из тонкой пленки мягкого ферромагнетика, напы- 
ленной на диэлектрическую подложку. Миниатю- 
ризация АМР-сенсоров легко достигается за счет 
их изготовления методами микроэлектронной тех- 
нологии (напыление тонких пленок, фотолитогра- 
фия, травление). На основе АМР-эффекта в по- 
следние годы был разработан ряд миниатюрных 
ЧЭ магнитной индукции, используемых в бортовых 
векторных системах ММА КА [3-5]. 

Однако чувствительность и разрешающая спо- 
собность современных АМР-магнитометров при- 
мерно на 1-2 порядка хуже, чем у наиболее ши- 
роко применяемого типа векторных магнитомет- 
ров ММА КА — флюксгейт-магнитометра. Из-за 
требований, связанных с повышением отноше- 
ния сигнал-шум и сохранением температурной 
стабильности, наиболее часто используемым для 
АМР-сенсоров материалом является пермаллой co- 
става №0 Ред јә, имеющий Ap/p ~ 1,8% [2,6]. 
Повышение чувствительности АМР-сенсора мо- 
жет быть достигнуто совершенствованием тополо- 
гии [4,5] либо повышением величины магнитосо- 
противления Ар/р исходной пленки. 

Второй путь представлялся более сложным, 
хотя еще McGuire и Potter [2] предложили два 
сплава — NiggjFegog и Соо 7F€o з, демонстрирую- 
щие в объемном образце Ар/р > 5,0%. Однако 
данных о тонкопленочных АМР-сенсорах на их 
основе до сих пор нет. Tumanski [6, с. 99] отме- 
чал, что «...обоснованное прогнозирование соот- 


ветствия между стехиометрическим составом фер- 
ромагнитного сплава, его остаточным удельным со- 
противлением ри Ар/р представляет собой зада- 
чу, не решенную даже на теоретическом уровне». 
Тем не менее ранее были получены тонкие плен- 
ки сплава № gFeg 18 Соо [7-9]. Эти пленки изо- 
морфны Nigg;Feg;o, но имеют меньшую величи- 
ну р и воспроизводимо демонстрируют значение 
Др/р > 5,0%. Из них были изготовлены ЧЭ, при- 
годные для применения в ММА КА [3,4]. В дан- 
ной работе выполнено исследование гальваномаг- 
нитных свойств тонких пленок Nip zgFeo 18 C00 og; 
описываются результаты эксперимента и анализи- 
руется возможный механизм влияния добавки Со 
на механизм токопереноса в пленке. 


Образцы и методика измерений 


В качестве опытных образцов на ситалло- 
вых подложках размером 60 x 48 мм были Ha- 
пылены 5 однослойных ферромагнитных пленок 
Nip 76Рео 1300 og. Напыление производилось в про- 
мышленной установке вакуумного напыления типа 
YPM3-279, оборудованной ионно-плазменным HC- 
точником и магнитной системой, обеспечивающей 
наведение в плоскости подложки магнитного по- 
ля напряженностью 200 Гс. Напыление производи- 
лось с мишени типа H2K18K6, содержащей соот- 
ветственно 76% Ni, 18% Fe и 6% Co. 

Толщина слоя Nig zgFeo 13 Соо og Контролирова- 


лась по времени напыления и составила 250 А. 
Для измерения магнитофизических параметров 
пленок (R, Н, H,, dR/R) использовался стенд 
ТИДН.411182.001, изготовленный в СКТБ ИРЭ 
РАН (рис. 1) [8], который является первой рос- 
сийской установкой, реализующей 6-зондовый ме- 
тод измерения магнитосопротивления, предложен- 
ный J. Cavallo [10]. Так же как и в традиционном 
4-зондовом, в 6-зондовом способе измерения сигна- 
лом является разность потенциалов У между дву- 
мя точками. В схему измерений 6-зондового ме- 
тода введена дополнительная пара потенциальных 
контактов (контакты 5 и 6 на рис. 1,6). То есть 
одна пара контактов (3 и 4) располагается в по- 
зиции максимума сигнала от магнитосопротивле- 
ния, а другая пара (5 и 6) — в позиции минимума 
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Рис. 1. a) Общий вид измерительного стенда ТИДН.411182.001, разработанного СКТБ ИРЭ РАН: / — катушки 

Гельмгольца, 2 — поворотный столик, 3 — образец, 4 — блок электроники и управления, 5 — 6-зондовая изме- 

рительная головка, 6 — привод измерительной головки, 7 — привод поворотного столика; б) схема расположения 
измерительных контактов зондовой головки 


сигнала. Линия расположения контактов под 45° 
относительно оси легкого намагничивания (ОЛН) 
пленки обеспечивает максимальный сигнал У от 
магнитосопротивления при прохождении положе- 
ния In(1 + Ар/р). При этом У между любой парой 
измерительных контактов будет описываться урав- 
нением 


у = IR E (1 + A Coos + sing?) = 
27 Po 


—]n (1 + AP (edd = sng?) (1) 
Po 


где ру — удельное сопротивление, перпендикуляр- 
ное M,; 

Др — магнитосопротивление; 

9 — угол между направлением намагниченно- 
сти и осью т. 

Шестизондовым методом непосредственно на 
поверхности ферромагнитной пленки можно изме- 
рить локальное значение поверхностного сопротив- 
ления R, (Ом/квадрат) и вычислить из него коэф- 
фициент магнитосопротивления dR/R;: 


В, max _ В, тіп 


dR[R, = Taxes (2) 


ѕ max 

Магнитное поле напряженностью +150 Э co- 
здается парой катушек Гельмгольца (позиция [ на 
рис. 1, а). На поворотный столик (2) устанавлива- 
ются исследуемые образцы (3). Блок электроники 


и управления 4 подает калиброванный ток на токо- 
вые контакты (1 и 2) 6-зондовой измерительной TO- 
ловки (5) и измеряет разность потенциалов на двух 
парах потенциальных контактов. Также блок элек- 
троники и управления управляет приводом измери- 
тельной головки (6) и приводом поворотного сто- 
лика (7). Измерительный стенд ТИДН.411182.001 
позволяет снимать ом-эрстедные характеристики 
пленок при различных углах поворота столика, 
а также в различных точках пластины. Коэрцитив- 
ная сила H, определялась как разность АН поло- 
жения пиков R max прямой и обратной ветви за- 
висимости А, = f(H). Положение ретроградного 
участка на нисходящей ветви этого графика соот- 
ветствует величине поля смещения H,. При этом 
поле магнитной анизотропии Н, представляет со- 
бой полуширину основания пика. 

По данным семейств ом-эрстедных характери- 
стик при разных углах поворота 9 для каждого об- 
разца были построены зависимости dR/R, = }(0), 


Н. = f(0), Н, = f(0). 


Результаты измерений 


Для всех образцов были сняты семейства 
ом-эрстедных характеристик в диапазоне углов по- 
ворота 0-360° с шагом 15°. За нулевое направление 
принималось технологически заданное направление 
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Др/р = 5,15% в пленке № „Ее уа Coo ов (образец Ш, 
таблица) 


ОЛН. Вид ом-эрстедной характеристики образца Ш 
при 0 = 0,0? приведен на рис. 2. 

По полученным характеристикам были вычис- 
лены основные магнитофизические параметры пле- 
нок и их изменение в зависимости от угла поворота 
(рис. 3). 
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Рис. 3. Изменение А, B в магнитном поле напряжен- 
ностью в диапазоне +30 Э при повороте столика от 0° 
до 90° (образец Ш, таблица) 


В таблице приведены все полученные при из- 
мерении магнитофизические параметры образцов. 

Величины Н., H, и dR/R указаны для на- 
правления, параллельного ОЛН (0 = 0,0?). Более 
наглядным и информативным является представ- 
ление величин dR/R, Н. и Н, в виде полярных 
диаграмм их зависимостей от угла поворота под- 
ложки. В этом случае обретает физический смысл 
собственная осевая симметрия описанной кривой, 


А. B. ГРАБОВ 


Таблица. 


Магнитофизические параметры образцов 
(отклонение от ОЛН — 0°) 


E EEA EEA 
m 15,2 |14,8 |14,4 | 
[s 


ее симметрия относительно полярной точки, CHM- 
метрия относительно ОЛН (0-180°) и оси трудного 
намагничивания (ОТН) (90?-270?). 

Диаграмма зависимости dR/R = f(0) подло- 
жек приведена на рис. 4. 


fée | 


=] 
-- II 
—Ш 
= IV 
э ү 


150 
165 


195 180 


Рис. 4. Диаграмма dR/R = }(0) 


Кривые имеют гантелеобразную форму с ожи- 
даемой перетяжкой в области ОТН, где d R/ R всех 
образцов падает практически до нуля. Площадь, 
занимаемая каждой кривой, четко соответствует 
максимальной величине магниторезистивного эф- 
фекта в пленке. 

В диаграмме основное внимание следует обра- 
тить на направление осей собственной симметрии 
кривых, которое характеризует отклонение направ- 
ления максимального магниторезистивного эффек- 
та от ОЛН. Приведенная диаграмма демонстрирует, 
что подавляющее большинство образцов имеет су- 
щественное отклонение максимума dR/R от ОЛН. 
Величина этого отклонения заключена, как прави- 
ло, в диапазоне от 2°—3° до 15°. Величина отклоне- 
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Рис. 5. Диаграмма H, = { (0). + - обр. I, A — обр. П, 
Ш — обр. Ш, e — обр. IV, = — обр. V 


ния не связана ни с величиной аб / В, ни с величи- 
ной ағ. Так, наименьшее отклонение имеют образ- 
цы Ш и IV, которые имеют соответственно наиболь- 
шую и наименьшую абсолютную величину d R/ R из 
исследованной выборки, а наибольшее отклонение 
дают образцы Ги П, которые в выборке имеют наи- 
большее и наименьшее значение ao. 

Наиболее характерным в диаграмме зависимо- 
сти H, = f(0) ферромагнитных пленок от угла 
(рис. 5) является четко выраженный минимум H, 
в области ОЛН. Причем понятие «минимум» дохо- 
дит вплоть до того, что в ряде образцов в области 
этого минимума наблюдается отрицательная вели- 
чина коэрцитивности, что на первый взгляд ка- 
жется лишенным какоголибо физического смысла. 
Кажущееся отрицательное значение H, в сплаве 
Nio;gFeg;&Cogog наблюдалось впервые и требует 
дополнительного анализа. Форма кривых Н, = 
= f(0) тоже может быть охарактеризована как ган- 
телеобразная (рис. 6). Ось симметрии «перетяжки» 
направлена вдоль ОТН. 


Обсуждение результатов 


С точки зрения современной теории магнито- 
сопротивление является гальваномагнитным транс- 
портным феноменом [11]. Перенос заряда про- 
исходит в ферромагнитном материале, состоящем 
из локализованных магнитных моментов, имеющих 


195 


180 165 


Рис. 6. Диаграмма Н, = f(0). + - обр. І, A — обр. П, 
Ш — обр. Ш, e — обр. IV, = — обр. V 


трансляционную симметрию. Магнитное поле за- 
кручивает носители заряда, а также поляризует их 
спины. Поляризация фиксирует направление спи- 
нов, и исчезает возможность рассеяния с поворотом 
спина. Магнитосопротивление в такой модели явля- 
ется следствием снижения полной вероятности рас- 
сеяния [12]. 

Использовавшийся В работе сплав 
Nip 760,18 C00 og представляет собой типичный 
твердый раствор замещения [13], составленный из 
трех различных 34-металлов. Он образован мат- 
рицей гранецентрированной кубической (ГЦК) ре- 
шеткой Ni с параметром ячейки а ~ 3,524 А 
и растворенными в ней атомами Ее и Со. Если ре- 
шетка немодифицированного М Feo; состоит из 
простой кубической решетки Ее и вставленной в нее 
октаэдрической решетки Ni, то B № 76 Ее 18 Соо о6 
в одной четверти узлов кубической подрешетки Со 
частично замещает Fe. При этом структура само- 
го Со испытывает искажение вследствие наруше- 
ния симметрии. Вид кристаллической ячейки сплава 
Nip 76 ео 18 Соб, ов приведен на рис. 7. 

Токоперенос в сплаве Мі, Ее Сок ocy- 
ществляется как по электронному, так и по дыроч- 
ному механизму. Для случая токопереноса двумя 
типами носителей известно следующее выражение 
для величины Ар/р [14]: 

Ap: — 000 — Роз)? B? (3) 
Po (а +05)? + (Rio? + Во)? В? 
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Рис. 7. ячейка сплава 
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где 0, и с» — собственные проводимости для элек- 
тронов и дырок, В; и Ё, — постоянные Холла для 
электронов и дырок соответственно, В — напря- 
женность магнитного поля. 

Следует обратить внимание, что в общем слу- 
чае в магнитном поле для носителей с двумя спи- 
новыми состояниями будет происходить расщепле- 
ние по эффективной массе m* [15]. Следовательно, 
процессы токопереноса в ферромагнитном сплаве 
будут характеризоваться 4 значениями Op — двумя 
электронными и двумя дырочными. Значение Op 
для каждого типа носителей заряда можно вычис- 
лить из известной формулы, связывающей длину 
свободного пробега данного типа носителей — lp 
и интеграл площади поверхности Ферми — 5 [16]: 


E 
с, = mx; | 148 
12121, 

S 


(4) 


Для анализа конкретного сплава B TOM чис- 
ле необходимо проанализировать форму и тополо- 
гию поверхностей Ферми каждого компонента. От- 
метим, что до широкого внедрения методов ней- 
тронного рассеяния основным способом исследова- 
ния поверхности Ферми в металлах и полупровод- 
никах было именно изучение зависимости Ар/р = 
= f(0) [14]. 

Анализировались поверхности Ферми Ni, Fe 
и Со для моделей, разработанных Вакохой и Яма- 
шитой (ВЯ) [17]. На рис. 8 приведена качествен- 
ная модель энергетических зон для одной ячей- 
ки сплава Nig;gFeg;gCogog, построенная в mph- 


А.Б. ГРАБОВ 


Рис. 8. Модель топологии поверхностей Ферми в кри- 
сталлической ячейке сплава № ео у Cog og 


ближении «жестких зон», т.е. исходя из пред- 
положения, что атомы элементов в сплаве в це- 
лом сохраняют форму поверхности Ферми, свой- 
ственную для элементарных веществ. После транс- 
ляции этой ячейки получается трехмерная сет- 
ка из ферми-поверхностей, анализируемая метода- 
ми анализа трехмерных топологий, разработанных 
И. М. Лифшицем и М. И. Кагановым (ЛК) [18, 19] 
и развитыми в работах В. Г. Песчанского [20], 
А.Я. Мальцева [12,21] и других [22]. 

Теория ЛК давно заняла место в арсенале ин- 
струментов теоретической физики металлов [16]. 
Однако в прикладных работах по АМР-сенсорам 
она используется редко, поэтому целесообразно 
привести основные положения ЛК, важные для ис- 
следования гальваномагнитных явлений в ферро- 
магнетиках: 

1) носители заряда могут перемещаться в кри- 
сталлах только по топологически разрешенным 
траекториям, обращаясь по ним с ларморовской ча- 
СТОТОЙ W; 

2) в зависимости от геометрии трехмерной 
сетки из ферми-поверхностей разрешенные траек- 
тории для носителей заряда могут быть как откры- 
тыми, так и замкнутыми. Открытые орбиты харак- 
теризуются зависимостью вида Лр/р о В?, a sa- 
крытые — Ap/p ~ B}; 
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3) разрешенные траектории могут иметь то- 
пологию двух типов: «гофрированный цилиндр» 
и «трехмерная решетка из трубок». Если поверх- 
ность Ферми цилиндрической топологии, в ней име- 
ется открытая траектория в том случае, если вектор 
индукции магнитного поля лежит в плоскости пер- 
пендикулярной оси цилиндра. Если плоскость на- 
правлена под углом — траектория замкнутая [18]; 

4) если магнитное поле наклонено к плоско- 
сти, нормальной к оси «гофрированного цилиндра» 
под малым углом 9, то длина траектории носителя 
заряда и период его обращения возрастает пропор- 
ционально 0-!. 

Проанализируем по ЛК транслированную 
ячейку сплава № 6 Ее 1 Соо о6с формой энергети- 
ческих зон, соответствующей моделям ВЯ в пер- 
вом приближении. В ГЦК-матрице Ni электроны 
спином «вниз» локализованы в двух больших за- 
мкнутых электронных карманах с центрами в точ- 
ке Г в 5-й и 6-й зонах [14, с. 229]. Поскольку 
карманы замкнутые, то, согласно [18], электроны, 
находящиеся в них, также имеют замкнутую тра- 
екторию и не участвуют в процессе гальваномаг- 
нитного токопереноса. Для электронов со спином 
«вверх» имеется многосвязанный лист поверхности 
Ферми с перешейками вдоль направлений (111). Пе- 
рекрываясь с соседними атомами, эти листы образу- 
ют топологию типа «трехмерная решетка», которая 
и должна играть главную роль в гальваномагнит- 
ном транспорте в никелевой матрице. Дырочные 
карманы у Ni располагаются в точке Г, 4-й зоны 
и B точке X 3-й и 4-й зонах и трехмерной тополо- 
гически связанной решетки не образуют [17]. 

Модель ВЯ энергетических зон ОЦК Ее топо- 
логически гораздо сложнее никелевой. Электрон- 
ные поверхности Ферми состоят из большой ок- 
таэдрической замкнутой поверхности в точке Г 
4-й зоны со структурой типа «валет», а также 
группы малых электронных линз, расположенных 
вдоль Г-Н и образующих перешейки между «ва- 
летами». Другими словами, в железе токопере- 
нос может осуществляться электронами и спином 
«вверх», и спином «вниз». Еще более сложную 
топологию имеют дырочные поверхности Ее. Они 
включают не только два дырочных кармана в точ- 
ке Н, но и сложную многосвязанную трубчатую 
дырочную поверхность, открытую вдоль направле- 
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ний (001), (111), (110). Перекрываясь с соседни- 
ми атомами, эта трубчатая решетка Ее, в отли- 
чие от Ni, создает открытую траекторию для ды- 
рок. Эффективная масса дырок больше электрон- 
ной, поэтому движение носителей по этой траекто- 
рии будет давать вклад в Н.. 

Зонная структура Со В сплаве 
№іо,лЕео18С0оо6в будет характеризоваться боль- 
шей неопределенностью [14, с. 232]. Чистый Со 
гексагонален, и модель ВЯ предусматривает для 
Со наличие 18 гексагонально-симметричных по- 
верхностей Ферми. Большую часть объема ячейки 
занимает сдвоенная электронная сфера в точке Г, 
три концентрические дырочные поверхности, 
а также три группы мелких электронных и ды- 
рочных карманов. При этом по ВЯ в Со больших 
дырочных зон, способных образовать структуру 
«трехмерной решетки», нет. Согласно модели ЛК, 
наличие зависимости вида Ap/pcoB? доказы- 
вает наличие открытых орбит в направлении, 
нормальном к В. Полученная гантелеобразная за- 
висимость Ар/р = f(0) (рис. 4) достаточно хорошо 
аппроксимируется известными уравнениями сину- 
соидальной спирали 2-го порядка или лемнискаты 
Бернулли в полярных координатах [14, 19]: 

Др 


—^ = B? . cos 20, 


А (5) 


которое говорит, что у сплава Nig Ее 1 С006 мы 
имеем дело с поверхностью Ферми, открытой вдоль 
ОЛН и закрытой вдоль ОТН для электронов во 
всех исследованных образцах (относительно тре- 
буемого для построения лемнискаты квадрата па- 
раметра Ар/р см. формулу (3)). 

Аналогичная аппроксимация Н. уравнением 
лемнискаты типа 


H2 = B? .соѕ(—920) (6) 


возможна только для образца IV (рис. 5), имею- 
щего наиболее низкое значение Ap/p = 3,45%. 
То есть этот образец отличается от осталь- 
ных наличием открытых дырочных траекторий 
вдоль ОТН. Наличие ретроградного участка на 
ом-эрстедной характеристике № в Еес у Соо ов 06- 
разца ІУ указывает на существование области 
преимущественно дырочного токопереноса. Между 
тем мы видим, что в образце Ш ретроградный 
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участок отсутствует, a H, при всех значениях 0 
не превышает 1,8 Э, причем форма Н, = f(0) 
не гантелеобразная, а «звездчатая», соответству- 
ющая зависимости Н, со В!, что свидетельствует 
о закрытости в образце III дырочных траекторий 
во всех направлениях и меньшей роли дырочного 
механизма в общем токопереносе. 

В связи с этим логично будет предположить, 
что, поскольку Со не имеет значительных дыроч- 
ных зон, атомы Со, внедренные в угловые по- 
зиции в подрешетке железа в структуре сплава 
№0 7gFe 1800,06 так, как это показано на рис. 7, 
как бы «разрывают» открытые Ее дырочные тра- 
ектории, не меняя формы траекторий электрон- 
ных. Вследствие этого вклад дырочного тока 
в суммарный перенос заряда в модифицирован- 
HOM Ni 7gFeo 15 C09 os существенно меньше, чем 
в немодифицированном М з Рес sg. Соответствен- 
но меньше и вклад «тяжелых» дырок из Ее в ко- 
эрцитивность пленки. А учитывая то, что Со да- 
ет дополнительный электрон в никелевую матри- 
цу, имеющую открытые электронные траектории, 
можно утверждать, что, по сравнению с немоди- 
фицированным пермаллоем, о, в № Ее 13Со0 os, 
наоборот, будет выше, что и наблюдалось 
(р. (№0 76Еео 18С0006) = 6,2-9,0 Ом/квадрат, 
р. (№0 g; Feo 19) = 22,0 Ом/квадрат [6]). 


Выводы 


На основании изложенного можно сформули- 
ровать следующие выводы: 

— повышение чувствительности АМР-сенсо- 
ров путем повышения Ар/р исходной пленки воз- 
можно, в частности, при использовании тройного 
сплава Nig 76 Рес Соо gg, имеющего Ap/p = 5,0%; 

— кристаллографический изоморфизм структу- 
ры Nio 26 Ред 18С0оов И Nio g1F€o 19 < высокой Bepo- 
ятностью выражается в замене атомами Со одной 
четверти атомов в подрешетке железа. Замена Fe 
на Со приводит к «разрыву» открытой спаренны- 
ми атомами Fe в М№ з ео 19 дырочной траектории, 
что снижает коэрцитивность пленки и улучшает ее 
гистерезисные свойства. 

Анализ взаимосвязи свойств гальваномагнит- 
ного транспорта с кристаллографической структу- 


А.Б. ГРАБОВ 


рой сплавов методами теории ЛК позволяет отно- 
сительно просто и наглядно моделировать свойства 
новых материалов магнитоэлектроники и спинтро- 
ники, что может оказаться полезным для разработ- 
ки новых высокочувствительных и малошумящих 
миниатюрных АМР-сенсоров, в том числе и для 
ММА КА. 
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